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' The momentum bands, energy dispersions, and velocities of the charge c fermions and spin-neutral 

I | two-spinon si fermions of a square-lattice quantum liquid referring to the Hubbard model on such 

q . a lattice of edge length L in the one- and two-electron subspace are studied. The model involves 

■ the effective nearest-neighbor integral t and on-site repulsion U and can be experimentally realized 
in systems of correlated ultra-cold fermionic atoms on an optical lattice and thus our results are 

, of interest for such systems. Our investigations profit from a general rotated-electron description, 

■ which is consistent with the model global SO(3) x 5*0(3) x U(l) symmetry. For the model in the one- 
' and two-electron subspace the discrete momentum values of the c and si fermions are good quantum 
, numbers so that in contrast to the original strongly-correlated electronic problem their interactions 

^v^j ■ are residual. The use of our description renders an involved many-electron problem into a quantum 

liquid with some similarities with a Fermi liquid. For the Hubbard model on a square lattice in the 
i i one- and two-electron subspace a composite si fermion consists of a spin-singlet spinon pair plus an 

infinitely thin flux tube attached to it. In the U/At — > oo limit of infinite on-site interaction the c 
I 1 fermions become non-interacting spinless fermions and the si fermion occupancy configurations that 

| generate the spin degrees of freedom of spin-density m = ground states become within a suitable 

CZ2 . mean-field approximation for the fictitious magnetic field B s \ e X3 brought about by the correlations 

■ of the original N electron problem those of a full lowest Landau level with N /2 degenerate one-sl 
fermion states of the two-dimensional quantum Hall effect. In turn, for U/At finite the degeneracy of 
the N/2 one-sl-fermion states is removed by the emergence of a finite-energy-bandwidth si fermion 

I 1 dispersion yet the number of si band discrete momentum values remains being given by B sl L 2 /*o 

and the si effective lattice spacing by a s i = l s i/V2n where hi is the fictitious-magnetic-field length 
and in our units the fictitious-magnetic-field flux quantum reads $o — 1- Elsewhere it is found 
that the use of the square-lattice quantum liquid of charge c fermions and spin-neutral two-spinon 
si fermions investigated here contributes to the further understanding of the role of electronic 
correlations in the unusual properties of the hole-doped cuprate superconductors. This indicates 
£f) . that quantum-Hall-type behavior with or without magnetic field may be ubiquitous in nature. 

l/-) ' PACS numbers: 71.10.Fd, 71.10.Pm, 71.27.a, 75.30.Ds 

oo : 

m ■ 

(N ; I. INTRODUCTION 

^r" ; 

The Hubbard model on a square lattice is the simplest realistic toy model for description of the electronic correlation 
effects in general many-electron problems with short-range interaction on such a lattice and therefore is the obvious 
starting point for the study of the role of such effects in the exotic physics of the hole-doped cupratesi~— . The model 
involves two effective parameters: the in-plane nearest-neighbor transfer integral t and the effective on-site repulsion 
U. Despite that it is among the mostly studied models in condensed matter physics, there is no exact solution and 
few controlled approximations exist for finite U /At values. 

In this paper we study the momentum bands, energy dispersions, and velocities of the charge c fermions and spin- 
neutral two-spinon si fermions introduced for the model on the square lattice in Ref* 9 -. The si fermions emerge from 
the spin- neutral two-spinon si bond particle o 9 ' 10 through an extended Jordan- Wigner transformation^— . Our study 
has as starting point the properties of the Hubbard model on the square lattice in the one- and two-electron subspace 
defined in Ref>i and profits from the general rotated-electron description introduced in that reference. 

There is consensus about the scientific interest of the Hubbard model on the square lattice as simplest toy model for 
describing the effects of electronic correlations in the high-T c superconductors^— and their Mott-Hubbard insulators 
parent compounds 15 ' 16 . However, many open questions about its properties remain unsolved. Interestingly, the model 
can be experimentally realized with unprecedented precision in systems of correlated ultra-cold fermionic atoms on 
an optical lattice^ and one may expect very detailed experimental results over a wide range of parameters to be 
available. 

The square-lattice quantum liquid of charge c fermions and spin- neutral two-spinon si fermions and the related 
general rotated-electron description of Refi£ are consistent with the global 50(3) x 5*0(3) x U(l) = [SO (A) x /7(1)]/Z 2 
symmetry found recently in ReL— for the model on any bipartite lattice. Such a global symmetry is an extension of the 
50(4) symmetry known to occur for the model on such lattice o 19 ' 20 . The extended global symmetry is related to the 
rotated electrons, which for on-site repulsion U > emerge from the the electrons through a unitary transformation 
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of the type considered in Refill, and to the local symmetries and unitary transformations studied in Refi^. 

The building blocks of the general description introduced in Ref£ are the ?7-spin-l/2 77-spinons, spin-1/2 spinons, 
and spinless and 77-spinless charge c fermions whose occupancy configurations generate the state representations of the 
?7-spin SU(2) symmetry, spin SU(2) symmetry, and U(l) symmetry, respectively, associated with the model global 
£0(3) x 50(3) x U(l) = [SU(2) x SU(2) x U(l)]/Z$ symmetry. Such three basic objects are well defined for U > 0. 
These state representations are found in Ref.— to correspond to a complete set of momentum eigenstates. The ??-spin- 
1/2 ry-spinons describe the 77-spin degrees of freedom of the rotated-electron occupancy configurations that generate 
such states involving doubly occupied and unoccupied sites, the spin-1/2 spinons the spin degrees of freedom of the 
rotated-electron configurations of the singly occupy sites, and the c fermions the charge excitations associated with 
rotated-electron motion that conserves the numbers of singly occupied, doubly occupied, and unoccupied sites. 

For U /At 3> 1 the c fermion, spinon, and 77-spinon operators become the quasicharge, spin, and pseudospin oper- 
ators, respectively, obtained from the transformation considered in Ref.- 3 , which does not introduce Hilbert-space 
constraints. Such quasicharge, spin, and pseudospin operators have for U /At 3> 1 the same expressions in terms of 
creation and annihilation electron operators as the c fermion, spinon, and 77-spinon operators for U > in terms of 
rotated-electron creation and annihilation operators. The unitary character of the transformation between electrons 
and rotated electrons then assures that the transformation that generates the c fermion, spinon, and 77-spinon opera- 
tors from the electronic operators also does not introduce Hilbert-space constraints. The vacuum of the theory is given 
in Eq. (| Al|) of Appendix A, where some basic information on the square-lattice quantum liquid one- and two-electron 
subspace defined in Refi2 is provided. There is a vacuum of general form given in that equation for each subspace 
with a constant finite value 2S C for the number of rotated electrons that singly occupy sites. S c is also the eigenvalue 
of the generator of the hidden charge U(l) global symmetry^. For both the model on the square and one-dimensional 
(ID) lattice such a S c > vacuum is invariant under the electron - rotated-electron unitary transformation. 

For the square-lattice quantum liquid associated with the Hubbard model in the one- and two-electron subspace 
only the charge c fermions and spin- neutral two-spinon si bond particles play an active rol o 9 ' 10 . In contrast to 
previous descriptions involving Jordan- Wigner transformation o 13 ' 14 or slave-boson representation al 24 ' 25 and referring 
to the model for large values of the on-site repulsion U or Heisenberg and t — J models, whose spinless fermions 
arise from spin-1/2 objects, the si fermions emerge from hard-core spin-neutral two-spinon composite objects and 
are well defined for finite U values of the on-site repulsion. As mentioned above, the spin-1/2 spinons are the spins of 
the rotated electrons that singly occupied lattice sites and emerge from a suitably electron - rotated-electron unitary 
transformation introduced in Ref* 9 -. Therefore, here the single-occupancy constraint is naturally fulfilled. 

For the Hubbard model on the square lattice in the one- and two-electron subspace as defined in that reference 
there are no olv fermions other than none or one zero-momentum s2 fermion and [S c — S s ] or [S c — S s — 2] si fermions, 
respectively, where S s denotes the spin. The interest of the description introduced in Ref* 9 lies in the fact that for the 
subspaces of the one- and two-electron subspace spanned by mutually neutral states where the si fermion operators act 
onto, the two si translation generators q s \ Xl and q s ix 2 m the presence of the fictitious magnetic field B s \ associated 
with the si fermion Jordan- Wigner transformation commute with each other and with both the Hamiltonian and 
momentum operator. In contrast, the Hubbard model on the square lattice in the whole Hilbert space does not 
commute with such translation generators 9 -. Therefore, for the square-lattice quantum liquid corresponding to the 
Hubbard model on a square lattice in the one- and two-electron subspace the si fermion discrete momentum values 
ijj = [qj Xl , Qj X2 ] are good quantum numbers. Here qj and qj are eigenvalues of the two si translation generators 
q s \ Xl and q s ix 2 , respectively. In addition, the description of Refi^ has been constructed to inherently the c fermion 
discrete momentum values being good quantum numbers for the whole Hilbert space. In turn, our method involves 
approximations to derive the shape of the si momentum band boundary and the form of the c and si fermion energy 
dispersions. The fulfillment of such tasks is the main goal of this paper. For the ID model the discrete momentum 
values of the c and si fermions in units of 2it / L are quantum numbers of the exact solution 9 -. For the model on the 
square lattice the nearest-neighboring components of the c and si discrete momentum values qj = [Qj Xl >Qj X2 ] are 
such that [qj — qj' .] — 2n/L where i — 1 or i — 2 and the indices j and j' refer to nearest neighboring discrete 
momentum values. The shape of the c band is that of the first Brillouin zone. Hence the main open issue is the 
shape of the si boundary line for x > and m = ground states, and that along with the c and si fermion energy 
dispersions is one of the problems studied in this paper. 

It follows from the results of Ref j 9 - that for the Mott- Hubbard insulator at hole concentration x = and spin density 
m = both the c and si bands are full. For the Hubbard model on a square-lattice the c band has for all energy 
eigenstates the same momentum area and shape as the electronic first Brillouin zone. For m = ground states with 
a small finite hole concentration x there arises a circular c Fermi line around —7? = — [ir,ir] of ratio ss y/xir 2, which 
encloses a c fermion unfilled momentum area xAtt 2 . Since the c effective lattice equals the original lattice and thus 
its spacing a c equals the original lattice constant a, the c fermion occupancy configurations that generate the energy 
eigenstates conserve translational invariance. 

In turn, consistently with the ground-state si effective-lattice occupancies of Ref. 9 , for x > and m — ground 
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states the si band is full. For the x = and m = absolute ground state that band has a momentum area 2tt 2 and 
is found in this paper to coincide with an antiferromagnetic reduced Brillouin zone such that \q Xl \ + |9x 2 | ^ which 
is enclosed by a boundary whose momenta qssi belong to the lines connecting [±7r, 0] and [0, ±7r]. Such a si band 
shape and momentum area are consistent with the result of Ref.— that at x — and m = the square si effective 
lattice has spacing a s i = \[2a. Indeed our studies confirm that its periodicity has increased relative to that of the 
original lattice owing to the appearance of a long-range antiferromagnetic order. Moreover, we find that the Fermi 
line has for the x = and m = ground state the same form as the si boundary line so that it refers to the lines 
connecting [±7r,0] and [0, ±7r]. For very small values of x, both the momentum areas enclosed by the si boundary 
line and the Fermi line decrease to (1 — x)2tt 2 . Then the latter line becomes hole like and both such lines remain near 
the lines connecting [±7r, 0] and [0, ±tt]. This is consistent with the recent experimental results of Ref^, which reveal 
that for small hole concentrations the momentum-space region near such lines plays a major role in the Fermi-line 
physics of the hole-doped cuprate superconductors. 

Upon further increasing the hole concentration x the si momentum band remains full for x > and m = 
ground states and encloses a smaller momentum area (1 — x)2ir 2 , alike for small x. For < x < x* and m = 
there is a short-range spin order. Here .x* G (0.23,0.28) for approximately uq < U/At < u\ where uq ~ 1-3 and 
iti ~ 1.6 is a critical concentration x„ introduced in Subsection IV-F, below which such an order prevails. For 
< x <C 1 it is a short-range incommensurate-spiral spin order, consistently with the spacing of the square si 
effective lattice reading a s i = \/2/(l — x) a. Such and related types of spin orders have been observed in the cuprate 
superconductor a 26 ' 27 . The schemes introduced in Refs i 28 ' 29 involve modified versions of the square-lattice quantum 
liquid investigated in this paper and contribute to the further understanding of the unusual properties of the hole- 
doped cuprate superconductors^—. 

Our study focus on the Hubbard model on the square lattice. The reason why often our analysis refers to the 
same model on the ID lattice as well is that in contrast to real-space dimensions D > 1 there is an exact solution for 
lD^ir— . In spite of then the model referring to a qualitatively different physics, our quantum-object description also 
applies to ID. For instance, in the U /At — » oo limit the c fermions become both for the model on a ID and square 
lattices non-interacting spinless fermions. In turn, the si fermion occupancy configurations that generate the spin 
degrees of freedom of spin-density m — ground states become in that limit for ID and the square lattice those of 
the spins of the spin-charge factorized wave function introduced both by Woynarovich^l and Ogata and Shiba^ 8 - and 
within a suitable mean-field approximation for a fictitious magnetic field brought about by the electronic correlations 
those of a full lowest Landau level with N sl = N 2 si = N/2 one-sl-fermion degenerate states of the 2D quantum Hall 
effect (QHE), respectively. Here iV 2 is the number of both sites of the square si effective lattice and si band discrete 
momentum values. 

Consistently with the validity of the c fermion, 2i/-?7-spinon r]v fermion, and 2^-spinon sv fermion description of 
Refj^ where ^ = 1,2, ... gives the number of ?7-spinon or spinon pairs, which for the model in the one- and two-electron 
subspace considered in this paper involves only the c fermions and composite two-spinon si fermions, it is shown in 
that reference that for ID and the limit iV a ^> 1, which such a description refers to, the discrete momentum values of 
the c fermions and uv fermions where a = f],s coincide with the quantum numbers of the exact solution. 

In contrast to the original strongly-correlated electron problem, the interactions of the c and si fermions are 
residual owing to their momentum values being good quantum numbers. Consistently, the non-perturbative and 
involved problem concerning the effects of the electronic correlations of the Hubbard model on the square lattice 
simplifies when expressed in terms of the c and si fermion residual interactions. Our results and those of Reb^ 8 . reveal 
that for the model on that lattice the residual interactions recombine the charge and spin degrees of freedom to such 
an extent that one cannot speak of a spin-charge separation as that occurring in ID. For instance, one cannot express 
the electronic spectral functions as a simple convolution of c fermion and si fermion spectral functions. Therefore, 
the concept of spin-charge separation does not apply to the model on the square lattice, at least with the meaning it 
has in ID correlated systems. 

The square-lattice quantum liquid studied in this paper is non-perturbative in terms of electron operators so that, 
in contrast to a 3D isotropic Fermi liqui d 39 i 40 , rewriting its theory in terms of the standard formalism of many- 
electron physics is an extremely complex problem. Fortunately, the description of the physics of such a quantum 
liquid simplifies when it is expressed in terms of the c and si fermion operators. In the x > studies of Refsi 28 ' 29 the 
effects of the residual interactions of the c and si fermions play a major role. At x = the residual interactions of the 
si fermions have no direct effects on the spin spectrum. Within the operator description used in this paper the study 
of the x = spin spectrum refers to an effectively non-interacting problem. In turn, in terms of electrons it is an 
involved many-body problem. The studies of Subsection VI-B confirm the validity of our description for the model on 
the square lattice: At x — our analytical expressions for the spin- wave dispersion, which corresponds to the coherent 
part of the x = spin spectrum, fully agree with the controlled numerical results of Ref. 41 , obtained by summing up 
an infinite set of electronic ladder diagrams. Moreover, an excellent quantitative agreement with the inelastic neutron 
scattering of the La2_ a; Sr a ;Cu04 (LSCO) Mott-Hubbard insulator parent compound La2Cu04 (LCO) 1 ^ is reached. 
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The paper is organized as follows. The model, its c and si fermion description and related rotated-electron 
representation, the limitations and advantages of such a description, and the emergence of the si fermions from the 
hard-core spin-neutral two-spinon si bond particles are the subjects of Section II. In Section III further evidence is 
provided that for the Hubbard model on the square lattice in the one- and two-electron subspace the c and si band 
discrete momentum values are good quantum numbers, consistently with the results of Ref. 9 on that issue. In addition, 
the Fermi line is expressed in terms momenta belonging to the c Fermi line and si boundary line and its anisotropy is 
investigated. The derivation of the c and si fermion energy dispersions of the square-lattice quantum liquid first-order 
energy functional associated with the ground-state normal-ordered Hamiltonian of the Hubbard model in the one- and 
two-electron subspace is the subject of Section IV. The c and si fermion velocities and analysis of their relation to the 
velocities associated with the one- and two-electron excitations is the subject of Section V. Section VI presents a brief 
discussion concerning the combination of the c and si fermion description with the exact Bethe-ansatz solution to 
study the dynamical and correlation functions of the ID Hubbard model. Moreover, in that section we study the spin 
excitations of the half-filling Hubbard model on the square lattice and find excellent agreement between the results 
obtained by use of the square-lattice quantum liquid of c and si fermions and the standard formalism of many-electron 
physics 41 . Finally, the concluding remarks are presented in Section VII. 

II. THE MODEL AND THE c AND si FERMION DESCRIPTION 

Here we address the problem of the description of the Hubbard model on a square lattice in the one- and two-electron 
subspace in terms of charge c fermions and spin-singlet two-spinon si fermions introduced in Ref. 9 . In that subspace 
the model refers to the square-lattice quantum liquid, which is expected to refer to a wider class of many-electron 
problems with short-range interactions on the square lattice belonging to the same universality class. We start by 
introducing the model, its global symmetry, and the rotated-electron representation and discuss the limitations and 
advantages of the corresponding c and si fermion description. 

A. The Hubbard model and its global symmetry, rotated electrons, and the limitations of our description 

The Hubbard model on the two-dimensional (2D) square lattice with torus periodic boundary conditions and the 
same model on the ID lattice with periodic boundary conditions, spacing a, = [N a ] D sites where D — 1 and 
D = 2 for the ID and square lattices, respectively, ~> 1 even and large, and lattice edge length L = N a a for 2D 
and chain length L = N a a for ID is given by, 

H = -t E <V> - + h - C -] + U i N a ~Q}/ 2 - (1) 

Here (rjfji) refers to nearest neighboring sites, the operator ct. creates an electron of spin projection a at the site 
of real-space coordinate fj, and the operator, 

<3 = E E ^^(l-'Vj.-tr), ( 2 ) 

where n?^ a — ct. u Cr j ,a an d — o~ =t (and — a =j.) for a =\ r (and a =y) counts the number of electron singly occupied 
sites. Hence the operator D — [N — Q}/2 counts the number of electron doubly occupied sites where N = N a 

N D 

and N a = Y^j=i n fj,cr- We denote the 77-spin (and spin) value of the energy eigenstates by S v (and S s ) and the 
corresponding projection by — — [N^ —N]/2 (and S* = — [Nf — N±\ /2). We focus our attention onto initial ground 
states with a hole concentration x = [N^ — N]/N® > and spin density m = [JV-j- — N^/N^ = and their excited 
states belonging to the one- and two-electron subspace defined in Ref£. 

The unitary operator V — V(U/4t) associated with the electron - rotated-electron unitary transformation plays 
a key role in the construction of the general description of the Hubbard model on a square lattice introduced in 
Ref£. Out of the manyfold of unitary operators of Refsi 18 ' 21 , the studies of that paper consider a unique choice 
for V — V(U/4:t) such that the states \^u/it) — V^l^oo) are energy and momentum eigenstates for U/At > 0. It 
corresponds to a suitable chosen set {l^oo)} of U/4t —> 00 energy eigenstates. The states \^u/it) — V^|*oo) (one 
for each value oiU/At > 0) that are generated from the same initial state l^oo) belong to the same V tower. The 



■5 



generator S c of the hidden U(l) symmetry of the global 50(3) x 50(3) x U(l) = [50(4) x U(l)]/Z 2 symmetry 
found for the U/At > Hubbard model on a square lattice in RefJ£ reads S c — S c V where S c = Q/2 and the 
operator Q is given in Eq. ([2]). Its eigenvalue S c equals one- half the number of rotated-electron singly occupied 
sites 2S C . The unitary transformation associated with the operator V"f maps the electron operators ct a and c? j>a 

onto rotated-electron creation and annihilation operators ct = ct V" and cV- )CT = V"t Cf?. )CT V', respectively. In 
terms of the latter operators the expression of the generator whose application onto the N = and S c = electron 
vacuum generates the energy eigenstates \^u/u} belonging to the same V tower are the same for the whole range of 
U/At > values: It has the same expression as the generator of the initial state |^oo) hi terms of electron creation 
and annihilation operators. (The N — and S c — electron vacuum corresponds to that of Eq. (|Alj) of Appendix A 
for N° = N® and N° = 2S C = 0.) 

However, for U/At > (and U /At — > oo) the expression of the generators of the energy eigenstates from the electron 
vacuum are very complex in terms of rotated-electron (and electron) operators. Indeed, those are not the ultimate 
objects whose occupancy configurations that generate such states have a simple expression. The studies of Rcf. 9 
considered a complete set {\§u/4t}} of general momentum eigenstates \&u/4t) = ^l^oo) given in Eq. (|A4[) of Appendix 
A. On the right-hand side of that equation the lowest-weight-state (LWS) momentum eigenstates \$LWS-,u/4t) are those 
given in Eq. (|A5|) of such an Appendix. These states can be generated from corresponding U/At — > oo momentum 
eigenstates as \$u/&t) = V^l^oo)- Moreover, \§LWS;U/it) — V^\$lwS;oo) where is the electron - rotated-electron 
unitary operator that also generates the energy and momentum eigenstates \^u/At) = ^l^oo) from the corresponding 
U/At — > oo energy and momentum eigenstates l^oo)- 

The generators onto the vacuum of Eq. (|Al[) of Appendix A of the LWS momentum eigenstates |$ Lws-.u/it) 
are Slatter-determinant products of c and av fermion creation operators. For the ID Hubbard model one has that 
\^u/4t) = V^u/At) so that the states \&u/4t) of Eq. (IA4[) of Appendix A are both momentum and energy eigenstates 
whereas for the model on the square lattice the energy and momentum eigenstates \^u/4t) — V^l^oo) are a superposi- 
tion of a well-defined set of momentum eigenstates \$>u/it) — ^l^oo) with the same momentum eigenvalue, values of 
S v , Sf , S s , Sg, S c , C v = J2 U v N v ,v> C s — J2 U v 'N s ,ui and c fermion momentum distribution function N c {qf . In some 
cases that set of states refers to a single state so that \^u/4t) — \&u/4t) and the state \$u/4t) under consideration is 
both an energy and momentum eigenstate. This is so for the momentum eigenstates \&u/4t) that span the one- and 
two-electron subspace as defined in Rcf^. That \^u/4t) — \^u/u) f° r the ID Hubbard model in the whole Hilbert 
space follows from its integrability being for N a 3> 1 associated with an infinite number of conservation laws^. Ac- 
cording to the results of Ref42 such laws are equivalent to the independent conservation of the set of numbers {N av } 
of av fermions, which for that model are good quantum numbers. 

The av fermion operators /t. av appearing in the Slatter-determinant products of the LWS momentum eigenstates 
\&LWS;U/4t) given in Eq. (|A5I) of Appendix A correspond to an independent problem in each subspace with constant 
values of S c , S s , and number — [N al/ + N^ v ] of both sites of the av effective lattice and discrete momentum 

values of the av band. The operators ft av and f^.au act on subspaces spanned by mutually neutral states, which 
are transformed into each other by av fermion particle-hole processes. As discussed in Ref.- , that assures that for the 
model on the square lattice the components q Xl and q X2 of the microscopic momenta q — [q Xl , qx 2 ] refer to commuting 
si translation generators q s i xt and q s ix 2 - hi turn, for the operators ft c appearing in the Slatter-determinant 
products of the LWS momentum eigenstates there are in ID two types of quantum problems depending on the even 
or odd character of the number [B v + B s ] — N au . Indeed, in spite of the c effective lattice being identical to the 
original lattice, both for the model on the ID and square lattices the c fermions feel the Jordan- Wigner phases of the 
av fermions created or annihilated under subspace transitions that do not conserve the number [B v + B s ] = Yl av Nov- 

The av fermion operators f t av and f^ av act onto subspaces spanned by neutral states. Creation of one av fermion 
is however a well-defined process whose generator is the product of an operator that fulfills small changes in the av 
effective lattice and corresponding av momentum band and the operator ft al , suitable to the excited-state subspace. 
In the Slatter-determinant products of Eq. (|A5p of Appendix A it is implicitly assumed that the av momentum bands 
are those of the state under consideration so that the corresponding generators on the vacua are simple products of 
ft av operators. The c effective lattice and number N^ 1 of c band discrete momentum values remain unaltered for 
the whole Hilbert space. In turn, concerning transitions between subspaces with slightly different discrete values for 
such momenta again creation or annihilation of one c fermion is a well-defined process whose generator is the product 
of an operator that fulfills the corresponding small changes in the c momentum band and the operator ft c or f^ c , 
respectively, appropriate to the excited-state subspace. 

Following the results of Refi^, the bad news is that for the Hubbard model on the square lattice the microscopic 
momenta carried by the av fermions are not in general good quantum numbers, in contrast to the integrable ID 
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Hubbard model. That results from the lack of integrability of the Hubbard model on the square lattice, which is 
behind the set of numbers {N av } of av fcrmions not being in general conserved, yet the numbers — v N^^ 
and C s — ~YLv v N s ,i; an d the c fermion momentum distribution function N c (q) are. It follows that for the model 
on the square lattice the generators onto the vacuum of Eq. (|Alj) of Appendix A of the energy and momentum 
eigenstates are not simple Slatter-determinant products of c and av fermion creation operators as those of Eq. (|A5|) 
of Appendix A and such states are in general different from the momentum eigenstates \$u/4t) — V^l^oo)- The good 
news is that for the model on the square lattice in the one- and two-electron subspace defined in Refi^ one has that 
N a v = except for the av = si and av — s2 fermion branches and owing to symmetries specific to that subspace 
N s i — [S c — S s — 2N S 2] and = 0, 1 are conserved so that the microscopic momenta of the si fermions are good 
quantum numbers and the states \Q>xj/4t) = \^u/4t) are energy eigenstates. 

Addition of chemical-potential and magnetic-field operator terms to the Hamiltonian ([TJ lowers its symmetry. A 
property of key importance follows from such operator terms commuting with that Hamiltonian. Such a property 
reveals that the c fermion and av fermion occupancy configurations and independent 77-spinon and spinon occupancies 
that generate the 4^ momentum eigenstates \&u/4t) = V^l^oo) of general form given in Eq. (|A4[) of Appendix A 
correspond to state representations of the SO(3) x SO(3) x U(l) group for all values of the densities n = (1 — x) and 
m, as confirmed explicitly in Ref.—. 

The unitary operator V commutes with the momentum operator P, three generators of the spin SU(2) symmetry, 
and three generators of the 77-spin SU(2) symmetry 9 -*^. These two SU(2) symmetries are contained in the model 
global £0(3) x 50(3) x U(l) = [SU(2) x ST/ (2) x U(l)]/Z% symmetry. Hence, the momentum operator and such six 
generators have the same expression in terms of electron and rotated-electron creation and annihilation operators. In 
contrast, the generator S c of the U(l) symmetry also contained in that global symmetry and the Hamiltonian ([l} do 
not commute with such a unitary operator. 

The general rotated-electron description introduced in Ref. 9 associates with any operator O a rotated operator 
O = V' O V , which has the same expression in terms of rotated-electron creation and annihilation operators as O in 
terms of electron creation and annihilation operators, respectively. Any operator O expressed in terms of electron 
creation and annihilation operators can then be written in terms of rotated-electron creation and annihilation operators 
as, 

6 = vdv^ = d + [d,s] + h[d,s], §] + ... = 6 + [d,s] + h[d, §],§} + ..., 



T+i - T_ 



+ 0(t 2 /U 2 ); S = -JJ T_,-T_, +Oit-/T-). (:-!) 



The operator S appearing in this equation is such that = e s , V = e~ s , and then [S, V] = and S — S. The 
expression of S involves only the kinetic operators To, and T_i so that that of S involves only the rotated 

kinetic operators Tq, T + i, and T-\. The former kinetic operators are related to the kinetic-energy operator tT = 
~tJ2(rjr ,) Scr=t j,[ c fv, a c »y .er + h.c] of the Hamiltonian ([1]), which can be expressed as f = Tq + T+i + T_i. (The 
expressions of the three kinetic operators are provided in Refi^.) The operator To does not change electron double 
occupancy. In turn, the operators T + \ and T_i do it by +1 and —1, respectively. For U /At > the operator S can be 
expanded in a series oit/U and the corresponding first-order term has the universal form given in Eq. ([3]). To arrive 
to the expression in terms of the operator S also given in that equation, the above property that S = S so that both 
the operators V and S have the same expression in terms of electron and rotated-electron creation and annihilation 
operators is used. This is behind the expansion S = —(t/U) [T+i — T_i] + 0(t 2 /U 2 ) given in that equation for the 
operator S whose higher-order terms involve only products of the rotated kinetic operators Tq, T+i, and T_i. 

Since the Hamiltonian H of Eq. ([1]) does not commute with the unitary operator V — e~ s , when expressed in 
terms of rotated-electron creation and annihilation operators it has an infinite number of terms. According to Eq. 
^ it reads, 

H = VHV^ = H+ [H, S] + i [[H, S], S] + ... . (4) 

The commutator [H, S] does not vanish except for U /At — > 00 so that for finite values of U /At one has that H ^ H . 
For U /At S> 1 the Hamiltonian of Eq. ((4]) expressed in terms of rotated-electron creation and annihilation operators 
corresponds to a simple rotated-electron t — J model. The higher-order t/U terms, which become increasingly 
important upon decreasing U/At, generate effective rotated-electron hopping between second, third, and more distant 
neighboring sites. In spite of the operators To, T + i, and T_i generating only rotated-electron hoping between nearest- 
neighboring sites, their products contained in the higher-order terms of S = —(t/U) [T+i — TLj] + 0(t 2 /U 2 ) generate 
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effective hoping between for instance second and third neighboring sites whose real-space distance in units of the 
lattice spacing a is for the model on the square lattice \/2a and 2 a, usually associated with transfer integrals t' and 
t", respectively^. 

In spite that when expressed in terms of rotated-electron operators the Hamiltonian has an infinite number of terms, 
as given in Eq. ((3]), for intermediate and large U /At values obeying for the model on the square lattice approximately 
the inequality U /At > uo ~ 1.3, besides the original nearest-neighboring hoping processes only those involving second 
and third neighboring sites are relevant for the square-lattice quantum liquid described by the Hamiltonian of Eqs. (JTJ 
and (HJ) in the one- and two-electron subspace^. Hence for such a range of U/At values, only the first few Hamiltonian 
terms on the right-hand side of Eq. Q play an active role in the physics of the Hubbard model on the square 
lattice in that subspace. The results of Ref<2 reveal that for such a model the usual large-?/ physics corresponding to 
energy contributions of the order of t 2 /U is only dominant for U/At 3> 19, so that the higher-order t/U terms become 
important for quite large U/At values. 

It follows from the above analysis that for intermediate and large values of U/At such a square-lattice quantum 
liquid can be mapped onto an effective t — J model on a square lattice with t, t' = t'(U/At) 7 and t" = t"{U/At) 
transfer integrals where the role of the processes associated with t' = t' (U /At) and t" = t"(U /At) becomes increasingly 
important upon decreasing the U /At value. For approximately U /At > uq the latter model is equivalent to the Hubbard 
model on the square lattice (|4]) in the subspace under consideration and expressed in terms of rotated-electron creation 
and annihilation operators. Indeed, the t — J model constraint against double occupancy is in that subspace equivalent 
to expressing the Hubbard model in terms of rotated-electron creation and annihilation operators. 

As discussed in Ref£, the general operator description introduced in that reference for the Hubbard model on the 
square lattice and used in the studies of this paper has two main limitations: 

i) For small and intermediate values of U/At the explicit form of the unitary operator V associated with the rotated- 
electron operators as defined in Ref£ remains an open problem. It is known that such a unitary operator has for 
U/At > the above-mentioned general form V = e~ s where the expression of S involves only the three kinetic 
operators To, T+i, and T_i. The finding of the explicit form of the [//4t-dcpcndcnt functional S in terms of the latter 
three operators, valid for the whole range of finite U/At values, is though a very involved quantum problem beyond 
the reach of our approach current status. Indeed, the operator description used in the studies of this paper has been 
constructed to inherently the solution of that problem being equivalent to the solution of the Hubbard model on the 
square lattice. 

ii) It turns out that the quantum problem under consideration is non-perturbative in terms of electron operators so 
that, in contrast to a 3D isotropic Fermi liqui d 39 ! 40 , rewriting the square-lattice quantum liquid theory emerging from 
the general description used in the studies of this paper for the model in the one- and two-electron subspace in terms 
of the standard formalism of many-electron physics is an extremely complex problem. Fortunately, such a quantum 
liquid simplifies when expressed in terms of the c and si fermion operators. The point is that their momentum values 
are good quantum numbers so that the interactions of these objects are residual. 

A third limitation beyond those discussed in Ref^ is related to the results obtained in this paper: 

iii) The phase factor arising from the extended Jordan- Wigner transformation considered below corresponds to that 
created by a gauge field whose effective vector potential generates long-range interactions between the si fermions 
emerging from the si bond particles. Alike previous theories by many authors involving a gauge theory formulation 
of the Hubbard or t-J model, our method cannot give a fully controllable approximation concerning the effects of the 
interactions brought about by the effective vector potential associated with the gauge field. 

We start by emphasizing that upon expressing the problem in terms of c and si fermions the higher-order con- 
tributions associated with the effective t' = t'(U/At) and t" — t"{U/At) transfer integrals are incorporated in the 
U/At dependence of the c and si fermion parameters studied in this paper and in Ref<2&. Furthermore, concerning 
limitation (ii), the microscopic processes corresponding to the effective t 1 = t'(U /At) and t" = t"{U /At) transfer inte- 
grals are important to characterize the type of order associated with the phases of the square-lattice quantum liquid. 
Fortunately, it is found in Ref£ that for general one- and two-electron operators O other than the Hamiltonian the 
leading operator term O on the right-hand side of Eq. © generates nearly the whole spectral weight. Hence in spite 
of the limitations (i) and (ii), such a description provides useful information about the physics contained in the model 
on the square lattice. Indeed, there are several reasons why, in spite of both the explicit form of the unitary operator 
V being known only for large values of U /At and the difficulties in rewriting the theory emerging from the description 
introduced here in terms of the standard formalism of many-electron physics, that description is rather useful to 
extracting valuable information on the quantum problem for values of U/At approximately in the range U/At > uq. 

Indeed, some of the effects of the unitary operator V are associated with the operator terms of the general expression 
(|3]) containing commutators involving the related operator S = S, which for one- and two-electron operators O are 
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found in Refi^ to generate very little spectral weight. Hence one can reach a quite faithful representation for one- and 
two-electron operators O by expressing the corresponding operator O in terms of the operators of the objects used 
in the studies of this paper whose occupancy configurations generate the energy eigenstates that span the one- and 
two-electron subspace. 

Also concerning the limitation (ii), in Subsection VI-B we consider one of the few physical limits where there are 
reliable and controlled studies of the Hubbard model on the square lattice by means of the standard formalism of 
many-electron physics 4 ^: The spin excitations at half filling. The results of that section reveal an excellent quantitative 
agreement between the results obtained within the square-lattice quantum liquid and those of Ref4^, which profit 
from the usual many-electron machinery and are reached by summing up an infinite set of ladder diagrams. 

Whether our description provides the correct x = and m — spin spectrum is a valuable checking of its validity 
and of the residual character of the c - si fermion interactions. Indeed, within such a description the derivation of 
the spin spectrum is in that physical limit equivalent to a non-interacting problem. That results from both the c and 
si bands being full for the initial x — and m — ground state. As found in Refi^, the spin excitations generate the 
emergence of two holes in the si momentum band but phase-space, exclusion-principle, and energy and momentum 
conservation restrictions prevent inelastic collisions between the two emerging si fermion holes. In contrast, in terms 
of the original electrons it is a complex many-electron problem whose solution requires summing up an infinite number 
of diagrams 4 ^. The excellent agreement reached in Section VI between the spin spectrum as expressed in terms of the 
si fermion dispersions and as obtained by the usual many-electron machinery then confirms the validity of the c and 
si fermion description associated with the square-lattice quantum liquid. 

Finally, concerning the limitation (iii), in our case the momentum band of the si fermions is full for the x > and 
to = ground states and has none, one, or two holes for the excited states that span the one- and two-electron subspace. 
The extended Jordan- Wigner transformation ft gl = e*^ sl gl. gl that maps the si bound-particle operator onto the 

si fermion operator involves the operator phase factor <j)j. s i — fp , s iff /,si ^j'.j^i, which corresponds to that 

created by a gauge field. Fortunately, in our case the operator phase factor 4>j,si can for the study of many properties 
be replaced by the trivial phase factor <^>° sl = Y^j'jkj <A?' ,j,si— • That results from replacing within the present N% — > oo 

limit the operators ft t sl fp.,^i by their average (ft f s x/fy,si) ~ 1- Indeed, for the sl fermion occupancies of the 
energy eigenstates that span the one- and two-electron subspace there are none, one, or two unoccupied sites and 
the total number of sites of the sl effective lattice is iV^ rj N/2 = (1 — x)N%/2. Moreover, the spin degrees of 
freedom of such states are generated by occupancy configurations in the sl momentum band with none, one, or 
two holes whereas the total number of discrete momentum values in that band is also N% rj N/2 = (1 — x)N%/2. 
Since the average value (ft ^ sl /fy, s i) can be expressed as a superposition of expectation values (ft r sl /<f,.,, s i), which 

involves summations over all momenta qf, one then finds (ft r sl /f\,, s i) ~ 1 within the present thermodynamic 

limit. As discussed below in Subsection II-C, replacing ft r sl f? jr , s i by (ft t sl fr jf , s i) ~ 1 renders the form of the 

effective vector potential associated with the extended Jordan- Wigner transformation ft,, sl — e l ^' sl gt. sl that of a 

Chern- Simons vector potential 44 . 

Moreover, phase-space, exclusion-principle, and energy and momentum conservation restrictions dramatically limit 
the effects of the long-range sl - sl fermion interactions generated by the effective vector potential associated with the 
above gauge field. Indeed, these restrictions follow from within the description used here the sl band being full for 
the initial ground state and displaying one and two holes for the one-electron and spin excitations, respectively. This 
prevents such long-range interactions leading to inelastic sl fermion - sl fermion scattering. Hence the direct effects 
of such interactions are somehow frozen. In turn, an indirect side effect of the latter interactions is the occurrence 
of residual interactions between the sl fermions and c fermions. Fortunately, the studies of Refsi 28 ' 29 reveal that the 
description of such interactions is within the range of the theory and that they play an important role in the scattering 
properties of the square-lattice quantum liquid perturbed by small 3D anisotropy effects. In turn, the direct c - c 
fermion interactions either vanish or are very weak. 

For U/At > 0, the investigations of Ref,- accessed the transformation laws and/or invariance of several operators 
and quantum objects under the electron - rotated-electron unitary transformation associated with the operator V, 
what provides valuable information about the physics of the Hubbard model on the square lattice. The suitable use 
of such properties and the combination of the general description introduced in that reference with different methods, 
to extract additional information on the quantum problem, allows the further study in this paper and Ref.~ of the 
square-lattice quantum liquid and its relation to and usefulness for the physics of real materials. 

That quantum liquid refers to the Hubbard model on the square lattice in the one- and two-electron subspace as 
defined in Refi^. For the study of some properties one may consider the smaller subspace spanned by the LWSs of 
both the ?7-spin and spin algebras, whose values of S a and are such that S a — —5* for a — ij, s. (The LWS one- 
and two-electron subspace is a subspace of the LWS-subspace.) In reference^ it is shown that the whole physics of the 
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model (HJ) can be extracted from it in the large LWS-subspace of the full Hilbcrt space. In Appendix A some basic 
information on the one- and two-electron subspace needed for the studies of this paper is provided. The quantum 
liquid of c and si fermions studied in the following is expected to play the same role for many-electron problems with 
short-range interactions on a square lattice as a Fermi liquid for 3D isotropic metal o 39 ' 40 ' 45 . 



B. The algebra of the c fermion operators and hard-core spin- neutral two-spinon si bond-particle operators 

Within the LWS representation, the c fermion creation operator can be expressed in terms of the rotated-electron 
operators as follows 9 , 

4,c = 4, ,t i 1 - + e "' ?J Cfj.t ^4 J "fj.cr = 4>CF JJ(J , (5) 

where here and throughout the remaining of this paper if denotes the momentum of components [tt, 7t] and ir for the 
model on the square and ID lattices, respectively, and e l7r rj is ±1 depending on which sublattice site fj is on. As 
mentioned in Subsection II-A, the rotated-electron operators are related to the original electron operators as, 

ct = yt ct V ; Cft. a = ^ Cft. J, (6) 

where V is the operator associated with the electron - rotated-electron unitary transformation. The unitary character 
of that transformation implies that the operators ct. a and cp j>(7 have the same anticommutation relations as ct. a and 
Oj>.,a- Straightforward manipulations based on Eq. (J5J) then lead to the following algebra for the c fermion operators, 

Let us introduce the corresponding momentum-dependent c fermion operators, 

4,c = -jm E e+i ^' 4- J ^ = 7^E e "^'/^- i = i>->^; * = (8) 

V a j' = i V a 

which refer to the conjugate variable (jj of the c effective lattice real-space coordinate Fj introduced in Refsi&iS. 

The si bond particles are more complex than the c fermions. According to the studies of these references they 
have an internal structure associated with a spin-neutral superposition of two-spinon bond occupancy configurations. 
Therefore, the si bond particle operators gt, sl and gp^si involve a sum of N„\ two-site one-bond operators, including 
N s i/2D of such operators per link family. The concepts of two-site link family and type, which characterize the set of 
two-site one-bond operators whose superposition defines a sl bond-particle operator, are introduced in RefAS.. Alike 
the remaining oa> bond particles and as confirmed in that reference, the sl bond particles have been constructed to 
inherently being hard-core objects so that upon acting onto the sl effective lattice the sl bond-particle operators 
anticommute on the same site of that lattice, 

{9^,81^9^,81} = I! {9^,8^9^,81} ={9f } ,OLV,g^,8l} =0, (9) 

and commute on different sites, 

[9^,81^9^,81] = [<?4, S i'5fy,si] = [9^,81,9^,1] =0; ■ (10) 

The anti-commutation and commutation relations of Eqs. (|9|) and (|10j) . respectively, follow in part from the algebra 
of the spinon operators si and s% 3 introduced in Ref. 10 , which are the building blocks of the two-site one-bond 
operators of that reference. Such spinon operators obey the usual spin-1/2 operator algebra. Indeed and except for 
the U /At ^> 1 limit, for finite values U /At rather than to electronic spins they refer to the spins of rotated electrons 
and according to the studies of Ref. 9 are given by, 

s fj n fj,c Qfj 3 Prj (1 Tlrj.c) Qrj 5 n rj,c 4n c f^j,c , I ^3 ■ (H) 

Here we provide also the expression of the ?7-spinon operator pi associated with the ry-spin SU(2) algebra. Within the 
LWS representation the rotated quasi-spin operators ql can be expressed in terms of the rotated-electron operators 
as follows, 

4=(4,t- e *^t)^ ^ = (4) + , q% = l-n rl ,i, (12) 
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and the c fermion operators are given in Eq. (J5J) . Straightforward manipulations based on Eq. (fT2"j) confirm the 
validity of the following usual algebra for the spinon operators si and s^ 3 , 

H>^.} = 1, {s±,s±}=0, (13) 

[4^j = [44 ]=0; (14) 

and, 




Hence the spinon operators s-. anticommute on the same site and commute on different sites. Consistently with 
the rotated-electron singly-occupied site projector expressed in terms of c fermion operators n? jtC appearing in the 
expression of the spinon operators sf. and s^ 3 provided in Eq. (|11[) . within the — > oo limit their real-space 
coordinates fj can be identified with those of the spin effective lattice. 

The hard-core character of the two-spinon si bond-particle operators associated with the algebra of Eqs. © and 
(fTU| is used in the following in the introduction of the si fermions. 



C. The si fermions emerging from the extended Jordan- Wigner transformation and relation to the 

two-dimensional quantum Hall effect 

1. The extended Jordan- Wigner transformation 

The present study focus mainly on the model on the square lattice, yet the corresponding ID problem is also 
often considered. For instance, the original Jordan- Wigner transformationii transforms ID spin-1/2 spin operators 
into spinless fermion operators. The extension to the square lattice of that transformation has been considered 
previously again for spin-1/2 spin operator o 12 i 13 . In turn, here we apply it to the hard-core spin-neutral two-spinon 
si bond-particle operators of Eqs. (0) and (flU|) . 

We follow the method of Ref»i2 for the spin-1/2 spin operators of both the ID and isotropic Heisenberg model on 
the square-lattice, which was used in Ref.— in studies of the t — J model. We recall that in terms of the creation and 
annihilation rotated-electron operators the spinon operators given in Eq. ([TTT) refer for U/At > to singly occupied 
sites. The corresponding two-spinon si bond-particle operators act onto the si effective lattice defined in Refs^^. 
Each of such objects corresponds to rotated-electron occupancies involving two sites of both the original lattice and 
spin effective lattice defined in RefA 

As a result of the si bond-particle operator algebra given in Eqs. ([§]) and ([TU]) for the si effective lattice of both the 
square-lattice and ID problems, one can perform a Jordan- Wigner transformation that maps the si bond particles 
onto si fermions with creation operators /I. sl and annihilation operators ff^si- Such operators are related to the 
corresponding sl bond-particle operators studied in Ref.— as follows, 

4*i = ^ < h,si = g Pj , sl , (16) 

where 

(j>j, s i = ^ 4>j',j,si ; ^j',j,si = arctan [ J 2 _ 32 J ; < <f)j>,j, s i < 2n . (17) 

jt^j 3 ' \ X j'l X 3\J 

Let Zj — Xj 1 + ixj 2 (with Xj 2 = for ID) be the complex coordinate of the sl bond particle at the sl effective- 
lattice site of real-space coordinate fj and j = 1, ...,N®. For the model on the square lattice the real-space vector 
fj = [xj 1 , Xj 2 ] of Cartesian coordinates Xj 1 and Xj 2 and the complex number Zj — Xj, + ixj 2 are alternative 
representations of the same quantity. The vector difference [fj> — fj] ~ [xj> 1 — Xj 1 , Xj> 2 — Xj 2 ] can be written as, 

[fj, - fj] = \fj, - fj\ e 0j , . isl (18) 
where <^y,j>i is the phase of Eq. (|17l) and here and in the remaining of this paper we denote by, 
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a unit vector whose direction is defined by the angle <j>. 

It follows from the Jordan- Wigner transformation (|16p that the operators ft sl and ffj.si have anticommuting 
relations similar to those given in Eq. (J7J for the c fermion operators, 

{fl^sl ' /Pj'.sl} = fijj' > {fp^sl i /fy.gl} = > /r 3 /,sl} = , (20) 

and the c fermion operators commute with the sl fermion operators. 

The sl fermion operators ft sl and fp j:S i act onto subspaces spanned by mutually neutral states^ corresponding to 
constant values of S c , S s , and number of sites of the sl effective lattice whose expression in terms of the number 
of sites of the spin effective lattice = (1 — x) is for the one- and two-electron subspace given in Eq. (jA17|) 
of Appendix A. Neutral states are defined below in terms of the occupancy configurations of sl fermions carrying 
microscopic momenta. 



2. The extended Jordan- Wigner transformation for the model on the square lattice and relation to the 2D quantum Hall effect 



For the model on the square lattice the problem is much more complex than for the ID model. The phase factor in 
the expressions given in Eq. (|16p corresponds to that created by a gauge field whose effective vector potential A s i(fj) 
and corresponding fictitious magnetic field B s \(fj) read, 



A sl (rj) = $ n V, s i X % —^fl ; B sl (fj) = V r - x A sl (rj) = $ ^ n ? r ,si 5{r r - fj) e Xs , (21) 



3 '3 

where we use units such that the fictitious magnetic flux quantum is given by $o = 1 and 



Pi',},* 



COS (pji jj'jsl 

sin^j/j^i 







are unit vectors perpendicular to and contained in the square-lattice plane, respectively. Moreover, 

- ft 



(22) 



(23) 



is the sl fermion local density operator. The effective vector potential (|21j) can be written in terms of 2D coordinates 
alone as, 



A 



"3 ,3,8 



i + 7r/2 = — arctan 



where the unit vector 



3' ¥=3 



J..1+T/2 



u 3l 



3 1 



b 32 



u 3'2 



has the general form (JTUJ) and is perpendicular to the unit vector e 4 



(24) 



of Eq. 



([18]) . This effective potential generates long-range interactions between the sl fermions. 

The two components of the microscopic momenta of the sl fermions are eigenvalues of the two sl translation gen- 
erators qsix! and q s ix 2 m the presence of the fictitious magnetic field B sl {rj) of Eq. (|2"Tj) given below. Due to the 
non-commutativity of such sl translation generators whose eigenvalues are the components q xx and q X2 , respectively, 
of sl fermion microscopic momenta q = [q xi , q X2 ] of the model on the square lattice one expects that it is impossible to 
classify the states in terms of such momenta. However, in a subspace spanned by mutually neutral states such trans- 
lation generators commute^^. Neutral states are transformed into each other by sl fermion particle-hole processes. 
The contributions of the sl fermion and sl fermion hole to the commutator of the two sl translation generators q s i Xl 
and q s ix 2 cancel in the case of such particle-hole excitations so that the two components q Xl and q X2 of sl fermion 
microscopic momenta q — [q Xl , q X2 ] can be simultaneously specified. The sl fermion operators are defined in subspaces 
spanned by mutually neutral states which correspond to constant S c , S s , and A^ values. Moreover, the values of the 
set of N® x discrete momenta qj — [qj Xl , qjx 2 ] where j — 1, N^ si are well defined for each such a subspace where 
the sl fermion operators of Eq. (|16[) act onto. The corresponding momentum-dependent sl fermion operators are 
given in terms of the operators labelled by real-space coordinates provided in that equation as, 



ft 

J qj,sl 



= E ft, sl ; 4, s i - -jL= e-^'V al ; j = 1, 2V£ ; L = a sl N a=1 



(25) 
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The si band microscopic discrete momentum values tjj where j = 1, ...,N^ play the role of conjugate variables of 
the si effective- lattice real-space coordinates fj which label the corresponding operators of Eq. (fT6)) . Often we omit 
the index j — 1, or j = 1, of the discrete momentum values (jj and denote them by q. According to the 
results of Refi^, the c and si translation generators read q c — ^2^qN c (q) and q sl — <f-/V sl (<7), respectively, where, 

N c (q) = fl c hc; N sl (q) = fl sl fasi, (26) 

are the momentum distribution-function operators. It follows that the si translation generators and q s i X2 

considered above are given by, 

S.i* 4 =X)?*A($ * = 1,2. (27) 

3 

The one- and two-electron subspace contains several subspaces with constant S c , S s , and values. It turns out 
that besides neutral states, which transform into each other by particle-hole processes generated by operators of the 
form f~ sl fq'.si or fq, s i f~, sl , also spin-singlet excited states generated by application onto the m = and x > 
initial ground state of the operator /o s2 fq,si fq'.si where q and q' are the momenta of the two emerging sl fermion 
holes are neutral states, which conserve S c , S s , and iV"^ . For the model on the square lattice the role of the s2 
fermion creation operator /g s2 is exactly canceling the contributions of the annihilation of the two sl fermions of 
momenta q and q' to the commutator of the sl band two sl translation generators g slll and q s \xi of Eq. (|27[) in the 
presence of the fictitious magnetic field B s \ of Eq. (|2"Tj) . so that the overall excitation is neutral. Since the s2 fermion 
has vanishing energy and momentum and the sl momentum band and its number iV^ of discrete momentum values 
remain unaltered^, one can effectively consider that the generator of such an excitation is /<f )S i /<f s i and omit the 
s2 fermion creation, whose only role is assuring that the overall excitation is neutral and the sl fermion microscopic 
momenta can be specified. It follows that for the one- and two-electron subspace the operators /<f, s i fq',si> fqi s \ f^ s ii 

ffai and fq,si f\',si g enerat; e neutral excitations. 

Transitions between states with different values for S c , S s , and/or involve creation or annihilation of single 

sl fermions. The sl fermion operators f~ sl and ,fq. s i refer to subspaces spanned by neutral states. However, as 
mentioned above creation or annihilation of one sl fermion is a well-defined process whose generator is the product 
of an operator that fulfills small changes in the sl effective lattice and corresponding sl momentum band and the 
operator fl sl or fy jS i, respectively, appropriate to the excited-state subspace. 

Upon replacing np^ s \ = ft. sl fp^ s \ by its average ~ 1 the effective vector potential of Eq. (|2ip becomes, 

This is a Chern-Simons like vector potential^^ so that each spin-neutral two-spinon sl fermion has on average 
under the Jordan- Wigner transformation a flux tube of one flux quantum attached to it. That the number of flux 
quanta attached to each sl fermion is odd is consistent with the sl bond-particle operators of the Jordan- Wigner 
transformation /I. sl = e l< ^ sl gt, sl of Eq. (fT6j) being associated with Bose statistics and the corresponding sl fermion 
operators obeying Fermi statistics, respectively. The field strength corresponds on average to one flux quantum per 
elementary plaquette of the square sl effective lattice. Thus for the Hubbard model on a square lattice in the one- 
and two-electron subspace a composite sl fermion consists of a spin-singlet spinon pair plus an infinitely thin flux 
tube attached to it. It follows that within our representation, each sl fermion appears to carry a magnetic solenoid 
with it as it moves around the sl effective lattice. The sl fermions interact with each other via the effective vector 
potential that they create. 

For the model on the square lattice in the one- and two-electron subspace the sl effective lattice spacing magnitude 
a s i = L/N asl « ^/2/(l — x) a is controlled by that of the magnetic length Z s i associated with the fictitious magnetic 
field of Eq. (f2~Tj). Indeed, for such a subspace one has that (tv^si) ~ 1 and the fictitious magnetic field reads 

B s i(fj) w (f>o<5( r j' — fj) £x 3 - It acting on one sl fermion differs from zero only at the positions of other sl fermions. 
In the mean field approximation one replaces it by the average field created by all sl fermions at position fj. This 
gives, 

B s \{rj) ~ $o n s i(fj) e X3 « ®o^jf- e X3 = e X3 ~ $ ° ^ ~ ^ i a li = 27r l li . (29) 
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where l s i is the above fictitious magnetic-field length. Hence the number N% of si band discrete momentum values 
equals [B s i L 2 ]/$ and the magnitude of the si effective lattice spacing a sl is determined by the fictitious magnetic- 
field length l s i as afj = 2%^. This is consistent with for such states each si fermion having a flux tube of one flux 
quantum on average attached to it. 

We find below that the si fermions have a momentum dependent energy dispersion and only in the U/At — > oo 
limit their energy bandwidth vanishes and the N% sl one-sl-fermion states corresponding to the N% = \B S \ L 2 ]/$o 
discrete momentum values are degenerate. Hence, in that limit N% sl = [B s i L 2 ]/$o plays the role of the number of 
degenerate states in each Landau level of the 2D QHE. For the one- and two-electron subspace the si filling factor 
v s \ = N„x/N% reads v s \ = 1 for the m = ground state and v s \ = 1, i/ sl = 1 — l/iV 2 si , or v s \ — 1 — 2/N% 1 for 
the excited states where N% sl ~ [(1 — x)/2] N%- Therefore, within the suitable mean-field approximation (|29]l for the 

fictitious magnetic field B sl of Eq. ([2"Tjl such a ground state corresponds to a full lowest Landau level and in the 
A 2 — > oo limit that our description refers to that applies as well to the excited states. Only for the U/At — > oo limit 
there is fully equivalence between the si fermion occupancy configurations of the one- and two-electron subspace of 
the Hubbard model on a square lattice and the 2D QHE with a full lowest Landau level. In spite of the lack of state 
degeneracy emerging upon decreasing the value of U/At, for finite U /At values there remains though some relation to 
the 2D QHE. The occurrence of QHE-type behavior in the square-lattice quantum liquid shows that a magnetic field 
is not essential to the 2D QHE physics. Indeed, here the fictitious magnetic field arises from expressing the effects of 
the electronic correlations in terms of the si fermion interactions. The si fermion description leads to the intriguing 
situation where the si fermions interact via long-range forces of Eqs. (|2~T1) . (IM|) . and (|2"5)) while all interactions in the 
original Hamiltonian (fl) arc onsite. That combination of our description with the small effects of 3D anisotropy and 
intrinsic disorder is shown in Refs i 28 i 29 to successfully describe the unusual properties of several families of hole-doped 
cuprate superconductors is an indication that QHE-type behavior with or without magnetic field may be ubiquitous 
in nature. 

A site of the si effective lattice occupied by a si fermion corresponds to two-sites of the spin effective lattice. 
Moreover, the = 2S C sites of the spin effective lattice describe the spin degrees of freedom of 2S C rotated electrons 
that singly occupy 2S C sites of the original lattice. In turn, the charge degrees of freedom of such rotated electrons 
are described by the 2S C sites of the c effective lattice that are occupied by c fermions^. 

Since the c fermions and si fermions describe the charge and spin degrees of freedom, respectively, of the same sites 
singly occupied by rotated electrons, the long-range interactions between the si fermions generated by the effective 
vector potential of Eqs. (|2"Tj) . and (|2"8)l are felt by the c fermions and lead to residual interactions between those 
and the si fermions, whose effective interaction energy is derived in Ref* 2 ^. For instance, according to Table HVl of 
Appendix A upon creation and annihilation of one electron, one c fermion is created and annihilated, respectively, and 
one si fermion hole is created. In contrast to the ID case, the residual interactions brought about by the extended 
Jordan- Wigner transformation are for x > behind inelastic collisions between for instance a si fermion going over 
to the momentum value of the si band hole created upon creation or annihilation of one electron and the c fermions 
with momenta near the c Fermi line introduced in this paper. 

3. The effects of the si fermion Jordan- Wigner transformation in ID 

The above discussion is different from ID where the quantum problem is intcgrable and due to the occurrence for 
N a 3> 1 of an infinite number of conservation laws, the c and si fermions have zero-momentum forward-scattering 
only. That the physical consequences of the extended Jordan- Wigner transformation through which the si fermions 
emerge from the si bond particles are different for ID and 2D is consistent with for ID the number Zj = xj 1 + ixj 2 
involving the coordinates Xj 1 and Xj 2 appearing in Eq. (|17l) being such that Xj 2 — and hence then reducing to 
the real-space coordinate Xj l of the si bond particle in its effective lattice. Therefore, for ID the phase 4>j',j, s i = 
arctan([xj' 2 — Xj 2 ]/[ x j'i ~ x j\\) m that equation can have the values <pj',j )S l = and fij'j.si — n only. Indeed, the 
relative angle between two sites of the ID si effective lattice can only be one of the two values. It follows that in ID 

e *(^+i,»i— = e "^j. sl f*j- sl _ s i fermion discrete momentum values Qj are eigenvalues of the si translation 

generator q sl = J2 q lN sl (q) in the presence of the fictitious magnetic field B sl of Eq. (|2T1) associated with the 
Jordan- Wigner transformation. Their values are given by the exact solution and can be expressed as a sum of a bare 
momentum of the usual form [2w/L] Afj 1 where Afj 1 — 0,±1,±2, ... are integers and a small deviation q®i/N s \ that 
for a subspace spanned by mutually neutral states has a constant value. This gives, 



qj = ^M 3 sl +q° sl /N sl - A^ sl =j-^ = 0,±1,±2,...; i = l,..,iV Qsl 



(30) 
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Concerning the value of q^ it turns out that since for ID the above phase 0j'j, s i can have the values 4>j',j.si — and 
4>j',j,s\ — 7T only, there are as well two types of subspaces only where it is given either by q° sl = or q° sl — n[N s i/L], 
respectively, for all j = 1, ...,N aBl discrete momentum values. The si effective lattice length L — N asl a s % equals that 
of the original lattice. Here a s i = L/N asl = [N a /N aBl ] a is the si effective lattice constant. It follows that in ID the 
transitions between subspaces with different si fermion discrete momentum values qj are associated with deviations 
given by 5q®i = ±tt[N s i/L]: Under such subspace transitions all the discrete momentum values qj of Eq. (|30[) are 
shifted by the same overall momentum Sq^/Ngi = ±tt/L. 

Each site of the si (and s2) effective lattice occupied by a si (and s2) fermion refers to the spin degrees of freedom 
of two (and four) sites of the original lattice singly occupied by rotated electrons whose charge degrees of freedom 
are described by two (and four) c fermions. In ID the sharing of the same rotated-electron sites by the si (and s2) 
fermions and c fermions, respectively, leads to residual interactions between such objects whose main effect is on the 
boundary conditions of the momentum values of the charge c fermions: Due to such a residual interaction the latter 
objects feel the Jordan- Wigner phase of the si and s2 fermions created or annihilated under a subspace transition so 
that their discrete momentum values can be written in the form, 

q 3 = ^M J c + q° c /N c ; AT? = j - ^ = 0, ±1, ±2, ... ; j = l,...,N a . (31) 

Again, there are only two types of subspaces where gj? is given either by qj? = or gj? = n[N c /L], respectively, for all 
j = 1, N a discrete momentum values of the c band. For the subspaces of the one- and two-electron subspace there 
are finite occupancies of c and si fermions plus one or none s2 fermion of vanishing energy and momentum so that the 
general equation for the phase factor e lqj L = e l7T l B v+ B ^} = e ™H av N av involving the qj c band momenta considered 

in Ref<^ simplifies to e lqi L = e m ^ N " 1+N " 2 \ It follows that the exact-solution quantum numbers 7? = [Af? + f^r-^H 
considered in that reference are integers and half-odd integers for subspaces spanned by states with [N„i + N^] even 
and odd, respectively, where consistently with the values of the number [2S S +2N S 2] given for the one- and two-electron 
subspace in Eq. (|A16|) of Appendix A one has iV s2 = for subspaces of spin S s = 1/2 or S s — 1 and N s2 — 0, 1 for 
subspaces of spin S s — 0. In Ref£ it is confirmed that such an effect results indeed from the phase <j>j jS i of Eq. (|17[) . 
which controls the si fermion Jordan- Wigner transformation. In turn, the energy is not affected, in contrast to the 
model on the square lattice. 

Analysis of the quantum numbers associated with the exact solution, which are related to the c and av description 
in Ref3, reveals that in ID the good quantum numbers are not the integer numbers Af? 1 — j — N asl /2 = 0, ±1, ±2, ... 
and Nj = j — N a /2 = 0, ±1, ±2, ... of Eqs. (|30|) and (|3"T1) . respectively, but rather the corresponding shifted numbers 

if = [A/ - / 1 + ^j^} and J| = [Af? + f^] where j = l,...,N asl and j = l,...,A a , respectively. For a given 
subspace the numbers Z| (and if ) are either consecutive integers Z| = J\f$ (and if = A^j 51 ) or half-odd integers 
Ij = [Afj — 1/2] (and if — [N"f — 1/2])- Such an effect is associated with boundary conditions controlled by the 
Jordan- Wigner phases. 

Hence in ID the eigenvalues qj — ([2tt/L] Afj 1 + q^/Ngi) of the si translation generator in the presence of the 
fictitious magnetic field are the good quantum numbers. Such conserving numbers are the si band discrete momentum 
values of Eq. ([3H]) . However, that as a side effect the si fermion Jordan- Wigner phase of a si (and s2) fermion created 
or annihilated (and created) under subspace transitions is felt by the c fermions so that all c band discrete momentum 
values are shifted by 5q®/N c = ±ir/L is not a trivial effect. (Note that creation or annihilation of two si fermions or 
creation or annihilation of one si fermion plus creation of one s2 fermion does not lead to changes in the c band discrete 
momentum values.) This confirms that for ID the Jordan- Wigner transformation that generates the si fermions from 
the si bond particles gives rise to zero-momentum forward-scattering interactions between the emerging si fermions 
and the pre-existing c fermions. 

That the Jordan- Wigner phases of the si fermions are felt by the c fermions under si fermion creation or annihilation 
through the residual interactions associated with their sharing of the same sites of the original lattice occurs for the 
model on the square lattice as well. In ID such a residual interactions involve only zero-momentum forward-scattering. 
In contrast, for the square lattice they lead to inelastic collisions involving exchange of energy and momentum, as 
confirmed in Refsj 28 i 29 . 



III. QUANTUM NUMBERS OF THE HUBBARD MODEL IN THE ONE- AND TWO-ELECTRON 
SUBSPACE: THE c AND si BAND DISCRETE MOMENTUM VALUES 

In this section we address the problem of the c and sl band discrete momentum values, including that of the relation 
of the Fermi line to the c Fermi line and sl boundary line introduced in the following. 
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The si bond particles studied in Refsi&i - correspond to well-defined occupancy configurations in a spin effective 
lattice that for the model on the square lattice is a square lattice. The si bond-particle description of RefJ^ involves 
a change of gauge structur e - 10 ' 24 so that the real-space coordinates of the sites of the si effective lattice correspond 
to one of the two sub-lattices of the square spin effective lattice. In Appendix B it is confirmed that for the ^> 1 
limit that our study refers to and for hole concentrations x such that the density n = (1 — x) is finite the two choices 
of si effective lattice lead to the same description. 

A. The states generated by c and si fermion occupancy configurations 

For the model in the one- and two-electron subspace the form of the general states \$u/4t) — V^l&oo) of Eq. (|A4[) of 
Appendix A whose LWSs \&LWS;U/4t) also appearing in that equation are those given in Eq. (|A5[) of that Appendix 
simplifies so that such a subspace is spanned by a x > and m = ground state and its one- and two-electron excited 
states of the form, 

(Qt)L s _i/2 Nal 2S <= 

\*U/4t) = ^ S V \*LW Si U/*t); \*LWS;U/4t) = K, K^U 4 ,1 ^ ' K)] iTI 4 J GS ^ ^ ' ( 32 ) 

Here the coefficient C s is provided in Eq. (|A4[) of Appendix A, according to Table HVl of that Appendix the number 
of —1/2 independent spinons is given by L s ,_i/2 = for the initial ground state and some of its excited states, 
^s, -1/2 — 1 f° r one-electron excited states, and L s _i/ 2 = 2 both for spin-triplet excited states and excited states 
involving addition or removal of two electrons with parallel spin projections, and J} 2 ^ /I. C \GS C ; 0) is a generalization 
of the state given in Eq. (|A2[) of Appendix A. Indeed, only for the spin density m = (1 — x) of that equation are the 
c fcrmions invariant under the electron - rotated-electron unitary transformation for all values of U /At. In turn, for 
m = such an invariance occurs only for U /At — » oo. 

As shown in Refi^, for the Hubbard model on the square lattice in the one- and two-electron subspace as defined in 
that reference the numbers N s \ of si fermions, N^ of sites of the si effective lattice and thus of discrete momentum 
values of the si fermion momentum band, and N S 2 of zero-momentum and vanishing energy spin-neutral four-spinon 
s2 fermions are good quantum numbers. They read N s i = [S c — S s — 2N S 2], N^ sl — [N®/2 + S s ] = [S c + S s ], and 
the s2 fermion number is given by N S 2 = for all excited states of spin S s = 0, S s = 1/2, and S s = 1 except for the 
S s — spin states for which N S 2 — 1. Hence the number of si fermion holes — [N^ — N s i] — [2S S + 2N S 2] is 
also conserved for the square-lattice model in such a subspace. For the ID model all such numbers are good quantum 
numbers in the whole Hilbert space. 

Moreover, according to the results of Ref3 the general states \§u/it) an d \^LWS;U/At) of Eqs. (|A4|) and (|A5[) . 
respectively, of Appendix A are for the model on the square lattice momentum eigenstates yet in contrast to ID they 
are not in general energy eigenstates. Indeed, for the model on the square lattice neither the set of av fermion numbers 
{N av \ nor the corresponding numbers {JV^ } of av band discrete momentum values are in general conserved. Also 
in contrast to ID, in general the corresponding set of av translation generators q av in the presence of the fictitious 
magnetic fields B av considered in Refi^ do not commute with the Hamiltonian, yet commute with the momentum 
operator. Since the components qj xl and <7 Jx2 of the av band discrete momentum values qj are eigenvalues of the 
corresponding av translation generators q a vx 1 and q a vx 2 i such momentum values are not good quantum numbers and 
thus are not in general conserved, again in contrast to ID. Consistently, the states of Eqs. (|A4I) and (|A5[) of Appendix 
A generated by momentum occupancy configurations of such av fermions are not in general energy eigenstates. 

In turn, since for the model on the square lattice in the one- and two-electron subspace the set of av fermion 
numbers {N au } are conserved and read N s \ = [S c — S s — 2^2], N s2 = 0, 1, N sv = for v > 3, and N VL , = for 
v > 1 and the two si translation generators q s i Xl and q s ix 2 °f Eq. (|27p in the presence of the fictitious magnetic 
field B s i of Eq. (|2"TT) commute with the Hamiltonian, the states (|3"2"j) are both momentum and energy eigenstates 
Consistently, it follows from the algebra of the c and si fermion operators of Eqs. (JTJ) and (|20l) and the relations 
of Eqs. ([5]) and (|25p that for subspaces with constant values of S c , S s , and A^ S 2 the set of states of form (|3"2"j) are 
orthogonal and normalized. In addition, both for the model on the square and ID lattices they span the subspaces of 
the one- and two-electron subspace for which the numbers (S c , S s ) are given by (N/2, 0) and N s2 — or A^ S 2 — 1, as 
well as its subspaces with N s2 = and either (N/2 ± 1/2, 1/2), (N/2, 1), (N/2 ± 1, 0), or (N/2 ± 1, 1) where N is the 
electron number of the initial m — ground state. In such subspaces the states of form (|3"2")) have been constructed 
to inherently referring to a complete basis. 

Creation or annihilation of one si fermion is a well-defined process whose generator is the product of an operator 
which fulfills small changes in the si effective lattice and corresponding si momentum band and the operator /t sl 
or fq t si, respectively, suitable to the excited-state subspace. In the Slatter-determinant products of Eq. (|3"2"j) it is 
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implicitly assumed that the si momentum band is that of the state under consideration so that the corresponding 
generators on the vacua are simple products of /t sl operators. The c effective lattice and number of c band 
discrete momentum values remain unaltered for the whole Hilbert space. Concerning transitions between subspaces 
with different values for such momenta, again creation or annihilation of one c fermion is a well-defined process 
whose generator is the product of an operator that fulfills small changes in the c band momentum values and the 
operator ft c or f^ c , respectively, suitable to the excited-state subspace. The operators of the vanishing- energy and 
zero-momentum +1/2 independent 77-spinons, +1/2 independent spinons, and s2 fermion do not appear explicitly in 
the expression of the generators onto the vacua of Eq. (|32| because the effects of such operators are take into account 
by the above-mentioned two operators that fulfill small changes both in the si effective lattice and corresponding 
si momentum band and in the c momentum band, respectively. Since in the case of the general states (|32p the 
operators are created onto the c and si momentum bands of the state under consideration such effects are accounted 
for implicitly. 

In the particular case of the one- and two-electron subspace as defined in Ref.— , the general momentum expression 
introduced in that reference simplifies to, 



Here N c (q) and N s ±(q) are the expectation values of the momentum distribution- function operators of Eq. (1261) . That 
for the Hubbard model on the square lattice in the one- and two-electron subspace the discrete momentum values qj 
of the c and si fermions are for U/4t > good quantum numbers implies that the interactions of such objects have 
a residual character. That allows one to express the excited-state quantities in terms of the corresponding deviations 
SN c (qj) and SN s i(qj) of the c and si fermion momentum distribution functions, respectively, relative to the initial 
ground-state values. Alike in a Fermi liqui d 39 ' 40 , one can then construct an energy functional whose first-order terms 
in such deviations refer to the dominant processes and whose second-order terms in the same deviations are associated 
with the c and si fermion residual interactions. Such a functional refers to the square-lattice quantum liquid whose 
c and si fermion energy dispersions are studied in this paper. 

The momentum area (2ir/L) 2 iV„ and discrete momentum values number N% 3l = [S c + S s ] of the si band are 
known. In turn, its shape remains an open issue. The shape of the c band is the same as for the first Brillouin zone, 
yet the discrete momentum values may have small overall shifts under variations of [N s \ + N S 2[ generated by subspace 
transitions. Indeed, alike in ID each si (and s2) effective lattice site occupied by one si (and s2) fermion corresponds 
to the spin degrees of freedom of the rotated-electron occupancy of two (and four) sites of the original lattice whose 
charge degrees of freedom are described by two (and four) c fermions. Such a sharing of the original-lattice sites 
by c fermions and si or s2 fermions, respectively, is behind the c fermions feeling the Jordan- Wigner phases (|17|) of 
the si or s2 fermions created or annihilated under subspace transitions. In the studies of this paper we use several 
approximations to obtain information on the c and si bands momentum values and the corresponding c and si fermion 
energy dispersions of the model on the square lattice in the one- and two-electron subspace, for which the c and si 
fermion discrete momentum values are good quantum numbers. In turn, the s2 band does not exist for N S 2 = 
subspaces such that . = and for the N S 2 — 1 subspace the s2 momentum band reduces to a single vanishing 
momentum value. The corresponding s2 fermion has both vanishing momentum and energy. 

The main difference relative to the usual Landau's Fermi liquid theor y 39 ' 40 is that upon adiabatically turning off the 
interaction U the c and si fermions do not become electrons. Indeed, the non-interacting limit of the present objects 
corresponds instead to t 2 /U — > 0. For the c and si fermions the second-order terms of the above energy functional are 
associated with residual interactions. However, the effects of the original electronic correlations are accounted for both 
by its first-order and second-order terms. In Section IV the ground-state normal-ordered Hamiltonian is described 
in terms of the square- lattice quantum liquid of c and si fermions with residual interactions and the corresponding 
energy functional is studied up to first order in the c and si fermion hole momentum distribution-function deviations. 
The second-order energy terms of that functional are discussed in Ref.-* for an extended square-lattice quantum liquid 
perturbed by small 3D anisotropy effects. The studies of this paper do not include the effects of the phase fluctuations 
associated with the superconducting order of that extended quantum problem, which are found in that reference to 
lead to a new first-order contribution. 



B. Subspace-dependent deviations of the c and si band discrete momentum values in ID and 2D 

Since according to the analysis of Ref. 9 the c and si band discrete momentum values are for the Hubbard model on 
a square lattice in the one- and two-electron subspace good quantum numbers they deserve further investigations. In 
contrast to the momentum values of the quasiparticles of a Fermi liquid, the number = [S c + S s ] of si fermion 




(33) 
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discrete momentum values may change upon state transitions involving variations of S c and/or S s . In addition to 
changes in the number of discrete momentum values, subspaces spanned by states with different AT? = [S c + S s ] 
values may have slightly different sets of discrete momentum values. The number of discrete momentum values 
of the c band remains unaltered for the whole Hilbert space. However, alike in ID such momentum values may also 
slightly change upon variations of the number [B s + B v ] = Y2 av N a v, which for subspaces contained in the one- and 
two-electron subspace simplifies to [B s + B v \ — [N s \ + N s2 ] — [S c — S s — N s2 ]. Alike in ID, for the Hubbard model 
on the square lattice the c fermions feel the Jordan- Wigner phases of the si (and s2) fermions created or annihilated 
(and created) under subspace transitions that do not conserve the number [B v + B s ] = ^2 av N av . The c and si 
band momentum shifts originated by transitions between such subspaces considered in the following refer both to 
the model on the ID and square lattices. In turn, the si band boundary line deformations, which in the present 
N% — > oo limit can be described by si band rotations, as well as the concept of si fermion doublicity studied in the 
ensuing subsection are specific to the one-electron excited states of the model on the square lattice and thus to the 
corresponding square-lattice quantum liquid. 

The c fermion operators ft c and f$. c refer to a c momentum band whose momenta are related to the real-space 
coordinates of the c effective lattice through Eq. © . The latter lattice is identical to the original lattice and in the local 
c fermion occupancies whose superposition gives c fermion momentum occupancy configurations that generate the 
states l|32p the c fermions occupy the sites singly occupied by the rotated electrons in the latter lattice. Consistently, 
for all energy eigenstates the momentum area (2ir/L) 2 N 2 and shape (and momentum width (2ir/L) N a ) of the c band 
are (and is) for the model on the square (and ID) lattice the same as for the first Brillouin zone. Out of the 
discrete momentum values qj where j = 1, 2V„ , N c = 2S C are occupied and Nf? — [N^ — 2S C ] are unoccupied and 
correspond to c fermions and c fermion holes, respectively. 

In turn, the si fermion operators /i sl and f^ s i are associated with a si momentum band related to the si effective 
lattice through Eq. As confirmed below, such a si momentum band is exotic. It is full for the x > and m = 

ground states considered here so that there is no si Fermi line and as given in Table ITVl of Appendix A for one- and 
two-electron excited states it has one hole and none or two holes, respectively. Moreover, for the model on the square 
lattice the shape of its boundary line is subspace dependent. That line is strongly anisotropic for that model, its 
energy dependence being d-wave like: As confirmed in Section IV, it vanishes for excitation momenta pointing in the 
nodal directions and reaches its maximum magnitude for those pointing in the anti-nodal directions. 

Let <fj sl (and ±?Bsi) be momentum values belonging to the boundary line (and being the two boundary points). 
For the model on the square (and ID) lattice such a line encloses (and such points limit) the si band Nj? discrete 
momentum values. The quantity d = ±1 in the boundary- line momentum <fj sl of the model on the square lattice is 
the doublicity. It refers to two alternative subspace-dependent shapes of that line, which for the subspace spanned by 
one-electron excited states result from well-defined si boundary line deformations relative to the initial ground-state 
si band, as discussed below in Subsection III-C. The si fermions whose occupancy configurations generate a given 
one-electron excited state have the same doublicity d. In spite of the c band not being deformed upon one-electron 
excitations, for x > the doublicity of c fermions at or near the c Fermi line is considered to be that of the si fermion 
holes created along with the creation or annihilation of the c fermion under consideration. Hence, the doublicity 
d is the same for the si fermion holes and c fermions created or annihilated upon an one-electron excitation. The 
doublicity d = ±1 and the corresponding two types of sl-band boundary line shapes occur only for one-electron 
excited states whereas both for x > and m = ground states and their two-electron excited states one finds in that 
subsection that q£ sl = q^si = QBsi- 

For the model on the square lattice the si boundary-line momentum values <fj sl are well defined for each subspace 
spanned by states with constant U /At, S c , S s , and thus N£ = [S c + S s ] values. The precise form of such a boundary 
line is in general an open problem, yet we know that it encloses a momentum area given by (2ir/ L) 2 N 2 1 . Fortunately, 
the problem can be solved for the x = and m = absolute ground state, as confirmed in Section IV-E. Below we also 
address it for some of the subspaces spanned by excited states generated from application of one- and two-electron 
operators onto x > and m = ground states. 

Alike in ID, one can in general decouple the components qj x and qj X2 of the discrete momentum values qj of the si 
fermions into a bare part of the form [2n/L] Mi where i = 1, 2 and A/i = 0, ±1, ±2, ... and a small momentum deviation 
whose value is the same for all momentum values. As in ID, the separation of qj in a bare component [2n/L] Mi 
and a small deviation (qj . — [2ir/L]Afi) is just a matter of convenience. The same separation can be fulfilled for the 
components of the c band microscopic momenta so that for the model on the square lattice the c and si band discrete 
momentum values may be written in following general form, 

qj = kj+tf/N^ 

<Usi= 3%*/^! A/7 = 0,±l,±2,...,; 7 = C ,sl; j = 1, ; i = 1,2, (34) 
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where N% c = N%. 

Since for all energy and momentum eigenstates (|32[) the momentum area and shape of the c momentum band are 
the same as for the original Brillouin zone, the integer numbers Aff associated with the Cartesian components of the 
bare momenta kj appearing in this equation read, 

M?=ji-^-=0,±l,±2,...,; j = {j 2 -l)N a +j 1 = l,..,,Nl; j t = 1, ...,N a ; i = 1, 2 , (35) 

and are such that \J\ff\ < N a /2. In turn, the shape of the si band is subspace dependent so that the limiting values 
of the integer numbers A^ sl associated with the Cartesian components of the bare momenta kj of Eq. ([34]) depend 
on the shape of its boundary line for each subspace and value of U /At. In turn, the number of discrete momentum 
values qj of that equation is always given by N% = [S c + S s ] so that j = 1, ...,N% for all subspaces with constant 
[S c + S s ] value. 

Subspaces with different values of S c , S s , and N% — [S c + S s ] have in general different si band discrete momentum 
values qj as well so that in addition to a change in the number N% — [S c + S s ] of such values the corresponding 
subspace transitions may involve effects such as si band overall shifts by the same small momentum Sq^/N^ and 
si band boundary line deformations, which in the present N 2 — > oo limit may be described by si band rotations. 
For momenta near the si boundary line such rotations are defined by the F angle introduced in Subsection III-C. 
Subspace transitions involving changes in the [S c — S s — N S 2] value may involve c band overall shifts by the same 
small momentum 5q® /N 2 . Both for the model on the ID and square lattices symmetry implies that the momentum 
shift Sq^i/Ngi originated by subspace transitions is such that, 

E£ = 0, (36) 

3=1 

both for the initial and final subspaces. Consistently, for the model on the square (and ID) lattice the sl band refers 
to a compact momentum area {2tt / L) 2 [N asl ] 2 (and momentum width (2tt/L) [N asl ]) centered at zero momentum. 
The physical meaning of Eq. ([56]) is that when a subspace transition leads to a deviation SN^ = S[S C + S s ] = ±1 
associated with addition (+1) or removal (—1) of a discrete momentum value qj to or from the sl momentum 
band, respectively, the corresponding momentum shift 5q~±/N s i is such that Sq^ — T'Tjo- That exact cancelation of 
the momentum ±Qj by the momentum shift generated by the subspace transition is a symmetry equivalent to the 
sum-rule (f3"6"|) always being fulfilled. 

The momentum <f c ° of Eq. (|3^|) is well defined for the subspaces of the model on the square lattice spanned by 
states with the same [S c — S s — N S 2] value. However, one cannot in general access it by the use of symmetry alone, 
except in some cases. For instance, for hole concentrations x > its shift <5<f c ° can be derived for transitions from 
x > and m = ground states to excited states generated by one-electron processes. For an excited state generated 
by addition of one electron onto such a ground state one has that 5q® = 0. In turn, for an excited state generated by 
removal of one electron from the same initial ground state it is found in the ensuing subsection that the value of the 
corresponding shift <5<f c ° is such that the Fermi line has the same form as for one-electron addition. 

C. Fermi line of the square- lattice quantum liquid 

1. Three hole concentration ranges according to the Fermi-line anisotropy 

According to Luttinger's theorem, the momentum area (and width) enclosed (and limited) by the Fermi line (and 
points), when centered at zero momentum, is for the Hubbard model on the square (and ID) lattice (1 — x)2ir 2 (and 
(1 — x)n). Given the non-perturbative character in terms of electrons of the Hubbard model on both the square and ID 
lattices, such Fermi line and points, respectively, could correspond to the generalized Luttinger line and points defined 
in Ref. 46 . That line (and points) is (and are) determined not only by poles in the one-electron Green's functions, but 
also by their zeros. It is argued below that for x > x c \ where x c \ ~ 0.13 for intermediate U/At values the Fermi line 
of the model on the square lattice obeys Luttinger's theorem. 

The non-perturbative effects are strongest in the case of the model on the ID lattice. Fortunately, for the ID 
Hubbard model one can profit from its exact solution to study the problem. Then we define the Fermi points as the 
momenta where the real part of the one-electron Green's function has strongest singularities at very small excitation 
energy ui. For < x < 1, < m < (1 — x), and U/At > such singularities occur at the [7 = Fermi momenta ±fci? CT , 
as given in Eqs. (25) and (26) of Ref4£. For the ranges < x < 1, m = 0, and U/At > considered here the two 
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exponents provided in Eq. (26) of that reference are identical and given in Table II of the same paper for magnetic 
field H — > 0. Then the singularities occur for U /At > at the U = Fermi momenta ±kp — ±(1 — x)ir/2 so that the 
Fermi points limit a momentum width (1 — x)tt for all values of U/At, as required by the usual Luttinger's theorem. 

For the ID Hubbard model the one-electron renormalization factor vanishes at the Fermi points as defined above. 
For general correlated 2D problems, if such a factor is finite for the whole Fermi line one defines that line in the usual 
way. In turn, if it is finite for some Fermi-line parts and vanishes for other parts, it could be defined as for ID in the 
latter parts. Finally, if the one-electron renormalization factor vanishes for the whole Fermi line it could be defined as 
the set of momenta belonging to a line in momentum space where the real part of the one-electron Green's function 
has strongest singularities at very small excitation energy Li, alike in ID. 

Let us consider a first Brillouin zone of the model on the square lattice centered at the momentum —7? = [— tt, — tt] 
and define the Fermi line in terms of a hole Fermi momentum kp related to the Fermi momentum kp as, 



kp = k F + tt = kp((f>) ; tt = 

4> G [4>AN, ? - &ANJ 5 (ij- + <t>AN,n ~ 4>AN^j 



K Fx 2 



= arctan [ ; h 2 J ; <f> 6 (0, 2ir) , < x < Xh 



K F Xl 



37T , \ / 3n 



'JAN, — <PAN 



— + <I>an,^-4>an), x h <x<x*, (37) 



2 r / V 2 

where <j> is the Fermi-line angle that defines the direction olkp and we assume that the hole concentration Xh above 
which the Fermi line is particle like belongs to the range x C 2 < Xh < x*. Here the hole concentration x C 2 obeys the 
inequality 2x^/3 < x C 2 < x* and is defined in Rcf. 28 upon application of a uniform magnetic field aligned perpendicular 
to the square-lattice plane yet it has physical meaning in the absence of such a field as well, the critical concentration 
x* in Subsection IV-F, and the angle 4>an = 4>an{x) > vanishes for x < Xh and is small (or vanishes if xt > x*) for 

x e (xh, x*). 

Consistently with the numbers and number deviations given in Table IIVI of Appendix A for one-electron excited 
states and results of Refi^, creation of one electron at the Fermi line leads to creation of one c fermion at the c Fermi 
line and one hole at the si boundary line so that the hole Fermi momentum kp of Eq. (|37p can be written as, 

^ = [fto-?li]=4Wef (38) 

Here the c fermion Fermi momentum qp^ and si boundary-line momentum <fj sl are given below and as confirmed in 
Section V the doublicity d — ±1 refers to two types of one-electron excited states with the same hole Fermi momentum 
kp and Fermi energy Ep but different Fermi velocity Vp. 

For < x < x* and intermediate U/At values obeying approximately the inequality uq < U/At < u\, where 
i*o ~ 1.302 is the U/At value at which the spinon pairing energy 2Ao introduced in Ref^ corresponding to < x <C 1 
has its absolute maximum magnitude and u\ ~ 1.6, the Fermi line level of anisotropy is controlled by the interplay 
of the si boundary line anisotropy and c Fermi line isotropy, as discussed in the following. The degree of Fermi-line 
anisotropy has impact on the form of the spectrum of the energy eigenstates. For instance, for the reduced one-electron 
subspace defined below the doublicity d = ±1 labels one-electron excited states with the same energy and momentum 
but different d-dependent velocities. However, we find in Subsection V-B that for hole concentrations where the 
Fermi-line anisotropy is small these two velocities are d independent and equal. Since the doublicity effects on the 
spectrum are for the electronic-operator description non-pcrturbative, that non-perturbative character is strongest 
for hole concentrations where the Fermi-line anisotropy is largest. We find below that the Fermi-velocity anisotropy 
is small for x > x c i where x c \ ~ 0.13 for the intermediate U/At values U /At € (uq, Ui). According to our scheme, that 
the Fermi-velocity anisotropy is small implies that for x > x c \ the Fermi line encloses the momentum area (1 — x)27r 2 
predicted by Luttinger's theorem. This is consistent with the experimental studies of Refi^, which for some hole-doped 
cuprates find that the Fermi line obeys Luttinger's theorem for approximately x > 0.10. In turn, we find below in 
Subsection IV-E that for the x = and m — Mott-Hubbard insulating phase the Fermi line has a shape independent 
of U/At and encloses the expected momentum area 2ir 2 , in spite of the Mott-Hubbard gap. However, our results are 
inconclusive concerning the validity of Luttinger's theorem for the range x G (0,x c i) where the Fermi- line anisotropy 
is largest. Indeed, our investigations about that Fermi line do not focus on that range. 

The anisotropy of the Fermi line results from that of the si boundary line. For si band momenta belonging to 
that line the absolute value of the si fermion velocity Vg sl = V^(qg sl ) derived in Section V is for < x < x* 
dependent on ip. The corresponding Fermi energy has the form Ep — fi + SEp where [i is the zero-temperature 
chemical potential and SEp = |A||cos2(/>| the anisotropic part arising from the si boundary-line energy SEp. It 
reaches its maximum magnitude max-j^E^} = |A| for hole Fermi momenta pointing in the anti-nodal directions and 
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vanishes when such momenta point in the nodal directions. Here |A| is the energy scale introduced in Ref^ whose 
< x <C 1 expression |A| « Ao(l — x/x*) given in that reference is argued in Subsection IV-F to be valid for the hole 
concentrations < x < x* in the approximate range uq < U/At < u n where x* increases smoothly from x» « 0.23 
at U/At = uq to x* w 0.32 at U/At = > u\. Whether for such a U/At range the hole concentration Xh at which 
that line becomes particle like is given by Xh ~ x C 2, is such that x C 2 < Xh < x», is given by Xh ~ x*, or is larger 
than x* remains an open question. In the studies of this paper and Refs . 28 ' 29 we assume that x C 2 < x/j < x*. For 
x > Xh the Fermi line is particle like. Hence one has that max{<5i?_F} = |A|| 008 20^^1 where 4>an = 4>an{x) is the 
angle of Eq. Q37p. which vanishes for x < Xh and is small for x € (x^, a;*). As discussed in Section V, we consider 
only corrections of first order in 4>an so that maxj^E'i?} = |A| for < x < x* and uq < U/At < u n , consistently 
with max {SE F } « |A|[1 - 2(j) 2 AN ] « |A| for x h < x < x*. In turn, v£ sl = V£(q£ sl ) reaches its maximum magnitude 
IAI/v/2 for momenta pointing in the nodal directions and for < x < x* vanishes for momenta pointing in the 
anti-nodal directions. If Xh < x* one finds instead that for Xh < x < x* the minimum value of V^ sl is small but finite 
and given by |A|| sm.2(j) an\ / ^/2 ss [4>an V2]\A\. 

The anisotropic Fermi energy part 8Ef is found in Section V to equal the </>-dependent energy of the si fermion 
dispersion at the si boundary line. The absolute value Vp 1 of the Fermi velocity is found to be a function of both 
the si fermion velocity Vg sl = VJi(q^ s -,) and c fermion velocity Vp c = V c (qp^) such that at the hole concentration 
range for which maxVg sl /VFc <C 1 one has that Vp ~ Vf = V c . Such an equality is not met the Fermi velocity Vp 
is dependent on both (f> and d = ±1 and the Fermi line is that of an anisotropic system. Such a dependence on cf> 
follows from that of the si boundary-line velocity Vg sl on that angle. In turn, for approximately x > x c \ one finds 
that ma,xVg sl /Vf c <C 1 and Vp ~ Vp — Vf c , where the c fermion Fermi velocity absolute value is <p independent 
and given by Vf c ~ [\/x7r 2/m*] so that the corresponding Fermi line is that of an isotropic system. Here m* is 
the c fermion energy-dispersion mass introduced below in Subsection IV-D. The Fermi- velocity anisotropy coefficient 
r/A = max V^^/Vfc provides a measure of such a Fermi- line anisotropy, that line fully becoming that of an isotropic 
system as tja — >• for x — > x*. 

A complementary coefficient, which also provides a measure of the Fermi- line anisotropy is the Fermi-energy 
anisotropy coefficient r) = \A\/W^. The zero-temperature chemical potential /i = [Ep — 8Ep] is for < x < x* and 
uo < U/At < u\ approximately given by /i m /i° + Wj? . This expression is exact for small hole concentrations x <C 1 
and a good approximation for the range < x < x* provided that U/At e (u ,«i). In addition to /i° = lim :E _ ! .o /i, 
it involves the energy scale = [8t - W%] where W/? is the energy bandwidth of the c fermion filled sea of Ref^. 
The energy bandwidth vanishes in the limit x — > and corresponds to a contribution to the zero-temperature 
chemical potential of the rotated-hole kinetic energy. Indeed, we recall that the c fermion holes describe the charge 
degrees of freedom of the rotated-electron holes whose concentration x equals that of the electronic holes. Hence Wj? 
is the isotropic Fermi energy measured from the upper level of the unfilled band. Indeed, the chemical potential range 
/i € (— /io, A*o) refers to half filling so that for x > one has that /i > /io- Therefore, for x > rather than | A|//j, with 
+ Wj?, the coefficient of interest to provide a measure the Fermi-energy anisotropy is r/o = \A\/W^. 

It is found in Section IV that in the units used here the mass m* appearing in the c Fermion energy dispersion 
reads m* — \/2r c t where the ratio r c — /m* involves the U/At oo bare c fermion mass = l/2t. For 
approximately the range uq < U /At < u± the magnitude of such a ratio is given by r c « 2r s where the spin ratio 
r s = Ao/AW®! is introduced in Ref£ and studied below in Subsection IV-F. For such a range of U/At values they 
read r c 2r s 7rx* where the critical hole concentration x* is studied in that subsection. The expressions of the 
anisotropy coefficients given in the following can be extended to the range U/At £ (uq, u^) where > u\ is the U/At 
value considered below in Subsection IV-F at which x* = 1/tt sa 0.32. However, here we provide such expressions for 
the approximate range U/At € (uo,ui) where their form simplifies. For < x < x* and the latter U/At range the 
coefficients are given by, 

fxT V£ sl \A\ l2x~A~(l 1\ 1 /A \ 2 

X ° = 4r c 2 f+ A X * ; Xcl = X * [Fc ~ ~ 1] 5 Fc = 1+ 7r 2(i + A C /4r c 2 2 ; 2X * ,Z ~ ~ ** ' ^ 

The hole concentrations xa and xq also given here are those at which ?7a = 1 and t]q = 1, respectively, and the hole 
concentrations x c \ that at which 7/a = 2xo- For hole concentrations above x c i the Fermi- velocity coefficient ?7a is 
considered small. In the present context the hole concentration x c2 defined in Ref. 28 is that above which both the 
anisotropy coefficients t]a and 770 are considered yo be small. The expressions Vf c ~ Wxir 2/m*] and ~ [x7t2/to*] 
used to derive the expression of such coefficients given in Eq. (|39[) are a good approximation for the range < x < x„ 
provided that Uq < U/At < u\. Indeed, for small U /At the c fermion mass m* introduced in Subsection IV-D becomes 
large and behaves as 1/m* — > for U/4t — > so that for < U/At -C 1 the expressions Vp c ~ [\fx~n2/m* c ] and 
W% « [x7r2/m*] are valid for x <§C 1 only whereas for U/At > u\ the hole concentration x* is larger than 0.28. 
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U/At 


XA 


Xo 


X C 1 






1.9 x nr 3 


3.0 x 10~' J 


0.12 


0.23 




1.1 x 10~ 3 


2.4 x 


0.13 


0.27 



TABLE I: Magnitudes of the hole concentrations xa, xo, and x c i given in Eq. f|39|) and critical hole concentration x, introduced 
in Subsection IV-F for U/At — Uo ~ 1.302 and U/At = u* = 1.525. The interplay of the Fermi-line anisotropy and the physics 
of the different hole-concentration ranges discussed in the text refers to intermediate U/At values uo < U/At < ui. 

The magnitudes of xa, xq, and x c \ and that of X* rs 2r s /ir rj r c /Tt are provided in Table I for U/At = uq 
and U/At — u*. The value U/At — u* — 1.525 is found in Refs i 28 ' 29 to be suitable for the description of several 
families of hole-doped cuprates with superconducting zero-temperature critical hole concentrations x = x c ~ 0.05 and 
x = x» « 0.27. The magnitudes of the coefficients 7/a and 770 are for U/At ~ u* provided in Table II for several values 
of the hole concentration. To reach those and the xa and xq values given above we used the magnitude Ao ~ 0.285 1 
found in Subsection VI-B for U/At s» w*. 

We emphasize that the hole concentrations xa, xq, and x c i do not mark sharp phase transitions. They refer instead 
to crossovers between different levels of Fermi-line anisotropy occurring for approximately the range Uo < u < u\. 
Whether x C 2 marks a crossover or refers to a critical point remains an open question 2 ^. The Fermi line anisotropy 
has important effects on the physics of the square-lattice quantum liquid and, therefore, these hole concentrations 
mark crossovers between regimes where some of the physical properties are different. Often we focus mainly our 
analysis on the value U/At ~ w* = 1.525 € (u 0l u{) of interest for the physics of several families of hole-doped cuprate 
superconductors 2 ^ 2 ^. For uo < U /At < u\ it is considered that the Fermi- velocity anisotropy is small when rjA < 2xo 
and thus x > x c \ where rjA ~ 0.048 and x c \ &s 0.13 for U/At ~ and that above the hole concentration x C 2 both 
the Fermi-velocity and Fermi-energy anisotropies are very small. The hole concentration x C 2 marks the beginning of 
the crossover to a particle-like Fermi line and the accurate value Xh > x C 2 at which the Fermi line becomes actually 
particle like remains an unsolved problem. If as assumed above Xh < x* then for < x < Xh the Fermi line angle 4> 
belongs to the range given in Eq. (|37p. where 4>an = <Pan{x) is a small angle that vanishes for x < Xh- 

The above analysis leads to three main ranges of hole concentrations where the system has different levels of Fermi- 
line anisotropy. The first range refers to hole concentrations smaller than x c i ~ 0.13. Then according to Table II 
both coefficients are larger than 2xq ~ 0.05 and the Fermi-line anisotropy is not small and plays an important role 
in the physics. This is the anisotropic Fermi-line hole-concentration range. It may be further divided into three 
subranges corresponding to < x < xa, xa < x < xq, and xq < x < x c \, respectively. For < x < xa both 
maxV^ sl > Vfc and SEp > W]} so that the Fermi-line anisotropy is huge. In the range xa < x < xq it remains 
large because SEp > Wj? yet one has that max V^ sl < Vf c . Finally, the range xq < x < x c i marks the transition of 
a system with large Fermi- velocity anisotropy at x « xo to one with small Fermi- velocity anisotropy at x x c \. The 
strong Fermi-line anisotropy of the range < x < x c \ renders the study of that line a complex problem, which we do 
not address here. 

Fortunately, the study of the Fermi line simplifies for the isotropic Fermi-velocity hole-concentration range x c % < 
x < x C 2, where the Fermi- velocity anisotropy is small and that line hole like. For that hole-concentartion range the 
coefficient is small whereas t]q is not as small. This range can be further divided into the subranges corresponding to 
x c i < x < x op and x op < x < x C 2 where the Fermi-energy anisotropy is not small and starts being small, respectively. 
Here x op — 0.16 for U/At « 1.525. The hole concentration x op that separates such subranges is found in Reb 2 ^ to 
play the role of optimal hole concentration, where the critical temperature of the superconducting phase of the related 
square-lattice quantum liquid perturbed by small 3D anisotropic effects is largest. (Such effects refer to weak plane 
coupling whereas the anisotropy considered here is that of the 2D Fermi line.) The Fermi velocity is for the isotropic 
Fermi-velocity hole-concentration range nearly independent of the Fermi-line angle <j> so that concerning quantities 
involving mainly that velocity the system behaves as isotropic. There remains though some level of anisotropy, which 
has effects on some properties. The Fermi-energy coefficient rjo provides a measure of such a remaining isotropy. 
Indeed, that coefficient decreases slower upon increasing x than r]A , yet both coefficients vanish in the limit x — > x* , 
where the Fermi line fully becomes that of an isotropic system. That nearly occurs in the third range considered here, 
x c i < x < x* , where both coefficients are small and the Fermi line is approximately that of an isotropic system. This 
is the isotropic Fermi-line hole-concentration range. (For U/At w 1.525 one has that x» w 0.27.) 

For intermediate U/At values uq < U/At < u\ the expressions introduced below for the Fermi line are valid for 
the isotropic Fermi-velocity hole-concentration range x c \ < x < x C 2 where the angles that define the directions of 
the momenta qp£ and <?J sl in expression Q38[) have a simple relation to the Fermi-line angle </>. As mentioned above, 
the experimental studies of Ref.— find that for some hole-doped cuprates the Fermi line does not obey Luttinger's 
theorem for approximately x < 0.10. For the hole-doped cuprate value U/At « 1.525 our study of the Fermi line is 
restricted to the range x c i < x < x C 2 where x c i w 0.13 so that we cannot judge on that interesting issue. In turn, 
our results agree with those of RefjS in that for approximately x > 10 _1 the Fermi line obeys Luttinger's theorem. 
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X 


X 


0.11 


0.13 


0.15 


0.16 


0.21 


0.23 


0.24 


0.27 




0.20 


0.06 


0.05 


0.04 


0.03 


0.02 


0.01 


0.01 





no 


1.00 


0.14 


0.11 


0.08 


0.07 


0.03 


0.02 


0.01 






TABLE II: Magnitudes of the Fermi-velocity anisotropy coefficient 77a and Fermi-energy anisotropy coefficient 770 of Eq. (|39p 
for different hole concentrations x and U/4t = u* = 1.525, where xo ~ 0.024. It is considered that the Fermi- velocity anisotropy 
is small when t]a < 2xo w 0.048 and thus x > x c i w 0.13. 

Anyway, for smaller hole concentrations the above issue concerns mostly how one defines the Fermi line. The zeros 
of the one-electron Green's function on the real axis signify poles in the self-mass operator and the derivation by 
Luttinger involves integration of it over frequencies. The trouble is that the result depends on how these singular 
integrals are regularized. The Fermi line is mostly a zero-temperature concept. Nevertheless the regularization used 
in ReL~ starts, before turning the temperature to zero, by considering finite temperatures and performs the Luttinger 
calculations in the Euclidean sector of the energy-momentum space where there are no poles in self mass. However, 
there are at zero temperature other possible regularizations as for instance that considered in Ref.— . 

Within our approach, for x > x c \ the Fermi line encloses the momentum area (1 — x)2t{ 1 predicted by Luttinger's 
theorem. In turn, the simple relation of the angles that define the directions of q/pf and <fj sl in expression (|38|) to the 
Fermi- line angle <f> is valid only for the isotropic Fermi- velocity range x c \ < x < x C 2- For such a x range the Fermi line 
is hole like and each si band momentum value q at or near the si boundary line referring to states belonging to the 
reduced one-electron subspace defined below has a simple relation to exactly one hole momentum value k = k + tt 
at or near the Fermi line. Indeed, q is the momentum of the hole emerging in the si band upon removal from a 
772 = initial ground state of one electron of hole momentum k h at or near the Fermi line. As discussed below, for 
that hole-concentration range each momentum of the Fermi line is generated from a si band momentum q at the si 
boundary line and a c band hole momentum q h at the c Fermi line. 

2. The reduced one-electron subspace and the doublicity 

An important exact property is that one-electron (and two-electron) excitations do not couple to excited states 
with an even (and odd) number of si fermion holes^. The transformation laws of the eigenvalue spectrum of the si 
translation generators q sl in the presence of the fictitious magnetic field B s i of Eq. (j21"]) under subspace transitions 
from a to = ground state to excited states involving creation of an odd number of si fermion holes are such that the 
si boundary line is deformed. We emphasize that the momentum shifts Sq^/N s i of the momenta q2 /N 1 of Eq. 
for the 7 — si band momenta under subspace transitions do not deform the si boundary line in the present N% — > oo 
limit. Indeed, the shift of each such momenta is of the order of l/iVf and vanishes, only the resulting macroscopic 
momentum q^ due to the shift of all N s ± occupied si band discrete momentum values being of physical importance, 
assuring that the sum-rule (|36|) is fulfilled. 

In turn, the changes in the eigenvalue spectrum of the si translation generators q sl considered here are of another 
type, leading to deformations of the si boundary line. Subspace transitions from a generating m — ground state lead 
to such si boundary line deformations only when the final subspace is spanned by excited states with an odd number 
of si fermion holes. Indeed, the transformation laws of the eigenvalue spectrum of the si translation generators q sl 
are such that the contributions from each of the emerging si band holes to the si boundary line deformation cancel 
when the number of such holes is even. In turn, when it is odd one of such contributions remains uncompensated 
so that the si boundary line is deformed. The latter subspace transitions involve creation of an odd number of c 
fermions or c fermion holes as well and when only one si fermion hole is created, the resulting excited states couple 
(and do not couple) to one-electron (and two-electron) excitations. 

Here we address the occurrence of such deformations of the si boundary line for the particular case of subspace 
transitions involving creation or annihilation of one electron. In the present N% — > oo limit, Tio < U/4t < u%, and 
for momenta near the si boundary line such deformations can be described by si band rotations, which for hole 
concentrations in the range x c \ < x < x C 2 preserve the absolute values of such momenta but change their angles. 
Similar deformations occur for x < x c \ but owing to the larger Fermi-line anisotropy their study is a more involved 
problem, which we do not address here. Similar mechanisms occur for x > x C 2 but again are nor studied here due 
to complications associated with the transition from a hole like to a particle-like Fermi line. For the present hole- 
concentration range x € (x c i, x C 2) such rotations occur in the square-lattice quantum liquid under subspace transitions 
from the initial m = and (Nf, N±) — (N/2, A^/2)-electron ground state to the subspace spanned by a (N/2 + 1, N/2)- 
electron or (N/2,N/2 — l)-electron ground state and its excited states of small momentum and low energy. All such 
states have a single hole in the si momentum band, as confirmed by the numbers provided in Table HVl of Appendix 
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A. They are one-electron excited states of the m = ground state, being generated from it by removal (or addition) 
of a spin-down (or spin-up) electron. Relative to the latter ground state such states may have finite momentum and 
finite energy. In turn, for two-electron excited sates of the initial m — ground state such rotation-like si band 
deformations do not occur. 

Most of the results of this subsection refer to the square-lattice quantum liquid in a well-defined subspace of its one- 
and two-electron subspace: The reduced one-electron subspace. Some of our results are valid only for the isotropic 
Fermi- velocity range x € (x c \,x C 2)- For each ground state with (N/2,N/2 — 1) electrons there is a corresponding 
initial m — and (N/2, iV/2)-electron ground state with one more spin-down electron. We say that such a ground 
state is the m — generating ground state of both the (N/2, N/2 — l)-electron ground state and its excited states of 
small momentum and low energy. For x > the latter states are generated from the (N/2, N/2 — l)-electron ground 
state by particle-hole processes in the c band. The reduced one-electron subspace is spanned by such a (N/2, N/2 — 1)- 
electron ground state plus its excited states of small momentum and low energy. (Similar results are obtained for a 
(N/2 + 1, iV/2)-electron ground state.) The reduced one-electron subspace is a subspace of the one-electron subspace 
where the square-lattice quantum liquid is defined. 

All sites of the si effective lattice of a m = generating ground state are occupied so that the si momentum band 
is full. In turn, the si effective lattice of both a (N/2 + 1, iV/2)-electron or (N/2,N/2 — l)-electron ground state 
and its excited states of small momentum and low energy have a single unoccupied site and thus the si momentum 
band has a single hole. In the present N 2 — > oo limit, the deformation of the si boundary line generated under the 
transition from the m = generating ground state to the reduced one-electron subspace can be described by a well- 
defined rotation of the si band about an axis perpendicular to the plane, crossing the zero momentum. Symmetries 
imply that the corresponding rotation angle is proportional to the hole concentration x so that the si band remains 
unchanged for the one-electron excited states generated from the x = and m — ground state. 

For finite hole concentrations one finds within a continuous representation of the si fermion discrete momentum 
values corresponding to (2tt/L) 2 — > within the N 2 — > oo limit that the si band momenta q of the reduced one-electron 
subspace are related to the momenta <fo of the corresponding m — generating ground state as follows, 

q = A d sl q + [2/(1- x)][5q^/N 2 a \; lim ^ = I ; d = ±1 . (40) 

x-s-0 

Here is a 2 x 2 matrix, I the unit 2x2 matrix, d = ±1 the doublicity, and [2/(1 — a;)] [Sq^/N 2 ] the small momentum 
shift Sq^/Nsi of Eq. ([34]) for 7 = si. The transformation (|40|) describes a deformation of the si boundary line under 
a subspace transition from the m = generating ground state to the reduced one-electron subspace. It refers to 
momentum values q at or near the si boundary line and for hole concentrations in the isotropic Fermi- velocity range 
x c i < x < x C 2 and uq < U/At < u\ it is described by a si band rotation by a well-defined angle <fip, as discussed 
below. 

Here we consider the particular case of states with none, one, or two unoccupied sites in the si effective lattice. 
The 2x2 matrix is well defined though for states with an arbitrary finite number of unoccupied sites. Due 
to symmetry it reduces to the unit matrix, A^ = I, for states with an even number of unoccupied sites in the si 
effective lattice. Therefore, such a matrix is the unit matrix for subspace transitions from a m = ground state to 
subspaces spanned by energy eigenstates with a finite even number of holes in the si band. In turn, the derivation 
of its form for states with an odd number of holes in the si band is in general a complex problem. For the one- and 
two-electron subspace of the square-lattice quantum liquid the states with an odd number of holes in such a band 
are those having a single hole and therefore here we limit our analysis to the reduced one-electron subspace, which is 
spanned by one-sl-band-hole states. 

The momentum shift [2/(1 — x)][Sq^i/N^] vanishes in the thermodynamic limit. The small c band momentum shift 
Sq® /N c — [1/(1 — x)][Sq^ /N 2 ] of the momenta of Eq. (f34|) has a similar limiting behavior. When one considers a 
single si band discrete momentum value q or a single c band hole momentum value q h one may ignore the shifts 
[2/(1 — xf^Sq^/N 2 ] or [1/(1 — x)][6q® /N%], respectively, since they refer to contributions of order 1/N 2 , which 
vanish in the thermodynamic limit. However, the corresponding excitation momenta Sq~i and <5<f c ° resulting from 
the contributions of all N s i si fermions and N c c fermions, respectively, are macroscopic momenta, which must be 
taken into account. For instance, for the si band the main role of the momentum Sq^ is to assure that the sl-band 
momentum sum rule (|36[) is fulfilled for the excited state under consideration. Therefore, the discrete momentum 
values are shifted by [2/ (1 — x^Sq^/N 2 ] after the si band transformation generated by the matrix A d sl is performed. 

The matrix A~ t of Eq. (|40f is for x — > given by lim^^o^fi = I- For hole concentrations x £ (x c \,x C 2) it can 
be evaluated for momentum values at and near the si boundary line. The generation of the excited states of small 
momentum and low energy from the above (N/2 + l,N/2) or (N/2,N/2 — 1) ground state involves only changes in 
the occupancy configurations of the si fermions at and near the si boundary line. 

For a m = generating ground state the si boundary-line momenta have the following general form, 



qBsi = 9s s i(^) e^ sl ; </> s i = + 7r. 



(41) 
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This result is valid for the whole range of hole concentrations < x < considered in the investigations of this 
paper. The term tt of <f> s i = <p + tt follows from the si band momenta belonging to the quadrant such that k Xl < 0, 
k X2 < 0, and < <f> < tt/2 of the Brillouin zone centered at —if = [— 7r, — tt] pointing in a direction defined by an angle 
(p s i in the range tt < <f>„\ < 3tt/2. 

Since the si boundary line encloses a momentum area (1 — x) 2tt 2 , the momentum absolute value qssi = qBsi((p) 
of Eq. (|4ip obeys the following normalization condition and has the following limiting magnitudes, 

2tt ji „ 



(I 



-^-7r[ tos i(</))] 2 = (1-x)2tt 2 ; q Bs i{<p) = , — . , x <C 1 ; ffu.ifa) = y/(l-x)2n, (1 - z) < 1 

Z7r cos © + sin © 



(42) 

It is confirmed below that at x — the si band coincides with an antiferromagnetic reduced Brillouin zone such that 
I foil + 1 9x 3 1 < 7r, consistently with the limiting behavior given in Eq. (|42|) for the corresponding boundary line when 
x < 1. 

That upon increasing the value of a; the si band remains full for x > and to = ground states yet encloses a 
smaller momentum area (1 — x)2ir 2 is associated with the short-range incommensurate-spiral spin order addressed 
in Ref.— . As discussed in Section I, this is consistent with the spacing a s i = y/2/(l — x) a found for the square si 
effective lattice in that reference such that the finite-energy real-space excitations involving si fermions with momenta 
pointing near the anti-nodal directions break the translational symmetry. 

As mentioned above, the momentum area enclosed by the si boundary line when centered at zero momentum is 
(1 — x)2ir 2 and by c Fermi line when centered at zero momentum is (1 — x)4tt 2 (and on the momentum —7? = — [tt, tt] 
is xAtt 2 ). Such properties are valid for the hole-concentration range < x < x* of interest for the square-lattice 
quantum liquid. However, that the momentum area enclosed by the Fermi line when centered at zero momentum is 
(1 — x)2ir 2 (and on the momentum —tt = —[tt, tt] is (1 + x)2tt 2 ) is a property that holds at least for x > x c \. These 
properties of such three lines together with (i) the Fermi- velocity anisotropy being small and (ii) the Fermi line hole 
like imply that for si band momenta q at or near the si boundary line the matrix A d x of Eq. (|40l) is for x c \ < x < x c i 
orthogonal and generates a rotation by a well-defined x dependent angle (p F . For each electronic hole momentum 
k h = k + tt at or near the Fermi line of Eqs. ([37| and (|38l) there are two values of q at or near the si boundary line 
corresponding to two different one-electron states whose si fermions have different doublicity d = ±1, as discussed 
below. 

For a si band momentum <fo = q(<p) e^ sl of the initial m — ground-state subspace at or near the si boundary line 
the corresponding momentum q of Eq. (I40|) of the reduced one-electron subspace si band is at or near that line as 
well. The derivation of the F angle relating these two momenta takes into account the emergence of a single hole in 
the si band of the reduced one-electron subspace. From the interplay of the hole momentum k h at or near the Fermi 
line of the removed or added electron, momentum of the c fermion removed or added, respectively, and momentum 
of the hole emerging in the si band one finds by imposing that the Fermi line, c Fermi line, and si boundary line 
enclose the momentum areas given above that for hole concentrations in approximately the isotropic Fermi-velocity 
range x c \ < x < x C 2 the si band momenta q at or near the si boundary line and corresponding momenta tfo of the 
initial generating ground state si band have the general form, 



q = q(4>) e^ si+K + [2/(1 -x)] ; <f = q(cb) ; 4(0) = darctan f ^^j gW -j ; d = ±1 , (43) 
where <pp((p) is the F angle. 

Alike the general expression given in Eq. (|40|. those provided in Eq. (|43|) are valid for the reduced one-electron 
subspace. For q at or near the si boundary line the expression provided in the latter equation reveals that the matrix 
of Eq. (|40|) is orthogonal for the range x c \ < x < x C 2 where the Fermi- velocity anisotropy is small and the Fermi 
line hole like. It reads, 



Ad 
Ap 



cos (p F — sin (p F 
sin (p d p cos <p F 



(44) 



Here the F angle given in Eq. (|43p can have two values, (pp 1 and <pp^. The two corresponding rotations refer to the 
doublicity d = — 1 and d = +1, respectively. The expression provided in Eq. (|43|) refers to si band momentum values 
at or near the si boundary line. The doublicity also labels the si fermions for all remaining momenta of the si band, 
which are also of the form (l4Tfl) yet the precise form of the matrix A d x is in general unknown. Therefore, all N sl si 
fermions whose occupancy configurations generate the same state have the same doublicity. 
The si boundary line momenta have the following general form, 

QbsI = lB S i{<P)e^ ■ (45) 
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Such a boundary line is rotated relative to that of the corresponding m — generating ground state with one more 
electron by the angle 4> F = 4> F ((j)) so that the si angle <j> s i of Eq. (|4H is for x c \ < x < x C 2 replaced by the angle 
on the right-hand side of Eq. (j45|) given by, 

& = [</>.! + <&] = [<!> + * + 4]- (46) 

In turn, for x > and m = generating ground states and approximately uq < C//4i < u w the c Fermi line is hole 
like for hole concentrations below i « r s > x, where x* = 2r s /n, as found in Section IV. Since the square-lattice 
quantum liquid scheme refers to the range x £ (0,a;*), for it the c Fermi line is hole like. Such a c Fermi line exists 
for x > and we define it in terms of hole c Fermi momenta q Fc centered at —7? = [— 7r, — ir] and given by, 

Qfc = QFc + 7? = q F c(<P) ; 4>c = (t>- (47) 

As mentioned above, when centered at —7? = [— 7r, — 7r] the c Fermi line encloses a momentum area x 4ir 2 so that for 
< x < x* the momentum absolute value q Fc — q Fc ((f>) obeys the sum rule, 

f =x4tt 2 ; x<x*. (48) 

As also mentioned above, for the hole-concentration range x c \ < x < x C 2 and uo < U/At < u\ one can label the 
c fermions created or annihilated at or near the c Fermi line upon one-electron excitations by the doublicity value 
d = — 1 or c? = +1 of the si fermion hole created upon the same excitation. According to the numbers and number 
deviations given in Table HVl of Appendix A, removal of one electron from a m = generating ground state leads to 
creation of a c fermion hole and a si fermion hole. In turn, addition of one electron to such a ground state leads to 
creation of a c fermion and a si fermion hole. The latter process does not change the c band discrete momentum 
values. Moreover, one can suitably choose the S c = N/2 value of the initial in = ground-state subspace so that 
q° = in Eq. ([34]). 

The reduced one-electron subspace si band momenta q of Eqs. (|4T)|) and (|43p are in these equations related to the 
corresponding momentum values <fo of the m = generating ground state. The changes in the si band originated 
by subspace transitions involving removal of one electron are simpler to describe than those involving one electron 
addition. Indeed, as confirmed by analysis of the numbers and number deviations given in Table HVl of Appendix A, 
for one electron removal the si band of the m — and (A/2, A/2) generating ground state has the same number of 
discrete momentum values as that of the (A/2, A/2 — 1) ground-state subspace, so that Eq. (|4*U| refers exactly to 
an one-to-one relation. In turn, the si band of the generating (A/2, A/2) ground state has one discrete momentum 
value less than that of the (A/2 + 1, A/2) ground state. Nevertheless, in the present thermodynamic limit the number 
of sites of the si effective lattice and thus that of discrete momentum values of the si band is N% w A/2 so that 
for x < x* such a difference is a effect and the momentum values q of Eqs. (|40| and (|43|) could as well be 

expressed in terms of those of the m = generating ground state with less one electron than the corresponding 
reduced one-electron subspace. 

We confirm below that within the c and si fermion description of the one-electron problem an identical hole Fermi 
line is generated by removal and addition of one electron from and to a m = generating ground state, respectively. 
Since the Fermi line generated by either process is identical, we start by considering addition of one electron. We 
emphasize that independently of the actual exact form of the momentum absolute value functions q Fc (<f) and (7Bsi(<^), 
for uo < U/4t < u\ and the hole-concentration range x € (x c i,x C 2) where the Fermi velocity anisotropy is small, 
Luttinger's Theorem holds, and the Fermi line is hole like, the hole Fermi line generated by creation of one electron 
and thus creation of one c fermion at the c Fermi line and one si fermion hole at the si boundary line encloses the 
correct area provided that the momenta of the two created objects are perpendicular to each other. It follows that 
for the model on the square lattice the Fermi line results in either case from si fermion and c fermion excitations 
that involve a well-defined relation between the angles defining the directions of the momenta belonging to the si 
boundary line and the c Fermi-line hole momenta, respectively. 

That for instance for electron addition the c fermion created under the subspace transition has doublicity follows 
from for the hole-concentration range x € (x c i,x C 2) the hole momentum of such a c fermion and the momentum of 
the si fermion hole created under the one-electron excitation at the c Fermi line and si boundary line, respectively, 
being perpendicular to each other. Indeed, depending on the doublicity of the si fermion hole, the c fermion hole 
momentum points in different directions. Those correspond to two alternative c Fermi momenta q F ^ where d = ±1 
has the same value as the momentum <fj sl of the si fermion hole under consideration given in Eq. (|45[) . Hence in spite 
of the c Fermi line remaining unaltered under the subspace transition associated with the one-electron excitation, for 
hole concentrations x c \ < x < x C 2 there are for d = ±1 two alternative c Fermi momenta given by, 



qk d = q h Fc (.<P d c)e^, 



(49) 
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where the momentum absolute value is defined by the same function q Fc — q Fc (<ft') as in Eqs. (jTf)) and but now 
with (ft 1 = (f>d rather than eft' — (ft where the angle (ftf. = <ftc(<ft) reads, 

#? = [0 - d?r/2 + ; d = ±l. (50) 

As confirmed in Section IV, for small hole concentrations < x <C 1 the hole-momentum absolute value q Fc is 
given by q Fc = y / xlr2. That expression remains a good approximation for x c \ < x < x c i and uo < U/At < u\, 

q F f rj \fxTr2e lj> d ; x c i < x < x C 2 ■ (51) 



5. The Fermi line for the hole- concentration range x 6 (a; c i,a;c2) 

For hole concentrations in the range x c \ < x < x C 2 and approximately uq < U/At < u\ the general si band momen- 
tum values q of states belonging to the reduced one-electron subspace are related to the corresponding momentum 
value of the m = generating ground state as given in Eq. (HOI) . The form of the matrix appearing in that 
equation is not in general known. The exception is for q at or near the si boundary line. Then it is given by A d x = A F 
where A F is the F rotation matrix (|44p . It follows from the value of the relative angle of the hole c Fermi momentum 
q F f and sl-boundary-line momentum <fj sl of the hole Fermi momentum k F of Eq. (|38p that the absolute value k F {4>) 
of the latter momentum reads, 

k h F (cf>) = ^[q'pM} 2 + {lBsi(4>)} 2 , (52) 
so that on using the integrals of Eqs. (l4"2l and (|4"5|) one finds, 

fj 7^ * 2 = fj^ * ( WfM] 2 + [tosiW] 2 ) = x 4tt 2 + (1 - x) 2tt 2 = (1 + x) 2n 2 . (53) 

Hence the Fermi line centered at —7? = — [tt, tt], which the hole momentum (|3"5|) belongs to, encloses the correct 
momentum area (1+x) 2ir 2 . Expression (|52[) is valid for the isotropic Fermi-line range x c \ < x < x C 2 and approximately 
Uo < U/4t < U\. For these ranges of x and U/At values the Fermi line has the simple relation to the c Fermi line 
and si boundary line described above such that the hole momentum qpf and momentum q^ sl of Eq. (j3"8)) are 
perpendicular. However, the general expression given in Eq. (|38[) is valid for the whole range < x < x* of the 
square-lattice quantum liquid provided that uq < U/At < yet for small hole concentration values we could not 
derive the angle between the directions where the si band momentum and c band hole momentum of the si fermion 
hole and c fermion, respectively, created under the one-electron excitation point. 

Moreover, one finds from the use of Eq. (|52|) that in the <j>p(4>) expression of Eq. (|43|) the quantity [k h ((j>) 2 — q(<j>) 2 ] 
is given by [k h (<j>) 2 — q(4>) 2 ] — q h {4 > ) 2 where q h (4>) is the absolute value of the hole momentum q h rj qpf at or near 
the c Fermi line of a c fermion created along with a si fermion of hole of momentum q ss <Tj sl upon creation of one 
electron of hole momentum k h w kp. As a result, for x c \ < x < x C 2 the F angle 4> F ((f>) defined in Eq. (|4^)l can be 
expressed as, 

m) = Cretan ( = rfarcsin ( jM ) = arccos ( ^) \ (g4) 

V9Bsl(0)y VV^cW] 2 + [iBsM)} 2 ) V \/ll h FM} 2 + [QBsim 2 ] 

Consistently with expression (l46l) for the angle cf)^, the matrix of Eq. (|44l) is such that for x c \ < x < x C 2 the 
momenta qssi of Eq. (|41|) and q^ sl of Eq. (l45]l are related as follows, 

QBsi = A F d q£ sl ; q£ sl = A d F q Bsl . (55) 

Here we have ignored contributions of order 1/N 2 . For x > and m — generating ground states and their two- 
electron excited states the si boundary- line momentum qssi is that given here and in Eq. (|41|) . In turn, for the 
states that span the reduced one-electron subspace it is instead the momentum qg sl of Eq. (|45|) and that of Eq. (|4"Tj) 
is called auxiliary momentum. 

The doublicity d — ±1 corresponds to two types of one-electron elementary processes that generate the same Fermi 
line. In such processes the hole Fermi momentum k F has an angle relative to the corresponding momentum —Qg s i or 
— QbsI ( ano - 9fc _1 or 9fc +1 ) gi yen by <j>p l or (ftp 1 = —(ftp 1 (and [(ftp 1 — n/2] or [(ftp 1 +tt/2}), respectively. It is confirmed 
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below in Subsection V-B that such two processes also reach the same Fermi energy but involve electrons with different 
velocities. Therefore, they refer to two degenerate ground states for the same electron numbers, which are obtained 
from the same m — generating ground state. Consistently, each value of the doublicity d = ±1 corresponds to a 
different reduced one-electron subspace. 

It follows that whether the F angle is given by (f>p^ or (ftp 1 = —(jf^ 1 depends only on the excited state and the 
reduced subspace it belongs to. Indeed and alike for the momentum shifts considered in the previous subsection, 
the value of the F angle depends for x c \ < x < x c i on the subspace transition doublicity and the shape of the si 
boundary line is different for one-electron excited states associated with the F angles given by (ftp 1 and (jip 1 = — ftp 1 , 
respectively. Note that since, 

[V] 2 ?B s i=fefi, (56) 

such shapes are transformed into each other under a rotation by the angle 2<j)p d . 

We now clarify the mechanism through which the c and si fermion description leads to the same Fermi line for 
one-electron addition and removal. In contrast to electron addition, electron removal does not conserve the values of 
the c band momenta of Eq. leading to a small overall finite momentum shift 5q® /N c . One can choose the hole 

concentration value of the initial m = generating ground state so that tf c ° = for the c band of that state subspace. 
In that case one finds that the hole Fermi momentum given in Eq. (|37[) can be expressed as, 

fe£ = - - <$ c d - qU = \-f c + q^ + qU , (57) 

where q c ° ^ refers to c band of the (N/2, N/2 — 1) ground state and we used that the excitation momentum of 
the one-electron removal process is — kp. For the reduced one-electron subspace the macroscopic momentum shift 
Sq^ = <f c ° and corresponding c band momentum shift 8q® /N c = <f c ° /N c of Eq. are such that the hole Fermi 

momenta given in Eqs. (1551) and (|57p associated with the one-electron addition and removal processes, respectively, 
are the same, 

5q c =q c = 2q sl = 2q Bs xW ; — = ^ ^ _ , . (58) 

The expressions given in Eqs. ([57| and (|58[) are valid for < x < x*. In turn, the particular form of the angle ^ of 
Eq. (I46|) refers only to the isotropic Fermi- velocity range x c \ < x < x c i. 

The inverse of the relation given in Eq. (|40l) plays an important role in the expressions of the si energy dispersion 
studied in Section IV. Neglecting the terms of order l/N% one finds, 

q Q = [A d sl ]~ l q; limf^]- 1 = I; d = ±1 , (59) 

where [A^-J -1 is the inverse matrix of A^. For qo at or near the si boundary line and x c i < x < x C 2 such a matrix 
reads [A^J - 1 = Ap d where A F is the orthogonal matrix given in Eq. (|44)) . 

As discussed above, the contribution to the si boundary line deformation of the two si band holes emerging in 
the case of subspace transitions from a m — ground state to two-electron excitations vanish so that the matrix Af t 
of Eq. (|4U|) reduces to the 2x2 unit matrix, A d x = I. In turn, for subspace transitions from such a ground state 
to excited states involving creation of one si fermion hole at the si boundary line that matrix is such that the si 
band hole momentum leads to the correct hole Fermi momentum kp through Eq. (|37|) . As also discussed above, the 
si boundary line is deformed under general subspace transitions from a generating m = ground state to excited 
states involving creation of an odd number of si fermion holes. In contrast, such a matrix reads A d x = I and the F 
angle vanishes for subspace transitions to excited states under which an even number of c fermions or c fermion holes 
and an even number of si fermion holes, respectively, are created. It also vanishes for subspace transitions to excited 
states involving creation of both an even number of c fermions or c fermion holes and of si fermion holes. 

According to the number and number deviations of Table IIVI of Appendix A, the charge and spin excitations are 
an example of subspace transitions under which none or two si fermion holes emerge so that the si boundary line 
remains unaltered and A d sl = I. Furthermore, also A^ = I for excitations generated by creation or annihilation of 
two electrons, as confirmed by the number and number deviations given in that table. 

Since one-electron (and two-electron) excitations do not couple to excited states with an even (and odd) number 
of si fermion holes, one has for m = generating ground states and their above two-electron excited states with 
N^ = 0, 2 holes in the si band that, 



gssi = qtsi = qBsi (4>) ; ?!> sl =0 + 7r; (for = 0, 2 states) . 



(60) 
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Here the angle <p s i = (f> + n is that of the auxiliary momenta belonging to the si boundary line. For x c \ < x < x C 2 it 
corresponds to the limiting value of the angle 4> d sl of Eq. ([46]) reached for tpp = 0. 

Finally, we emphasize that the doublicity d = ±1 introduced in this subsection is only well defined for the square- 
lattice quantum liquid in the reduced one-electron subspace. Furthermore, our above results and expressions refer to 
the N% — > oo limit whereas for a finite system the si boundary line deformations studied here cannot be described 
by simple si band rotations. 

IV. THE SQUARE-LATTICE QUANTUM LIQUID OF c AND si FERMIONS 
A. The general problem of expressing the model in terms of c and si fermion operators 

Expressing the Hamiltonian H of Eq. ([I]) in the one- and two-electron subspace in terms of c fermion and si fermion 
operators is a very complex problem. First one should write it in terms of rotated-electron creation and annihilation 
operators. Since for finite values of U/4t that Hamiltonian does not commute with the operator V associated with the 
electron - rotated-electron unitary transformation, its formal expression in terms of such operators is given in Eq. (U) 
and has an infinite number of terms. Fortunately, however, following the discussions of Subsection II-A only a small 
finite number of such terms is for U/4t > uq ~ 1.302 relevant to the physics of the square-lattice quantum liquid. 

Next one should invert the expressions given in Eqs. ([5J) and ([12")) to express the rotated-electron operators in 
terms of c fermion operators and the rotated quasi-spin operators, respectively. The obtained expressions refer to the 
LWS-subspace and read, 

4t = 4- (l +q h)+ ***** (} - 4) ; 4a = % (4- - e "' rl M ■ ( 61 ) 

The use of these expressions on the right-hand side of Eq. ((4]) leads to an expression of the Hamiltonian in terms of 
c fermion operators and rotated quasi-spin operators. 

An expression of that Hamiltonian for the subspace without rotated-electron doubly occupancy where the one- 
and two-electron subspace is contained is simply obtained by replacement in the Hamiltonian expression in terms of 
c fermion operators and rotated quasi-spin operators of the latter operators gi. by the spin operators si given in 
Eq. (flTT) where I = ±, z. Indeed, since the occupied sites of the c effective lattice correspond to the rotated-electron 
singly-occupied sites, the operator n?^ c = ft. c fpj,c appearing in that equation counts the number of the latter 
sites. Therefore, the spin operator s\*. — np j<c (fy. and the 77-spin operator pi = (1 — n? jtC ) q^. of Eq. (fTTj) are such 
that q l p . = si -(- pi and refer to rotated-electron singly occupied sites and rotated-electron doubly and unoccupied 
sites, respectively. Since for the subspace under consideration all Hamiltonian terms involving the 77-spin operator 
pi. = (1 — np. c ) gi. must be eliminated, replacement in the above Hamiltonian expression of the rotated quasi-spin 

operators gi = si + pi by the spinon operators si leads indeed to the desired expression. 

Hence expression of the Hamiltonian in terms of c fermion operators ft and ff. jC and spinon operators si where 
I = ±, z is a well-defined yet quite involved procedure. However, its ultimate expression in terms of c fermion and 
si fermion operators requires further expression of the contributions associated with its spin degrees of freedom in 
terms of two-spinon si bond-particle operators followed by the above extended Jordan- Wigner transformation of the 
latter operators. That problem is addressed in Refi^ for a particular class of important terms of the Hamiltonian 
rotated-electron operator expression which control the quantum fluctuations of the square-lattice quantum liquid 
and hence play a major role in its physics. Such Hamiltonian terms are consistent with the result of Rcf. 9 ' that the 
x > and m = ground states are spin-singlet states containing N s i — N/2 si bond particles. The two spinons 
of each of such spin-neutral composite objects refer to rotated electrons that singly occupied sites. Therefore, in the 
subspace without rotated-electron doubly occupancy where such ground states are contained and that the square- 
lattice quantum liquid refers to there is an energetic preference for the formation of spin-singlet rotated-electron bonds. 
Consistently with the number of independent bonds being N/2, the corresponding simplest Hamiltonian terms are of 
the form, 

N/2 

jjbonds = J2jjbonds + ( h . c .) . £ bonds = y A^v, [ct , >t gt 4 - ct , ^ gt >t ] , (62) 

j=l j' -const] 

where the complex gap function Aj/j» replaces the corresponding pairing operator, the summation Y^f j"[j-const] 
is over all N/2 spin-singlet rotated-electron bonds centered at a point of real-space coordinate [fy + fj"]/2 near 
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fj = [ff + fj"]/2 — I [a s /2] e Xd , and j — 1, N/2 — N s \. Here the two-site bond indices I ± 1 and d = 1, 2 arc those 
used within the notation of RefJ^ and a s = a/y/1 — x is the spacing of the spin effective lattice. 

After expression of (|62|) in terms of c and si fermion operators one finds in Ref.— that the action describing the 
quantum fluctuations of the square-lattice quantum liquid perturbed by small 3D anisotropy effects is controlled by 
contributions from the summation ^ j>> over nearest neighboring sites. However, in order to arrive to that result 
one must take into account the contributions from all possible bonds involving rotated electrons at arbitrarily distant 
sites associated with the summation j"\j-consi\ °f ^q. (|62|) . 

The expression of (|62j) in terms of c fermions and si fermions is performed in Ref. 28 *. The expression in terms of 
such operators of other Hamiltonian contributions is a very complex task, which we do not address in this paper. 
Fortunately, there is another path to handle the problem for some of these contributions, as described in the following. 
That allows to take into account the effects of the quantum fluctuations of the square- lattice quantum liquid associated 
with the short-range spin correlations, whereas those associated with the long-range superconducting order of the 
related extended quantum problem perturbed by small 3D anisotropy effects are addressed in Ref. 28 . 

For the study of the Hubbard model in the one- and two-electron subspace the quantum problem can be described 
by the Hamiltonian in normal order relative to the initial x > and m — ground state. Such a problem simplifies 
when expressed in terms of c and si fermion operators owing to the discrete momentum values q*j of such objects 
being good quantum numbers 9 . Indeed, that implies that the interactions between them have a residual character. 

B. The general first-order energy functional 

The general energy spectrum of the ground-state normal-ordered Hamiltionian of the Hubbard model on the square 
lattice in the one- and two-electron subspace describes the square-lattice quantum liquid of c and si fermions studied 
in this paper and in Refj2£. Here we study such an energy spectrum up to first order in the c and si fermion hole 
momentum distribution- function deviations for the U/At range uq < U/At < for which x* < 1/tt. That excludes 
the first order contributions found in Refi^ 8 - to emerge for the hole-concentration range x c < x < x* where x c rs 10~ 2 
from the superconducting fluctuations of the extended square-lattice quantum liquid perturbed by small 3D anisotropy 
effects. As justified below, it is convenient to consider the c and si fermion hole momentum distribution-function 
number operators rather than the related operators of Eq. (|26p . The former operators read, 

N' c l (q h ) = fo, c 4„ iC ; N^(q) = fl sl , (63) 

where for a c band centered at the momentum —7? the hole momentum values q h are given by q h = [q + vf]. Let 
N' c l (q h and N^q) denote the expectation values of such operators. We limit our study to the model in the one- 
and two-electron subspace defined in Ref£ for which the hole momentum values q h and momentum values q are 
good quantum numbers so that the corresponding hole momentum distributions N^(q h ) and N^q) are eigenvalues 
of the operators given in Eq. (f63l) . which read 1 and for unfilled and filled, respectively, hole momentum values 
and momentum values. For the model in such a subspace the energy dispersion e c (q h ) [and —e s i(q)} associated with 
creation of one c fermion (and si fermion hole) of hole momentum q h (and momentum q) onto the initial ground 
state is well defined. One of the goals of this paper is deriving such energy dispersions. 

For the subspace with a constant number N c — 2S C of c fermions, which the vacuum given in Eq. (|A1[) of Appendix 
A refers to, the value of the number t — [S c + S s ] of discrete momentum values of the si band provided in Eq. 
(|A17|) of that Appendix depends on the spin S s and thus corresponding number L s = 2S S of independent spinons. 
Therefore, for subspace transitions leading to finite deviations 5S S so that SN^ si = 5S S ^ the usual relation 
SNgi(q) — —5N s i(q) is not fulfilled. In turn, the number of discrete momentum values of the c fermion band is given 
by so that the usual relation 5N^(q h ) = —SN c (q h ) holds for the whole Hilbert space. 

Some of the reasons why it is convenient to express the energy functional associated with the ground-state normal- 
ordered Hamiltionian in terms of the c and si fermion hole momentum distribution- function deviations 5Nj}(q) and 
SNgi^q) rather than the usual momentum distribution- function deviations SN c (q) and 5N s i(q), respectively, are that 
for x > and m — ground states the si band is full, for their one- and two-electron excitations that band has one 
and none or two holes, respectively, and in addition the concentration of c fermion holes x in the c band equals that 
of holes. 

The ground-state normal-ordered c and si fermion hole momentum distribution- function deviations read, 

6Nt(q h ) = [N*(q h ) - N^{q h )] I 8N*@ = [N*@ - N^(q)\ , (64) 
where N^°{q h ) and N£ °(q) are the initial-ground-state values. 
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In the following we consider a range of hole concentrations x € (0, x*) where is the critical hole concentration 
found below above which within the square-lattice quantum liquid there is no short-range spin order at zero temper- 
ature. To first order in the c and si fermion hole momentum distribution-function deviations of Eq. (|64p the energy 
functional has the following general form, 

8E = -Y, e c (q h )5N h c {q h ) - Y, £sl(q)5N> s \(q) . (65) 
q h q 

Such an energy functional refers to the model in the one- and two-electron subspace for which the discrete hole 
momentum values q h and discrete momentum values q are good quantum numbers so that the first-order terms (|65|) 
and those found in Refi^ for an extended quantum problem perturbed by small 3D anisotropy effects refer to the 
dominant contributions. 

The corresponding excitation momentum spectrum is linear in the c and si fermion hole momentum distribution- 
function deviations and reads, 

5P = 5$ + J2 [* - tW*{q h ) - Y qsKi (0) > (66) 

q h q 

where 5q£ is the deviation in the c band momentum of Eq. (|34[) . Fortunately, the general excitation momentum given 
in Eq. (|66p does not involve the si band momentum deviation q~y of the si band momenta of Eq. (|34|) . The reason is 
that its role is to assure that the excited-state si momentum band obeys the sum-rule ^2^q = of Eq. (|3l)|) . Hence 
the net excitation momentum arising from the si momentum band is given by — X^f^^YsiCs) onr y- 

Another advantage of expressing the momentum functional (|66[) in terms of the si fermion hole momentum 
distribution- function deviation SN^ (q) rather than in the corresponding deviation 6N s i (q) is that the absence of 
the momentum deviation q^ from that momentum functional occurs only when one expresses it in terms of SN^q). 
That follows from the above exotic property according to which SN^q) ^ —SN s i(q) for state transitions leading to 
finite spin deviations SS S . 

C. The c band Fermi line and si band boundary line 

Based on several symmetries and approximations one can reach useful information about the c and si energy 
dispersions e c (q h ) and e s i(q), respectively, on the right-hand side of Eq. (|65l) . As a result of a symmetry such that 
the unitary operator V preserves nearest hopping only for rotated electrons, as occurs for electrons within the model 
([T]), the hole momentum (and momentum) dependence of the energy dispersion e c (q h ) (and e s \{q)) must for the model 
on a square lattice be fully determined by that of related elementary functions e c (q^ ) and e c (q X2 ) (and e s i(goa;i) and 
e s i(qox 2 )) whose arguments are the Cartesian components of the corresponding hole momentum q h — [q Xl ,q X2 ] (and 
auxiliary momentum <fo = [^fi] -1 ? = [Qoxn ?to 2 ])- The variables of the elementary functions e c (<jg ) and e c (gj 2 ) 
are the components of the hole momentum q h . Those of the elementary functions e s i(qo Xl ) and e s i(qo X2 ) are the 
components of the auxiliary momentum <fo of Eq. (|59j) of the m = generating ground-state associated with the 
momentum q under consideration. (For x > and m — ground states and their two-electron excited states one has 
that q = qo.) 

We often consider hole momenta q h at or near the c Fermi line and momenta q at or near the si boundary line 
such that = Aj? and [A^] -1 = Ap d where the 2x2 matrix A^ x and the 2 x 2 F rotation matrix Aj? are those of 
Eqs. (j4H]) and (|4"4"]) , respectively. We recall that for one-electron excited states the F angle (f>p(4>) vanishes at x = 
so that the corresponding si band rotation occurs for subspace transitions from a m = ground state to a reduced 
one-electron subspace provided that the hole concentration x of the initial ground state is finite. Furthermore, for 
subspace transitions to 7Vj\ = 0, 2 states the momentum values whose components appear in the argument of the 
si elementary functions e s i(qoxi) and e s i(qo X2 ) are not rotated since the F angle vanishes for such states so that 
qoxt — q Xl and qo X2 = q X2 . However, for the sake of generality we use the components qo Xl and qo X2 in the arguments 
of the elementary functions e s i(go^i) and e s i (qo X2 ), respectively, including when qo Xl = q Xl and qo X2 — q X2 . 

For the model on the square lattice the expression of the energy dispersions e c (q h ) and e s i(q) involves auxiliary 
energy dispersions e° c {q h ) and t° s i{q} given by, 

e c (<&) + e c (q X2 ) ~ eoG&Lj " e c (q h F d CX2 ) ; e° sl (q) = e°J ([A^q) = e°J(q ) , 
e s i(go^i) + e s i(<j(te 2 ) - e s i(q Bs ixi) - e s i(<j_B S i X2 ) . (67) 



e»i(9o) = 
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For one-electron excited states the Cartesian components [qBsix 1 ,qBsix 2 ] appearing here correspond to the auxiliary 
momentum qBsi generated from the momentum <fj sl belonging to the si band boundary line by the transformation 
([55]) . In turn, for a x > and m — ground state and its two-electron excited states such components refer to 
momenta qBsi = q B si belonging to the si boundary line. The shapes of the c Fermi line and si boundary line are 
fully determined by the form of the auxiliary energy dispersions defined in Eqs. (|67p as follows, 

qpc € hole c Fermi line <^=> e° c {q^) = , 

<fj sl £ si boundary line e° sl (qg sl ) = e°J(q Bsl ) = . (68) 

The c and si band discrete momentum values whose occupancy configurations generate the states belonging to 
the same V tower remain the same for the whole U/At > range. Since, as discussed below, for that range of U/At 
values the x = 0, fx = 0, and m — absolute ground state corresponds to a single V tower, its c and si band discrete 
momentum values remain the same for U /At > 0. Such a state refers to full c and si momentum bands with numbers 
N c = 2S C = N% and N s i = N% sl = N%/2 and thus N% = N' s \ = 0. These properties are consistent with an exact result 
according to which the momentum eigenvalues given in Eq. (|33[) are the same for all energy eigenstates belonging to 
a given V tower. 

On the contrary, the energy eigenvalues of states belonging to the same V tower are in general U /At dependent 
both for the model on the ID and square lattices. For the model on the square lattice and different U/At values 
the ground state of a canonical ensemble associated with a given hole concentration x > and spin density m = 
does not belong in general to the same V tower. This is consistent with the shape of the c fermion auxiliary energy 
dispersion e°(^- ) and thus that of the hole c Fermi line of Eq. (|6"8"|) depending on U/At, as confirmed below. 

The expression of the 7 = c, si fermion energy dispersions in terms of the elementary functions e 7 (g) on the right- 
hand side of Eq. (f6"T)) plays an important role in the square-lattice quantum liquid. In general we denote the ID 
variable of the arguments of the elementary functions e~ ( (q) and elementary velocities v 7 (q) — de 7 (q)/dq by q. For 
N% 3> 1 such elementary functions and the corresponding elementary velocities v 7 (q) have universal properties. The 
momentum P given in Eq. (|33p being additive in the 7 = c, si band momentum values q = [q Xl ,q X2 ] together with 
the general properties of the matrix of Eq. ([40]) reveals that due to symmetry the following relations hold, 

e 7 (g) = e 7 (-q) ; v 1 (q) = -Wy(-g) ; w 7 (0) = 0; sgn{u 7 (g)} = sgn{g} . (69) 

These symmetry relations imply corresponding symmetries for the c Fermi line and si boundary line of Eq. (|68p . For 
instance, the shape of the c Fermi line (and si boundary line) for the quadrant of the c band (and si band) whose 
hole momentum values q h = [Qa, j <7a; 2 ] ( an d auxiliary momentum values [^-J -1 q = qh = [qo Xl , qox 2 ]) are such that 
9xi — and gj 2 > (and qo Xl < and qo X2 < 0) provides full information about its shape in the remaining three 
quadrants. Therefore, we are free to restrict the expressions of momentum dependent quantities to a single quadrant. 

Among the si boundary-line momenta qg sl of the model on the square lattice it is useful to consider those whose 
corresponding auxiliary si boundary-line momenta qBsi related to them by the transformation Eq. (|55[) point in 
the nodal and anti-nodal directions. Indeed, by definition the nodal and anti-nodal si boundary-line momenta q^x 
and qg^ N , respectively, are those whose corresponding auxiliary momenta 9^ sl and have for instance for the 
quadrant such that q$ Xl < and qo X2 < the following Cartesian components, 
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Here <7^ sl and q^\ are the absolute values of both the auxiliary momenta q/g sl and qg^[ and corresponding momenta 
qg^[ and <?eji ', respectively. As confirmed below, for U/At > 0, m = 0, and x — > one has that, 

3b2=tt; q% s i=ir/V2; x^0. (71) 

The minus signs in the components of the expressions of the momenta q^ sl and qg^[ given in Eq. (|70[) refer to 
(f> = it/ A and <fr — 0, 7r/2, respectively. Indeed, the angle of such momenta is that provided in Eq. (pO)) . which reads 
4>si — 4> + tt. It corresponds to a limiting value of the angle (f)^ of Eq. (T^)) reached at 4>f = 0. 

As a result of the limiting values qp c — > plus those reported in Eq. (fTTj) . all reached for x — > 0, together with 
the si momentum band being full for x > and m — ground states, one finds that the maximum ranges of the 
momentum components in the argument of the elementary functions e e (q) and e s i(q') appearing in Eq. (|67|) are for 
rn = given by q G (— tt, tt) and q' G {—q B ^i,qBsi)i respectively Moreover, for m = and x > the zero-temperature 
chemical potential fi is fully determined by the charge c fermion energy dispersion and is given by the last two terms 
of the e° c {q h ) expression provided in Eq. (j67j), 

= i-eMLj ~ e c (q h F i 2 )} , (72) 
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where within our LWS representation we use the convention that for x > the chemical-potential sign is that of 
the hole concentration x. This general expression refers to x > 0. According to the results reported in Appendix A 
for x — and m = the range of the chemical potential is /i 6 (— /z , /j, ) where fi° = lim. E _i.o A 1 equals one half the 
Mott-Hubbard gap, whose magnitude is finite for U/At > 0. For < x < 1 and m = the chemical potential (|72|) 
is an increasing function of the hole concentration x such that /i° < fi(x) < fi 1 where /i 1 is given in Eq. (IA11[) of 
Appendix A and fiP has the limiting behaviors provided in Eq. (|A12[) of that Appendix. Hence fi° — > as U/At — > 
and fj, — > oo for U/At ^> 1. (This also applies to the ID model.) The energy scale [i° appears here associated 
with the charge degrees of freedom yet in the studies of Refi^ and below it is found to play an important role in the 
square-lattice quantum liquid spin degrees of freedom as well. 

The ground-state filled c fermion energy dispersion bandwidth W£ — [^c('Ifc) ~ e c(^)}i corresponding ground-state 
unfilled c fermion energy bandwidth W^, 1 — [ec(0) — ^(Qf c)]> ano - s l fermion auxiliary energy dispersion bandwidth 
W s i — [c®i{qBsi) ~ e Ssi(0)] can f° r the model on the square lattice be written in terms of the elementary functions 
as, 

Wf/2 = [e c (q h Fc )/2 - e c (n)} ; W*/2 = [e c (0) - e c (q h Fc )/2] ; W sl /2 = [e sl (±q% sl /V2) - e sl (0)] . (73) 
Moreover, the relations, 

e c {q h Fc ) = [e c (q h Fc /2) + At - W*/2] = [e c (0) - 8i + W%\ , (74) 

for the hole-like c band Fermi line where we have used that [W c p + W*] = St and, 

e sl (±q^) = [e. 1 (±qg al /V2) + W sl /2] = [e sl (0) + W sl ] , (75) 

for the si boundary line result from the relations of Eq. (|68p . which define the c Fermi line and si boundary line, 
respectively. 

According to Eq. (|A13j) of Appendix A, the ranges of the energy e c for addition onto a x > and m — ground 
state of one c fermion and the energy — e c for removal from that state of one c fermion are controlled by the energy 
bandwidths W' c l and W?, respectively, of Eq. (73J). Here W' c l = [ADt - W?\ where D = 1, 2. According to Eq. 
of Appendix A such ranges are defined by the inequalities < e c < Wj? and < — e c < W£, respectively, for both 
the model on the square and ID lattice. Moreover, one has that Wg = ADt for x — and W£ = for x = 1. Indeed, 
for U/At > the energy bandwidth Wj? & (0,ADt) (and W% € (0,ADt)) decreases (and increases) monotonously for 
increasing values of hole concentration x € (0, 1). 

According to Eq. (|A14[) of Appendix A, the range of the energy — e s ± for adding to a m = ground state one 
si fermion hole is for the square- lattice model controlled by the si fermion auxiliary energy dispersion bandwidth 
W s i or maximum si fermion pairing energy per spinon |A| according to the inequality < — e s \ < max{W s i, |A|}. 
The limiting value = lim^^o W s \ of that energy bandwidth plays an important role in the physics of the square- 
lattice quantum liquid. According to Eq. (|A15|) of Appendix A, upon increasing the value of U /At it decreases from 
W°i = 2Dt at U/At — to w 2Dwt 2 /U for U/At -t oo. In turn, as discussed in Ref^, the si fermion pairing 
energy per spinon |A| has a singular behavior at x — and for < x <C 1, being given by fjP/2 and Ao = \im x ^.Q |A|, 
respectively, where Ao < /io/2 for U/At > 0. Such a behavior is due to the sharp quantum phase transition occurring 
at x — 0, such that there is long-range antiferromagnetic order and short-range spiral-incommensurate spin order at 
x = and for < x <C 1, respectively. 

D. The c momentum band and c fermion energy dispersion 

As discussed in the previous section, for the model on the square lattice and U/At > the c band hole momentum 
values q h are good quantum numbers. The energy e c (q h ) [and — e. c {q h )\ associated with addition to (and removal 
from) the ground state of one c fermion of hole momentum q h has for U/At > 0, x > 0, m = 0, and both the model 
on the square and ID lattices the general form, 

e c {q h ) = sgn(e°(q"))/ + e° c (q h ) . (76) 

Here £° c {q h ) is the c fermion auxiliary energy dispersion given in Eq. (167[) and the x = chemical-potential zero level 
corresponds to the middle of the Mott-Hubbard gap whose magnitude is 2/i°. (We recall that its limiting behaviors 
are given in Eq. (|A12[) of Appendix A.) For U/At <C 1 and U/At > 1 they are such that 2^° -> and 2fi Q oo, 
respectively, both for the model on the square and ID lattices. 
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It follows from Eq. (|75|) that e c (q h ) — t° c {q h ) for x > and m = where consistently with the inequalities provided 
in Eq. (|A13[) of Appendix A, 

e° e (0) = [4Dt-W*\; D = l,2; e° c (q^) = 0; = -W* . (77) 

For U/At 1 the electrons that singly occupy sites are insensitive to the on-site repulsion and behave like the 
spinless c fermions. It follows that the corresponding c fermion occupancy configurations that generate the charge 
degrees of freedom of the energy eigenstates have a non-interacting character. Therefore, in such a limit for which 
the rotated electrons become electrons the c band hole momenta q h are associated with electron hopping and the 
corresponding c fermion hole-momentum occupancy configurations are behind the whole kinetic energy. One then 
finds that in such a limit the elementary function e c (q) has an non-interacting form and reads e c (q) = —U/2 + 2tcosq. 
From its use in the general expression provided in Eq. (|67[) one finds that the c fermion auxiliary energy dispersion 
e®{q h ) defined in that equation is for U/At^> 1 given by, 

D 

e° c (q h ) = 2tJ2 [cosfe* ) - cos(^)] ; U/At » 1 . (78) 

i=l 

(For ID one has that q^ l = q h and q FcXl = q Fc -) For small values of Ig^l 2 — \qp£ | 2 (and (q h ) 2 — {q Fc ) 2 ) and < i « 1 
this gives, 

«,(fV-^M, (79) 

for the model on the square lattice and U/At > 1 (and e c (q h ) = -[(q h ) 2 - {q Fc ) 2 ]/2mf for ID and U/At > 1) where, 

m~ = l/2t. (80) 

For the limit U/At — > oo, which expressions (|79|) and (|80[) refer to, the c fermions are invariant under the electron - 
rotated-electron unitary transformation. Then they are the non-interacting spinless fermions that describe the charge 
degrees of freedom of the electrons. Also for maximum spin density m — (1 — x) and arbitrary U /At values the c 
fermions are invariant under that transformation and describe the charge degrees of freedom of the electrons of the 
fully-polarized state (|Al|) of Appendix A. In that limit their energy dispersion is for x > of the form given in Eq. 
(|A7|) of that Appendix. Transformation of the momenta of that dispersion into hole momenta gives expression (|78|l . 
For m = (1 — x) and arbitrary V j At values the c Fermi line equals the Fermi line. 

For spin density m = the expression (|79|) can be generalized to finite values of U /At by the use of the relation 
between the U/At 1 and U/At > physics. That relation is controlled by the electron - rotated-electron unitary 
transformation associated with the operator V. For < x <C 1 and U/At > the unitary character of such a 
transformation preserves for the model on the square (and ID) lattice the value of the hole momentum qpf of Eq. 
({51) (and hole momentum values ±q Fc — ±xtt) and that of e c (q h ) (and e c (q h )) with the bare mass to£° = l/2t 
replaced by a renormalized mass m*, 

, c(f V-^Mf, (81) 

2m* 

(and e c (q h ) = — [(q h ) 2 — (qp c ) 2 ]/2m*). From the form of the c fermion energy dispersion (|8ip one confirms that for 
the model on the square lattice, < x <C 1, and U/At > the absolute value qp c ((f>) of the c hole Fermi momentum 
<f^ c d of Eq. is indeed independent of the angle (p and given by q Fc {4>) — ^/xtt2, as provided in Eq. (|5"T|) . It follows 
that \qpc\ 2 = xAir and for U/At > and < i <C 1 the hole momentum values that belong to the c Fermi line of the 
model on the square lattice have indeed the form given in Eq. (|51[) . The expression of the c fermion energy dispersion 
provided in Eq. ([8lj) is valid for some finite energy bandwidth |e c (<f' t )| measured from the c Fermi level. For the 
model on the square lattice at intermediate U/At values U/At £ (u ,ui) and < x < x* that expression is a good 
approximation for |e c (<f' l )| < W^\ x =x, where Wj? = [8t — Wg] (and W/!) is the energy bandwidth of the ground-state 
unfilled (and filled) c fermion sea. 

Moreover, from the use of Eqs. (|72|) . (l76l) . and ([81]) one finds that for U/At > and hole concentrations < i « 1 
the zero-temperature chemical potential dependence on x is given by, 
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where /i° is the chemical potential /jP = lim^^o M of Eq. (|A12[) of Appendix A. For intermediate U /At values 
U/At € (uo,Ui) the expressions given in Eq. (|82|) remain good approximations for < x < x*. 

Alike the quasiparticle mass ratios in Fermi liquid theor y 39 ! 40 , the £//4£-dependent charge mass ratio r c — m™/m* 
involving the c fermion mass appearing in the expression provided in Eq. (|81|) and the spin ratio r s = Aq/AW^ 
considered below play an important role in the quantum liquid physics. However, in contrast to the former liquid here 
the mass = l/2i of Eq. (|80[) . which refers to the limit of infinite on-site interaction plays the role of that of the 
non-interacting Fermi liquid. Moreover, in contrast to the Fermi-liquid quasi-particles the c fermions do not evolve 
into electrons upon turning off adiabatically the interaction U. Instead, upon adiabatically turning off the parameter 
At 2 /U they evolve into the spinless fermions that describe the charge degrees of freedom of the electrons that singly 
occupy sites within the energy eigenstate configurations. In that limit all finite-energy electronic configurations that 
generate the energy eigenstates involve such electrons only. In turn, the to = ground state si fermion occupancy 
configurations that generate the spin degrees of freedom of such electrons are degenerate and for the ID model become 
in that limit those of the spins of the spin-charge factorized wave function introduced independently by Woynarovich 37 
and by Ogata and Shiba 3 ^, respectively. For the model on a square lattice such si fermion occupancy configurations 
are within the suitable mean-field approximation (|2"Tj)) for the fictitious magnetic field B s i of Eq. (|2"Tj) those of a full 
lowest Landau level with N s i = N% = N/2 one-sl-fermion degenerate states of the 2D QHE. 

For the latter model the ratio to£°/to* increases smoothly for increasing values of U/At and has the following 
limiting behaviors, 

to 00 

r c = — c - = , U/At -> 

777* 

ps 2 c- 1 « 0.736, f7/4i« u ps 1.302, 
= 0.27tt w 0.848 , U/At ps = 1.525 , 

= 1, U/tt^oo. (83) 

This charge mass ratio and the si fermion energy dispersion given below for x = and to = are well-defined for 
U/At > 0. In turn, the x dependence of the energy parameter | A| appearing in si fermion energy dispersion for x > 
and 777 = derived below is a good approximation for u Q < U/At < u\ and thus the studies of the remaining of this 
paper are often restricted either to such a U/At range or to U/At > u . 

In contrast to the quasiparticles of a Fermi liquid, for the c and si fermions of the square-lattice quantum liquid the 
non- interacting U/At — > limit is not trivial. The values obtained in that limit for some physical quantities are well 
defined. However, in some cases one must be careful, since the physics associated with the U/At — > limit and that of 
the U / At = non-interacting system may not coincide. For instance, for the model on the square lattice the expression 
(|8ip of the c fermion energy dispersion for small \q h \ and < x <C 1 and the corresponding chemical-potential formula 
([8"2]l are valid for finite values of U/At > but do not apply at U/At = 0. Indeed, they apply provided that 1/m* 
is finite at x — so that expressions (|5Tj) and (1521 arc the dominant first-order terms in the hole concentration x of 
expansions in x valid for < x -C 1. However, for U/At — > one has that 1/m* = 0. 

Study of the spectrum of the model (JXJ) on the square lattice at U /At = reveals that for hole concentrations x«l 
the chemical potential and its derivative with respect to x are given by fi(x) ps x 2tt/w(x) and d/j,(x)/dx ps 2n/[w(x) — 
2/(ni)], respectively, where the inverse function of w(x) reads x(w) — [8wt/iv] e~~Tr h 2 so that \im x ^ow(x) = 00. 
Since for U/At — > one has /i° — > 0, the U/At = chemical-potential expression (i(x) ps x2ir/w(x) is to be compared 
with the general expression (i(x) ps x2ir/m* c of Eq. (|82|) for U/At — » 0. From such a comparison one finds that the 
U/At — quantity w{x) plays the role of the effective mass 777*. For < x <C 1 and U/At = the x dependence 
of w{x) is defined by the equation w{x) = x[n /8t 1 \e w ^ t ^~^ so that it is large but finite for small values of x. It 
becomes the mass to* in the limit x — > for which one finds that \im x -+ow(x) = 00, consistently with the limiting 
behavior r c = m™/m* c ^ as U/At ->• of Eq. {S3J where according to Eq. dSDJ) = 1/2*. The result 

lim[//4 t _yo to* = lxm x -^o w(x) = 00 is consistent with the U/At — > limiting behavior of the ratio r c given in Eq. 
([83]). For < x -C 1 and U/At = the U/At > chemical-potential formula ([82]) must be replaced though by 
n{x) ps x2tt/w(x), where rather than w(0) = to* = 00 one must use the large but finite value of w(x) corresponding 
to the small hole concentration x under consideration. This confirms how careful one must be concerning the interplay 
of the physical quantities magnitudes at U /At = and for U/4t —> 0, respectively. 

For the model on the ID lattice the mass to* is that called m* h in Ref. 49 . Its exact expression is given in Eq. (48) 
of that reference and reads to* = Cg/2cii for all values of U/At so that, 

r c ^^-^; Cj = l-2 ckj^L-, j = 0,l; D = l, (84) 
m c c o Jo 1 + e 2* 

where Jj(uj) are Bessel functions. In ID the charge mass ratio l/r c = rn* c /m^° has a behavior opposite to that of the 
model on the square lattice: For the square (and ID) lattice it decreases (and increases) smoothly upon increasing 



35 



U/At from m*/m^° = oo (and m* c /mf = 0) for U/At ->• to m*/m£° = 1 as U/4* » 1. Such a qualitatively different 
behavior is related to the different value of the chemical-potential derivative d/j.(x)/dx\ x =o at U/At = 0. For J7/4t = 
one has that fi(0) — both for the model on the square and ID lattices. In contrast to ID, where the derivative 
dfj l (x)/dx\ x —o — irt is finite, for the model on the square lattice it vanishes, dfi(x)/ 'dx\ x —o — 27r/[w(0) — 2/(irb)\ = 0. 

E. The si momentum band and si fermion energy dispersion at x = and m — 

1. Invariance under the electron - rotated- electron unitary transformation of the = and m — absolute ground state 

The expression kp = [qp^ — QBsl] given in Eq. (|38| for the hole Fermi momentum kp of Eq. (|37|) refers to electron 
addition and is valid for U/At > 0. (That the directions of and q^ sl are perpendicular is a result that applies to 
x € (x c i,x C 2) and approximately U/At <G (uq,u\) and otherwise the relative directions of such momenta remains an 
unsolved problem.) Furthermore, according to Eq. (|57)) . for one-electron removal one has kp = [—q£ + qp^ + qg sl ] 
so that <f c ° = 2(f/ 1; as given in Eq. (|55|) . 

At x — the i* 1 angle provided in Eq. vanishes for subspace transitions from the m = ground state to 

one-electron excited states so that <fj sl = q Bs \ in Eq. (|55|) . Moreover, expression (j3"5j) for the Fermi hole momentum 
kp simplifies at x — and m = since then qp c — so that, 

kp=-qBsi] k F = {-n - q Bsl ] , x = 0, m = 0, (85) 

where the &f expression follows from Eq. (|37[) . 

Importantly, the x = 0, /j, = 0, and m = absolute ground state is the only ground state that for U/At > 
belongs to a single V tower. Hence for its subspace the si boundary-line momenta q Bs i are independent of U/At. 
This is consistent with the above x = and m = expressions given in Eq. ([55} and the requirement of the particle- 
hole symmetry specific to x = 0, fi = 0, and m = that fc^ = kp. The latter requirement is fulfilled provided 
that for the absolute ground state the si band coincides with an antiferromagnetic reduced Brillouin zone such that 
I fell + 1 9xa I ^ i") which is enclosed by a boundary line whose momenta <fs s i have Cartesian components qBsix! and 
qBsix 2 obeying the equations, 

qBsl G si boundary line gssij;! ± ^si^ = tt or q Bs i Xl ± fesi^ = -7f . (86) 

This implies that the absolute value q Bs i{4>) has the form given in Eq. (fi2"|) for x <C 1. Hence the si boundary line 
refers to the lines connecting [±7r, 0] and [0,±7r]. Such a si band shape and area are consistent with the result of 
Refi£ according to which at x = and m — the square si effective lattice has spacing a s i = y2a. Its periodicity 
increases relative to that of the original lattice owing to the appearance of the long-range antiferromagnetic order, 
consistently with the results of that reference. 

This result reveals that not only the boundary-line momenta q Bs \ are the same for the whole range of U /At > 
values but also in the U/At — > limit are associated with the correct U/At = Fermi line. It follows that the shape 
of that line is for the Hubbard model on a square lattice independent of U/At and equals that of its U /At — Fermi 
line, consistently with the absolute ground state corresponding to a single V tower of energy eigenstates. 

Only for x = and m — is the Fermi line invariant under the electron - rotated-electron transformation: its 
Fermi momentum values have then the form kp = kp = —qBsl where q Bs i is U/At independent and has Cartesian 
components obeying Eq. (|86[) . Our results then provide strong evidence that the x = 0, /i = 0, and m = absolute 
ground state with numbers Nj; — 0, N? — N% — 2N s i, and 2S C = N% is invariant under the electron - rotated- 
electron transformation. The theory vacuum of Eq. (jAll) corresponds to a fully polarized state and is indeed invariant 
under that transformation. 

The form of the x = and m = Fermi line found here together with the property that for U/At > the c and si 
bands of the absolute ground state are full is consistent with the discussions of Refj^ according to which the Hubbard 
model on a square lattice is for x = 0, fx = 0, and m = a Mott-Hubbard insulator with long-range antiferromagnetic 
order for all finite values of U /At. This is for instance consistent with the behavior reported for the Mott-Hubbard gap 
2/jP in Eq. (|A12I) of Appendix A, such that it is finite for finite values of U/At and behaves as 2fi° — > for U /At —> 
and 2/i° — > oo for U/At 3> 1. The studies of Ref^£ on the quantum fluctuations of the present quantum problem 
confirm the occurrence of zero-temperature long-range antiferromagnetic order for x = 0, m — 0, and all finite values 
of U/At. 

Interestingly, our results reveal that the singular physics associated with the x = 0, fi = 0, and m — absolute 
ground state and the corresponding Mott-Hubbard insulator phase is related to that state being invariant under the 
electron - rotated-electron unitary transformation. However, since for finite values of U /At the Hamiltonian and the 
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generator of the global U(l) symmetry are not invariant under that transformation the physics of the Mott-Hubbard 
insulator phase is not trivial. 

2. The si fermion energy dispersion for x — and m = 

The expression of the si fermion energy dispersion e s i (q) involves the auxiliary energy dispersion e®i(q) of Eq. (|6"T)) 
whose momentum dependence defines the shape of the si band boundary line through Eq. (|68|) . Since for m = 
ground states the si band is full it follows from Eq. ([731 that [e°x(?Bsi) — e si(0)] = W*i f° r x — ^ and m = 0. However, 
the dispersion e®i(q) does not include the si fermion spinon-pairing energy. According to Eq. (|68l) . ejx(9B*i) = so 
that the energy that is left for momentum values belonging to the si boundary line is the si fermion pairing energy 
per spinon. Hence for x = and m — the si fermion energy dispersion e s i(o) reads for momenta at the boundary 
line, 

e s i(q Bs x) = -\A sl {q Bsl )\ ; < \A sl {q Bsl )\ < M °/2 . (87) 

Here |A s i(<jb s i)| denotes the pairing energy per spinon of a si fermion of momentum qssi, which may be different 
for different si boundary-line momentum directions and may vanish for some of these directions so that it obeys the 
above inequality. 

The maximum si fermion spinon-pairing energy equals the excitation energy below which the long-range anti- 
ferromagnetic order survives for x = 0, m = 0, and zero temperature T = 0. In reference^ evidence is provided that 
the Mermin and Wagner Theorem^ applies to the half- filling Hubbard model on a square lattice for U/At > 0. This 
is consistent with the existence of long-range antiferromagnetic order only at x = 0, rn = 0, T = 0, and U/At > 0, in 
agreement with the numerical results of Refs i 51 ' 52 , and its replacement by a short-range spin order both for (i) x = 0, 
T > 0, and U/At > and (ii) < x < 1, T > 0, and U/At > 0. 

The electron occupancy configurations associated with the long-range antiferromagnetic order refer to a subspace 
that for U /At finite exists below the energy scale fjP, which equals one half the magnitude of the Mott-Hubbard gap. 
Indeed, the Mott-Hubbard gap 2fi°, which refers to the charge degrees of freedom, also affects the spin degrees of 
freedom. In reference^ it is argued that for the model on the square lattice the occurrence of long-range antiferro- 
magnetic order requires the absence of a 77-spin effective lattice and the spin effective lattice being identical to the 
original lattice. That happens only for a subspace with numbers N% = and N% = N% so that indeed there is no 
?7-spin effective lattice and the spin lattice has as many sites as the original lattice and is identical to it. In such a 
subspace there are neither rotated-electron doubly occupied nor unoccupied sites, so that the Hubbard-model spin 
effective lattice is identical to that of the spin- 1/2 isotropic Heinsenberg model on a square lattice. 

For U/4t^> 1 the spin degrees of freedom of the half- filling Hubbard model are described by the Heisenberg model 
and the above subspace for which N% = N% describes the whole finite-energy physics, since /i° — U/2 — > 00. However, 
for finite U/At the spin occupancy configurations behind the long-range antiferromagnetic order that refer to a spin 
effective lattice with N% — N% sites exist only for excitation energy below fjP. This justifies why for the Hubbard 
model on a square lattice /z° is the energy below which the long-range antiferromagnetic order survives for x = 0, 
/x = 0, m = 0, and temperature T = 0. 

In the absence of any spin order the si fermion spinon-pairing energy would vanish and the energy dispersion e s i(q) 
of the si fermions would be given by the auxiliary energy dispersion e®i(q) of Eq. (|67j) . An important point for the 
derivation of the si fermion energy dispersion e s \{q) in the presence of such an order is that for excitation energy 
below the energy scale fjP one can ignore the amplitude fluctuations of the order parameter and the problem can be 
handled by a suitable mean-field theory where the occurrence of that order is described for the x = and m — 
problem by a si energy dispersion of the general form — \/l e si('?)| 2 + l^si(<7)| 2 where |A s i(g)| is the si fermion pairing 
energy per spinon appearing in Eq. (|87p . Indeed, the si boundary- line condition of that equation imposes that the 
gap |A s i(g)| is the pairing energy per spinon of a si fermion of momentum q. 

There are both bad and good news. The bad news is that according to Eq. (|A12[) of Appendix A for small U/At 

values the energy scale /i°/2 becomes small and given by /i°/2 s» 16 1 e^V 4 */^ so that the amplitude fluctuations 
of the order parameter cannot be ignored and thus the energy dispersion — e^i (<?) | 2 + |A s i(<t)| 2 may not be a 
good approximation. Indeed and as discussed above, within the c and si fermion description the smal\-U/At physics 
corresponds to a non-trivial problem. This applies to the derivation of the si fermion dispersion e s ±(q) for small 
values of U/At, which remains an unsolved issue. Fortunately, the above energy dispersion is expected to be a good 
approximation for U/At > uq. The good news is then that according to the studies of Subsection VI-B the relation 
of the present quantum problem to the unusual physics of the Mott-Hubbard insulator parent compounds such as 
LCO refers to a value U/At s» u* = 1.525 > uq so that the description introduced here is of interest for the study 
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of such materials. The U/At value uo corresponds to the maximum magnitude of the energy scale 2Ao whose U/At 
dependence is investigated in RefA 

We recall that for transitions from the x = and m = ground state to all subspaces of the corresponding one- 
and two-electron subspace the matrix equals the 2x2 unit matrix I so that the elementary functions can be 
expressed as e s i(qo xl ) — esifei) an d e s i(qo X2 ) = e s ±(q X2 ). The same symmetry arguments that imply that the 
si fermion energy dispersion e° s i(q) has for the model on the square lattice the general form provided in Eq. (|67[) 
impose that |A s i(<f)| is a superposition of the elementary functions e s i(q xi ) and e s i(q X2 ). Since a function of the form 
|A s i(g)| oc [a*°/2] \e s i(q xi ) + e s i(q X2 )\ would lead to a mere change in the si energy bandwidth, the only remaining 
possibility for the model on the square lattice is, 

\^si(q)\ oc — |e sX (g Xl ) - e s i(q X2 )\ , (88) 

where /i°/2 is the maximum si fermion pairing energy per spinon. 

That for the x = 0, fj, = 0, and m = absolute ground state the components of the boundary-line momenta qssi 
obey the relation given in Eq. ([86)) and for all values of U /At that line coincides with the limits of an antiferromagnetic 
reduced Brillouin zone is consistent with a elementary function e s i(q) ss — cos q whose U /At dependence occurs 

through that of the energy bandwidth W^. Hence combining all above results one finds, as confirmed in Subsection 
VI-B, for approximately U /At > uo, excitation energy ui below and temperatures T below ui //cb, that the following 
si fermion energy dispersion is a good approximation for x — and m = 0, 

/ W° 

tsl(q) = -VI&(S)l a + |A.l($)l a i e° sl (q) = -^[cos qxi +cosq X2 }, 

\^si(q)\ = -^F sl (q); F sl (q) = . (89) 

Here the maximum gap magnitude ^°/2 of the si fermion energy dispersion equals the maximum si fermion pairing 
energy per spinon and corresponds to the si band momentum values belonging to the si boundary line and pointing 
in the anti-nodal directions, 

/i /2 = |A 5l (^)| = |A sl (7r,0)| = |A sl (0 j7 r)|; x = 0, m = 0. (90) 

According to Eq. (|AT2]) of Appendix A, fjP/2 sa 16 te-^V^/u for jjj^ < 1 w h er eas fi°/2 w [U/A-2t] for U/At > 1. 
However, only for approximately U/At > uq is is the expression for the energy dispersion e s \(q) provided in Eq. (|89j) 
valid and thus can one identify the si fermion anti-nodal energy gap |A s i(gg^)| with /x°/2. Indeed, owing to the 
effects of the amplitude fluctuations of the order parameter, for small U /At values such a si fermion energy dispersion 
is expected to be different from that provided in Eq. (|89|) . 

We recall that for excited states of the x = 0, n = 0, and m = absolute ground state the one- and two-electron 
subspace considered here corresponds to excitation energy to < . For excitation energy below such an energy scale 
there is both no rotated-electron and rotated-hole double occupancy. Indeed, for the x = 0, /x = 0, and m = 
absolute ground state with the chemical-potential zero level at the middle of the Mott-Hubbard gap the energy scale 
/i is the smallest energy required for creation of both one rotated-electron doubly occupied site and one rotated-holc 
doubly occupied site. Hence that the si energy dispersion is valid for excitation energy below reveals that 
provided that U/At > uo it is fully consistent with our general c and si fermion description, which refers to the one- 
and two-electron subspace. 

Moreover, that according to Eq. (I89|) the si fermion spinon-pairing energy is at x — 0, m — 0, and U/At > uq 
given by 2|A s i(<f)| = \p /2] | cosg Xl — cos<? X2 | reveals that the si fermion spinon pairing has <i-wave symmetry. An 
important point is that at half filling one may have for the model on the square lattice a ground state with both 
long-range antiferromagnetic order and si fermion spinon pairing with <i-wave symmetry. According to RefJ^, for 
values of the two-site bond length £ g given in that reference not too small the absolute value \h g \ of the coefficients 
h g of the si bond-particle operators of Eq. (|A18[) of Appendix A falls off as \h g \ C (£ 3 )~" sl ■ Provided that for 
the x = and m = absolute ground state the value of the exponent a s \ is approximately in the range a s \ < 5, 
the corresponding spin-singlet spinon si bond pairing of the N%/2 two-spinon si bond particles associated with such 
operators can for U /At > uq have <i-wave symmetry yet the corresponding spin occupancy configurations have long- 
range antiferromagnetic order 53 . The upper value 5 is that obtained from numerical results on spin-singlet two-spin 
bondsi^i^ 5 .!^ 3 -. 

The construction of the si bond-particle spin configurations fulfilled in ReL 1 - takes into account the result found 
here that for x = and m = the si fermion energy dispersion has for U/At > u n the form 2|A s i(<f)| = \fi°/2] | cosq^ — 
cosga; 2 | given in Eq. (|59")l so that the si fermion spinon pairing has ci-wave symmetry. Let us denote by sl the 
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primary bond operator obtained by restricting the summation in the expression for the si bond particle operator 
gp , s i of Eq. (|A18|) of Appendix A to g — only. The corresponding si fermion operator is according to Eq. (|T6]) 



— P -i<t>j 



gpj.si where the phase 4>j lS i is provided in Eq. (TTTl) . By combining Eq. (|A18|) of Appendix 



given by fp jtSl = e 

A with Eq. (|16l) one can then introduce a primary si fermion operator /2 sl 
' tic results of Ref. 10 that for the model on the square lattice such a primary sl fermion operator reads, 



-i<j>j 



One then finds from 



d=i i=±i 

2 



V2 



d=l l = ±l 



2 + S rj+2?li 



(91) 



(Rather than doublicity d = ±1, here d = 1 and d — 2 refer according to the notation of Ref.— to two types of 
bonds.) The spin operators of the second expression given here are those of Eqs. (jlll) and p^|) -([T5 |) and the index 
values d = 1 and d = 2 refer to the horizontal and vertical primary two-site one-link bonds, respectively, with the 
d = 1 (and d — 2) two-site bond whose index I reads I = — 1 and / = +1 being that whose link connects the site of 
real-space coordinate fj with the nearest-neighboring site on its left-hand and right-hand side (and above and below 
it), respectively. The summation and phase factor J2d=i J2i=±i(~ °f Eq. (|9"Tj) then imply that for the square 
lattice the primary sl fermion operator /2 sl involves a spinon pairing with cZ-wave symmetry, consistently with the 
corresponding sl fermion energy dispersion having the form e s i(q) 



2|A s i(g)| = [n°/2] | cosq Xl - cosq x 



- Vl^i^l 2 + \&si(q)\ 2 given in Eq. where 



F. The sl momentum band and sl fermion energy dispersion for x > and m = 

1. The general sl fermion energy dispersion for x > and m = 

Alike for x = 0, for < x < a;*, m — 0, and approximately U /At > ug the auxiliary energy dispersion does not 

contain the sl fermion pairing energy per spinon associated with the short-range spin order. Again since according to 
Eq. (|68l) the corresponding auxiliary energy dispersion vanishes at the sl boundary line £®i(<1b s i) = ^i(qBsi) = 0, 
the energy this is left for momentum values belonging to the sl boundary line is the sl fermion pairing energy per 
spinon so that the following result holds for x > and m = 0, 

e sl (q£ sl ) = -\A sl (q£ sl )\ = -\Al(q Bsl )\ ; \A sl (q)\ = lA^K]- 1 ^ = |AjJ x (gb)| ; < |A^(g Bsl )| < |A| . (92) 

Alike in Eq. (1871) . here |A s i((^)| = (A^go)] is the pairing energy per spinon of a sl fermion of momentum q. 

As discussed for the x = and m = problem and except that for x > the F angle is finite for transitions from 
a x > and m = ground state to the reduced one-electron subspace yet vanishes for transitions to two-electron 
excited states, in the absence of spin order our analysis in terms of the elementary function e s i(go) would lead to a 

sl fermion energy dispersion of the general form e s i(q) — e° sl (q) = e^"(q*o) = ^si^ox^ + s s i(qo X2 ) + constant. The 
amplitude of the order parameter of the short-range spin correlations is frozen below the energy 2|A|. Hence the 
occurrence of the corresponding short-range spin order leads for m = initial ground states corresponding to the 
hole concentration range < x < x* and below an excitation energy and a temperature for which the amplitude 
fluctuations of the corresponding order parameter can be ignored to an effective mean field theory associated with a 
sl fermion energy dispersion of the general form, 

e sl (q) = - ^i(<z)l 2 + |A s i(<?)| 2 , (93) 

where e^^q) is the sl fermion auxiliary energy dispersion given in Eq. (|67[) and the gap function |A s i(g)| is derived 
below. 

Similarly to /i°/2, also the energy scale Ao becomes very small and reads Ao ~ /io/2 ~ 16t e~ 7r "^ 4t/ ' c/ for small 
values of U /At so that then the amplitude fluctuations of the order parameter cannot be ignored and the sl fermion 
energy dispersion is not expected to have the form given in Eq. (I93p . Since according to the investigations of Ref. 9 
the energy scale A also vanishes for If/At ^> 1, for hole concentrations < x < x* the range of validity of the sl 
fermion energy dispersion expression (|93[) corresponds approximately to U/At > u$. Indeed, all spin energy scales 
vanish for U/At 3> 1 so that Aq vanishing in that limit has less severe consequences than it being small for small 
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values of U/At. Consistently, the si fermion energy dispersion provided in Eq. remains valid provided that its 
energy bandwidth is finite. A qualitatively different 2D QHE physics is reached in the limit U /At 3> 1 when the 
spectrum of the two-spinon si fermions becomes dispersionless and all one si fermion states are degenerate. Provided 
that the energy scale At 2 /U is finite that energy bandwidth is also finite and the above description associated with 
the si fermion energy dispersion of Eq. applies. 

Hence provided that the value of U/At is larger than uq ss 1.302 and finite the si fermion energy dispersion is of 
the form given in Eq. (|93j) and the physics is that described by the square-lattice quantum liquid of c and si fermions 
considered in this paper. However, that for x £ (0,x*) the magnitude of the energy parameter 2|A| decreases linearly 
in a; is a result found below for the smaller range U/At £ («o,«7r) corresponding to x* values obeying the inequality 
x* < l/n. Therefore, our expressions of physical quantities involving such a linear- a; behavior of 2|A| are valid only 
for that U/At range. Moreover, often we consider the smaller range U/At £ (uo,Mi) for which r c w 2r s sa 7rx». 
Fortunately, the studies of Ref^S reveal that alike for the parent compounds discussed in Subsection VI-B, the values 
of the effective parameters U and t that lead to good agreement between the square-lattice quantum liquid perturbed 
by small 3D anisotropy effects considered in Refj2& and experiments on hole-doped cuprates with superconducting 
zero-temperature hole concentrations x c ~ 0.05 and x* as 0.27 is reached at the value U/At a u* = 1.525 belonging 
to the range U/At £ (uq,Ui). This is consistent with the critical hole concentration x» being found below to read 
x* « 0.27 for U/At » 1.525. 

Furthermore, the symmetry arguments implying that the half-filling si fermion energy gap |A s i(g)| has the general 
form provided in Eq. (|55)l . indicate here that for hole concentrations < x < x* the gap |A s i(q)| has a similar form 
but with [/i°/2] replaced by |A|. The use of arguments similar to those considered for the x = and m = problem 
suggests that for < x < x*, m — 0, and finite values of U/At > uq the gap function |A s i(g)| appearing in the si 
fermion energy dispersion e s \{q) of Eq. (j<?3"|) and related quantities are given by, 

\^M\ = 1^(^)1 = |A|Fll(?o); F' l 1 (gb)= |e ' l(gbgl) " e ' l(gtea)l ; < x < x, , m = 0. (94) 

Note that |A sl (<fJ£)| = (A^Of^)] = and \e sl (q£») - e sl (0)| = W sl and thus, 

F^(q )£ (0,1); F^(q Bsl ) = | cos20| ; F ] } 1 (qg 1 )=0; ) = 1 , (95) 

where the expression F^qBsi) = | cos2c/>| is found below for < x <C 1 and is expected to be a good approximation for 
< x < x* provided that U/At £ (uq, u„). Hence the maximum energy gap magnitude of Eq. (f94|) of the si fermion 
energy dispersion equals the maximum si fermion pairing energy per spinon |A| and corresponds to momentum 
values belonging to the si boundary line whose auxiliary momenta of Eq. (|55|) point in the anti- nodal directions. 
Such behaviors are consistent with the d-wave symmetry of the x = and m = problem surviving at least in part 
for finite values of the hole concentration x in spite of the x — long-range antiferromagnetic order being replaced 
by a short-range spin order. 



2. Shape of the si band for x > and x dependence of the energy parameter 2| A| 



For x > and m = ground states the si band is full and enclosed by the si boundary line whose momenta are 
defined by the relation given in Eq. (fr?5|) . Only for the x = and m — Mott- Hubbard insulator does such a line 
contain the nodal momentum — [7r/2,7r/2] and anti-nodal momenta — [tt, 0] and — [0,7r]. For x > and m = these 
momenta are replaced by the nodal momentum and anti-nodal momentum related to the auxiliary momentum values 
qg al and of Eq. (J7DJ) with components q^ sl = -[q^ sl /V2,q^ sl /V^] and q^ x = — fos^O], respectively, as given 
in Eq. For x > and m = one has according to Eq. (I7T1) . q^ sl < tt/^/2 and q^^ < tt, consistently with the 

momentum area enclosed by the si boundary line being given by (2tt/L) 2 N% s1 = (2ir/L) 2 [(1 — x)/2]N% — 2tt 2 (1 — x) 
and thus decreasing upon increasing x. This deformation and contraction of the si band boundary line upon increasing 
x is found in Ref.— to play an important role in the square-lattice quantum liquid physics. For instance, strong evidence 
is given in that reference that it is behind the positions in momentum space of the incommensurate low-energy peaks 
observed by neutron scattering in LSCO^. 

Fortunately, owing to symmetry arguments the x dependence of the absolute values g^ sl and qg^[ of the nodal 
and anti-nodal si boundary momenta, respectively, can be estimated for approximately x < x c \ and U/At > uq. It is 
useful to consider the si band auxiliary nodal arc momentum qff arc such that, 



Hal. arc 



tt/2 
tt/2 



?b.i (96) 
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and qg sl = -[7r/2,7r/2] - q^ arc - 

The starting points of the study of the x dependence of q^ sl and q^i are that the area enclosed by the si 
boundary line is exactly given by 27r 2 (l — x), the si band momentum values of states belonging to the same V 
tower are independent of U/At, the < x <C 1 and m = ground state belongs to the same V tower except 
for very small U/At values alike the x = and m = ground state does for the whole range U/At > 0, and the 
momenta belonging to the ground-state si boundary line are determined by the si fermion auxilary energy dispersion 
e° sl (q) = e^dA^j^q) = e J(q ) through the relation given in Eq. (1551) . According to Eq. (|57)) . such an energy 
dispersion can be expressed in terms of the elementary function e s i(g) = e s i(—q) whose ID momentum q belongs 
to the domain q £ (— g^fi' Qbsi)- At a; = and m — such an elementary function is found above in Subsection 
IV- to be given by e s i(q) ~ — [Wji/2] cosg. Concerning physical quantities involving the elementary- function sum 
e s i(9oxi) + e si(Qox 2 ) such as the dispersion e^g) of Eq. (p7|) and the related expressions found in the following for 
IbsI ano - Qb^ii ^ i s a good approximation to consider that e s i(g) ~ — [W r s ° 1 /2] cosg for < x <C 1, m = 0, and the 
whole q range q £ (-q^i, q£%). 

For approximately U/At > uq and x < x c \ the shape of the m = ground-state si boundary line and thus the 
magnitudes of both g^ sl and q^i are nearly independent of U/At. Indeed, that such a shape and the eigenvalues 
of the momentum operator are the same for states belonging to the same V tower implies that there is an interplay 
between the U/At dependence of the ground-state c Fermi line and si boundary line. Specifically, if for some region 
of parameter space one of such lines does not depend on U /At implies in general a similar behavior for the other. In 
turn, that the shapes of such m = ground-state lines depends on U /At means that upon varying U /At the ground 
state belongs to different V towers. 

Consistently, if for some range of hole concentrations x and U/At values the c Fermi line is U/At independent 
provides strong evidence that the ground state refers to the same V tower for such a parameter-space region. If so the 
si boundary line is U/At independent as well. It turns out that for U/At £ (uj(x),oo) and small x where Ui(x) ~ 
for < x <C 1 and Ui(x) ~ uq for hole concentrations x x c \ and smaller than x c \ the c Fermi hole momentum reads 
q Fc ~ ^/xtt2 so that the ground state belongs to the same V tower or neighboring V towers and the c Fermi line 
shape does not depend on U/At. Then that strongly suggests that the same occurs for the si boundary line. 

We recall that the x = and m — ground state refers exactly to a single V tower for the whole range U/At > 0. 
Consistently, the physics of such a behavior is then for very small hole concentrations < x <C 1 that remains 
true for U/At £ (iii(x), oo) where Ui(x) « and upon further increasing x that remains true as well for some 
range U/At £ (uj(x),oo) where we estimate that u.^x) £ (0,Uq) for approximately x £ (0, x c i). For the range 
U/At £ (ui(x), oo) for which the x > and m = ground state refers approximately to a single V tower the ground- 
state c Fermi line and si boundary line shapes are nearly independent of U/At. It is expected that for x < x c \ 
and U/At > uq the ground state belongs to neighboring V towers in Hilbert space so that both the c Fermi and 
si boundary line are independent of U/At or nearly independent of U/At. Concerning some properties, it is a good 
approximation to consider that q Fc ss y/xlr2 remains nearly independent of U/At for uo < U/At < u\ and x £ (0, £*). 
Note however that for x £ (x c i,x*) the ground-state hole- like c Fermi line may not be a perfect circle anymore, its 
shape and that of the ground-state si boundary line having some U/At dependence. 

The energy parameter | A| of Eq. (|94p is the maximum energy bandwidth of the si boundary line of the si fermion 
dispersion e s i(g) of Eq. (|93[) and thus refers directly to that line. We recall that for the si fermion auxiliary dispersion 
€gi(q) of Eq. (|67p such a line has vanishing energy bandwidth, yet this is not so for the si fermion energy dispersion 
e s i(g), owing to its d-wave-like structure. For U/At — > the energy parameter Ao = lim x ^o |A| equals one half the 
chemical potential no/2 — lim^-j.o A*/2. For approximately U/At > Uo and < i « 1 the si fermion-pairing-energy 
per spinon shift ||A| — Ao| plays concerning the si boundary line a role similar to that of the chemical-potential shift 
|/i — fio\ = [— e c (qpd Xi ) — e c{qFcx 2 ) ~ t 1 ^] for the Fermi line. Indeed, consistently with ||A| — A | being the shift of 
the maximum energy bandwidth of the si boundary line upon slightly increasing x, for a small hole concentration 
< x <C 1 such a shift controls the si boundary line energy of the si fermion auxiliary dispersion £®i(q) and si 
boundary line nodal energy of the corresponding dispersion e s x(q), respectively, as follows, 

[e s i(gs s ixi) +e s i(q Bs ix 2 )] = -Sfi a i ; #Msi = j - A | , < x < 1 , U/At>u , (97) 

where the factor 1/4 is justified below and is related to the U/At — > oo ratio W^/Aq = 1/4. The si boundary energy 
5fi s i = — [e s i(qBsixi) + e s i(qBsix 2 )] of Eq. (|9"7| vanishes for x = and is small for < x -C 1. For approximately 
U / At > uo and < x -C 1 one can estimate its magnitude. The chemical potential shift |jtx — A*o| obeys the inequality 
\fji — jtio| = [—e c (qpc X1 ) — e c{<lFcx 2 ) ~~ Mo] < xtt[W c /2] where W c = \Wf + W//}] — 8t is the c fermion energy dispersion 
bandwidth, the chemical potential \x is given in Eq. (|72)) . and the equality [— e c (q F /; Xi ) — e c (g^f X2 ) — Mo] = xir[W c /2] 
is reached for U/At — > oo. Similarly, symmetry implies that the si boundary-line energy 5fi s \ = [— e s i(qBsixi) ~ 
£si(qBsix 2 )} obeys a corresponding inequality <5/i s i < X7r[W/ 1 /2] and is given by 6/j, s i = x ir [Wg ± /2] for U/At — > oo. 
Furthermore, that the si boundary line is nearly U/At independent for approximately U/At > uq and x < x c i requires 



41 



that 5{j, s i ~ xir[W%i/2] for such a parameter-space region so that the only U/At dependence of 8[i s i occurs through 
that of the energy bandwidth W®± . Hence, 

W\ 

[esiiqBsixJ + e s i(qBsix 2 )] = —SfJ, a i ; Sfi s i ~ xn— ^- , x < x c i , £7/4^ > it , (98) 

so that according to the relation of Eq. (|97|). ||A| — A |/4 w x irW^/2 for < x <C 1 and approximately t//4i > uq. 

The fulfillment of Eq. (|9"5|) implies both that 2e s i(— g^ sl ) = — 5/i s i for qozi = 9(te2 = — 9bs1 an< ^ [ e si(0) + 
e si( — Sssi)] = — ^Msi for the auxiliary si boundary-line momentum pointing in the nodal and anti-nodal directions, 
respectively. Such conditions are met provided that, 

QBsi-V2^(l-x); < orc «V2|x; « [tt - v^F] , x < z c i , (7/4t > u , (99) 

where gjj arc is the absolute value of the si band auxiliary nodal arc momentum q^[ arc of Eq. (|96p . The Cartesian 
components of auxiliary momenta belonging to the si boundary line and pointing in or near the nodal directions then 
obey the following equations, 

qssix-L ±£Bsix 2 = [f-Sg^] or q Bs i Xl ±Qbs1x 2 = -W — Sg^]; Sq^[ = V2q^ larc , x < x cl , U/At>u Q . (100) 

Importantly, the U/At independent expressions given in Eqs. (1991) and (|100l) for momenta belonging to the si boundary 
line confirm that the above expression 8fi s i = xttIW^/2] is consistent with the requirement of the si boundary line 
being U/At independent for approximately U/At > uq and x < x c i- Since 5/j s i = ||A| — Ao|/4 for < x <C 1 where 
the factor 1/4 is confirmed below, such a requirement is fulfilled provided that ||A| — Ao| ~ x2ttW^ 1 for small hole 
concentrations and approximately U /At > uq. 

If |A| = there would be no short-range spin order and the si boundary line would be a vanishing energy line for 
all its momenta. However, that only holds for auxiliary momenta pointing in the nodal directions, the energy scale 
|A| being the maximum energy bandwidth of the si boundary line reached for auxiliary momenta pointing in the 
anti-nodal directions. The si fermion auxiliary energy dispersion e®i(q) of Eq. (|67[) defines the momentum shape of 
the si boundary line through the second relation provided in Eq. (|68j) . In turn, the si fermion pairing energy per 
spinon |A s i(g)| given in Eq. (|94[) defines the energy bandwidth of that line within the dispersion e s i(q) of Eq. (|93l) . 
Furthermore, that for small x the shift ||A| — Ao| of its maximum magnitude |A| reads ||A| — Ao| = A5fi s i ~ x 27rW~ 1 
assures that the si boundary line is U /At independent for approximately U /At > uq, as Eqs. (|99)) and <|100[) confirm. 

For < x <C 1 such a maximum magnitude | A| of the si fermion pairing energy per spinon given in Eq. (|94p can 
be expressed as, 

|A| = |A^(g^)| « A ^ = ^-|e sl (0) - e sl (-q£»)\ ; < x « 1. (101) 

si si 

Here the energy parameter W s \ — 2 [e s i(— q^ sl /V2)~ e s i(0)] of Eq. (|73l) . which according to Eq. (|75|) can be expressed 
both as W s i = 2 [e a i(— qs^i) ~ e si{— Sssi/V^)] and W s \ = [e s i(— q^x) — e s i(0)], refers to the energy bandwidths 
W sl = [e°i(<fj sl ) - e° sl (0)] and W sl = [e sX (qg s f) - e sl (Q)] of the si fermion auxiliary dispersion e° sl (q) of Eq. (J67]) 
and between the nodal energy level and zero-momentum energy level of the si fermion dispersion e s i(q) of Eq. (|93l) . 
respectively. 

Within the square- lattice quantum liquid of c and si fermions the U/At dependence of the charge mass ratio 
r c = to£° /m* of Eq. (|83|) controls the effects of electronic correlations on physical quantities related to the charge 
degrees of freedom, including the U/At dependence of the c fermion energy dispersion of Eq. (|81[) . In turn, the 
spin ratio r s = Ao/4W s ° 1 controls the effects of such correlations on quantities associated with the spin degrees of 
freedom. According to the results of Refi^, upon increasing U/At the energy parameter Ao interpolates between 

A « fjP/2 w 16te- n V^/u f or jj/At -> and A « AW° t « n(At) 2 /U for U/At > 1. Hence both fi°/2 and A vanish 
for U/At —> 0, whereas for U/At 3> 1 the energy parameter Ao becomes the energy scale 4Wj\ ~ Tr(At) 2 /U associated 
with antiferromagnetic correlations. Within the square- lattice quantum liquid scheme, the spin ratio r s = Ao/4W s ° 1 
is often parametrize as r s — e~ Aa where A s = | ln(Ao/4W s ° 1 )|. For U/At > it controls the interpolation behavior of 
the energy scale A , which for the whole range of U/At values reads, 

Ao - r s AW° sl = AW° sl ; X s = | in(A /4W°)| , (102) 

where A s has the limiting behaviors^, 

A s = wy/At/U, [7/4* <1; X s ^Atu a /U, u q < U/At < u x ; \ s = , U/At -> oo , 
uoo ~ (uo/tt) 2 ~ 0.171 ; w ~ 1-302; u x w 1.600. (103) 
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The Mo magnitude provided here is obtained taking into account that according to the investigations of Section VI, 
r s = 0.135 7r s» 0.424 for U/At s» 1.525. Indeed within the approximate expression r s m e ~ 4tu °/ u it is then uniquely 
determined by the equation 0.135 it = er u °' 1 - 525 and reads uq ~ 1.302, as given in Eq. (|102[) . That for approximately 
uoo < U/At < u\ the U/At exponential dependence of A s involves t/U rather than y/t/U, as it does for U/At <C 1 when 

r s = e~ 7I V 4t / c/ j is consistent with for large U/At values r s being of the form r s ss [1 — C t/U}. For large U/At only the 
zeroth order term of the expansion e ~ 4tu o/ u sa [1 — Atuo/U] is exact, 4mo ~ 5.2 being an approximate value for the 
unknown coefficient C > of the first order term —Ct/U. Nevertheless, the expression r s sa e —Uuo/u j s expected to 
be a good quantitative approximation for the range uoo < U/At < u\ of Eq. (|103j) . 

It follows from the above analysis that for x < x c \ and the approximate range U /At > Uo the momentum shape of 
the ground-state si band boundary line is independent of U/At. The factor 1/4 of the relation 5/i s i — ||A| — Aq|/4 
provided in in Eq. (|97[) is confirmed from properties specific to the U/At — > oo limit. Indeed, symmetries associated 
with that limit imply that the x = elementary-function expression e s i(q) ~ — [W^/2] cosq is for < x <C 1 and 
U/At — > oo a good approximation for both [e s i(qo Xl ) + e s i(qo X2 )] and [e s i(goa;i) — e si{qox 2 )]- This is in contrast to 
finite U/At values and small finite hole concentrations < x <C 1, for which it is a good approximation concerning 
the quantity [e s i(qo Xl ) + e s i(qo X2 )] but not [e s i(qo Xl ) — egi^oa^)]- It then follows that for < x -C 1 and U/At — > oo 
the energy bandwidth ratio W sl /W^ = \e sl (0) - e a x{-q£&)\/W* x appearing in the |A| expression of Eq. ()101|) reads 
lim^/ 4t ^ TO W^i/W^ « [1 — cos(7r — ?s^i)]/2 = (1 — xir/2). This implies that the si fermion auxiliary dispersion- 
energy bandwidth deviation [W al — W s i] has for < x <C 1 the limiting behavior lim;//4t->ao[I i ^si — W a x] = ^A*si- 
Furthermore, since r s = 1 for U/At — > oo one finds from the use of Eqs. (|101l) and (|102l) that Hmi7/4t_>.oo[Aci — |A|] = 
lini!7/4i->oo l|A| — Ao| = 45/x s i so that the relation S/i s i — ||A| — Ao|/4 of Eq. (1^71) holds in that limit. As confirmed 
above, that the si boundary line is nearly U /At independent for approximately U /At > uo and x < x c \ implies that 
for < x <C 1 and approximately U/At > uo the ratio 5/x s i/||A| — Ao| = 1/4 is independent of U/At and x so that 
the relation <5/i s i/||A| — Ao|/4 remains valid for such a range, as given in Eq. That relation refers to m = 

and very small but finite values of x and assures that the si boundary-line momenta expressions provided in Eqs. 
Eqs. and (|100|) are independent of U/At. Indeed, since Sfi s i oc ||A| — A |, if the ratio 5/x s i/||A| — A | was not 
independent of x and U/At such a symmetry requirement would not be fulfilled. Hence its value 1 /A found here for 
U/At — > oo is valid as well for approximately U/At > uq. 

We recall that for m = the si fermion pairing energy per spinon |A| of the model on the square lattice has a 
singular behavior at x = 0, with the above energy parameter Ao = lim^^o |A| being for U/At > different from 
/xo/2 = jAU^^o > Ao. This is due to a sharp quantum phase transition between the x — and m = ground state 
with long-range antiferromagnetic order and the m = and < i « 1 ground state with a incommensurate-spiral 
short-range spin order with strong antiferromagnetic correlations^. At m = and both for small temperature T > 
and x = and for T > and < x <C 1 the system is driven into a renormalized classical regime where the T = 
and x = long-range antiferromagnetic order is replaced by a quasi-long-range spin order as that studied in Ref. 54 
for simpler spin systems. 

Combination of expressions ([57]) . (|10ip . and (|102[) reveals that for small hole concentrations < x <C 1 and the 
approximate range U/At > uo the spin ratio r s = e~ As where A s is given in Eq. (|103[) controls the dependence on x 
of the energy scale |A| and energy bandwidth W a i, which read, 

|A| « 4r s |e sl (0)-e sl (-^f 1 )|«fl-4)2Ao; x° « ^ = ^± , 

W.i = |e sl (0) - e sl (-^)l » (l - > < x < 1 , C7/4t > w , (104) 

respectively. Since the spin ratio r s is a monotonous increasing function of U/At, for the range U/At > u the 
parameter x° is an increasing function of U / At as well. We emphasize that the U / At dependence of the energy scale 
2 1 A | and parameter x® does not contradict the U /At independence of the si boundary line shape in momentum space. 
Such an independence is related to the invariance under the electron - rotated-electron unitary transformation of the 
momentum operator. Indeed, except for small values of U/At the m — ground states refer for small a; to a single V 
tower. In turn, the Hamiltonian is not invariant under that transformation. Therefore, the si boundary line energy 
bandwidth |A| depends on U/At. In turn, that as given in Eq. (|97l ) [Ao — |A|]/<5/i s i = 4 is U/At independent for 
0<i<1 and the range U /At > uq following indirectly from the U/At independence of the si boundary line shape. 

In Appendix C it is shown that except for U/At — > oo the elementary- function expression e s i(g) w — [W^/2] cosg, 
which is a good approximation for 0<i<1 concerning physical quantities involving [e s i(qo Xl ) + e s i(qox 2 )], does not 
lead to the correct expression of those involving the quantity [e s i (qBsi xt ) — e s i (qBsi x 2 )] such as the energy parameter 
2|A| and energy bandwidth W s \ of Eq. (|104p . The point is that except for U/At — > oo there are extra terms in e s i(q) 
that do not contribute to [e s i(QBsixi) +e s i(9ssix 2 )]- I n that Appendix a corresponding e s i(q) expression valid up to 
first order in x and given in Eq. (|Cip of that Appendix is evaluated. 
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As discussed in the following, for the range U/At > uo considered in our studies the short-range spin order prevails at 
zero temperature for finite hole concentrations < x < x» where the critical hole concentration x* is for approximately 
uq < U /At < given by the parameter x® of Eq. (|104[) , consistently with the central role plaid in the quantum liquid 
physics by the spin ratio r s — Ao/4M / s ° 1 = e _As . 

3. The critical hole concentration x t 

For the square-lattice quantum liquid the linear dependence on (1 — of the maximum magnitude of the si 

fermion spinon-pairing energy 2|A| given in Eq. (|104|) is obtained here for small hole concentrations 0<i<1. The 
studies of Rcf. 2 - reveal that the fluctuations of a phase denoted in that reference by B\ and the corresponding amplitude 
gi = |(e t9l )| play an important role in the physics of the square-lattice quantum liquid. For zero temperature that 
amplitude is denoted by g\ and for U/At intermediate values in approximately the range U/At € {uq^u^) a consistent 
picture is reached provided that the linear- a; decreasing of the energy scale 2|A| found above for < x <C 1 refers 
to the whole range < x < x» . This implies that for U/At e (uq,^) the critical hole concentration x* equals the 
parameter x° of Eq. (|104l) . Indeed, the investigations of that reference lead to g\ « (1 — x/x») 2l/ for 0<(i,-i)«1 
where x* is a critical hole concentration above which there is no short-range spin order at zero temperature, = 1 the 
dynamical exponent, and owing to symmetry arguments the exponent v is expected to read v = 1 both for < i « 1 
and < (x* — x) <C 1, respectively. This is fulfilled provided that x» ~ xj and v = 1. The fulfillment of symmetry 
requirements confirms that g\ (1 — x/x») is a good approximation for < x < x*. In turn, for U/At > u n and 
thus 1/tt < x" < 2/ir ss 0.64 either 2|A| decreases linearly or faster upon increasing x so that x» < x®. Based on the 
available information one then finds, 

2r 2e~ Xs 

x* « x° = — = , for u < U/At <u^- x* < x° , for U/At > u K , (105) 

where whether x* x° or x* < x* for U/At > u v depends on 2|A| decreasing linearly in x or faster, upon increasing 
x. It follows that the critical hole concentration x* has for U/At> uq the limiting values, 

x* « 0.23, U/At = u « 1.302, 
w 0.27, U/At = u* = 1.525, 
w 0.28, U/At = m ~ 1.600 

= - w 0.32 , U/At = u^> ui 

< - w 0.64, £/ /4i — ^ oo . (106) 

The fluctuations of the x dependent phase 6\ refer to short-range spin fluctuations. According to the investigations 
of Ref.— , such fluctuations become large for x — > x*. For a temperature range 0<T<T*,0<x<x*, and 
approximately uq < U/At < where T* w giAo/ks is the pseudogap temperature the fluctuations of the phase #i 
remain small so that the amplitude g\ = |(e l8l )| remains finite. This confirms that the phase Q\ and corresponding 
amplitude g\ are directly related to the short-range spin correlations, which prevail provided that g\ > 0. For 
approximately U /At > uq that occurs for < x < x* and T < T*. 

Since for U/At € (ito, %) the prefactor rj (1 — x/x*) of 2|A| = (1 — x/x*) 2Ao obeys for the hole concentration 
range < x < x* the amplitude properties gi = |(e iei )| < 1 and g\ — \ — > for x — > x», the studies of Refi^ 

identify the zero-temperature magnitude $i 2Ao of 2|A| = g± 2Ao with the maximum magnitude of the si fermion 
spinon pairing energy for the whole range of hole concentrations < x < x* of the short-range spin ordered phase. 
According to such an identification, for uq < U/At < the corresponding linear behavior in (1 — x/x*) of that energy 
scale remains dominant for m = and a hole-concentration range < x < x* where 2|A| = (1 — x/x*)2Ao > 0. 

The expression of Eq. (|101[) is only valid for < x <C 1 and for larger values of x a linear behavior of the energy 
scale 2|A| does not imply a similar behavior for the x dependence of the energy bandwidth W s \ given in Eq. (|104[) 
for < x -c 1. That the dominant processes behind the linear decreasing of the magnitude of the si fermion 
spinon-pairing energy 2|A| given in Eq. (|104j) remain dominant for uq < U/At < u n and larger values of the hole 
concentration x in the range x € (0,x*) plays a central role in the physics of the square-lattice quantum liquid: It 
implies that for approximately uq < U/At < u n such a pairing energy decreases linearly with x until it vanishes, as the 
limit reached. Indeed, extension to a larger hole-concentration range of the linear dependence on (1 — x/x„) 

given in Eq. (|104[) for that energy scale would lead to negative values of 2|A| for x > x*. The results of Ref»^ confirm 
that <7i = for uq < U/At < and x > x* so that the zero-temperature expression 2|A| = <7i2Aq is only valid for 
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< x < x*. Consistently with the energy scale 2|A| being associated with the short-range spin order, it follows that 
the hole concentration x* of Eq. (|105[) corresponds to the critical hole concentration above which the short-range spin 
order disappears at zero temperature for the intermediate U/At range U/At e (mo,?/^). Consistently with both gi = 
and 2|A| = for x > x*. Often our analysis focuses on the smaller range U/At S (uo,ui) for which the charge and 
spin ratios obey approximately the relation r c 2r s nx* and r s ss e~ Uu °/ U , Indeed for finite U/At > values 
one has that r s < r c < 2r s and such a relation is not fulfilled. In addition, for such U/At values the approximate 
expression r s ss e ~ 4tu o/u [ s no ^ valid anymore and thus we could not access the accurate value of u-^ > u\ at which 
the x* magnitude reads x* = 

For m = the si fermion spinon-pairing energy 2Ao = lim^^o 2| A| is the energy below which the short-range 
spin order with strong antiferromagnetic correlations survives for hole concentrations < x -C 1. The short-range 
spin order that emerges for small hole concentrations and temperatures T > is similar to that occurring for x = 
and T > and was studied previously as for instance in Refi^. The results of Refi2£ indicate that if one adds 
to the square-lattice quantum liquid considered in this paper and in Ref3 a small perturbation associated with 3D 
anisotropy effects brought about by weak plane coupling, for hole concentrations x c < x < x* and temperatures 
T < T c w gAo/2kB the short-range spin order coexists in such a quasi-2D system with a long-range superconducting 
order. Here g = [{x — x c )/{x* — x c )]g\ for uq < U/At < u\ and the critical hole concentration x c reads x c ~ 1CP 2 
for the small 3D anisotropy effects considered in that reference. The maximum energy gap magnitude of Eq. (|94|) of 
the si fermion energy dispersion equals the maximum si fermion pairing energy per spinon |A| and corresponds to 
momentum values belonging to the si boundary line whose auxiliary momenta of Eq. (1551) point in the anti-nodal 
directions. For approximately uq < U/At < it reads, 

|A| = |A sl ( g -^H| - |A^(^)| « A ( 1 - — ) ; < x < x* , m = , u < U/At < . (107) 



V. c AND si FERMION VELOCITIES AND THE ONE- AND TWO-ELECTRON EXCITATIONS 

In this section we derive the c and si fermion velocities from the corresponding energy dispersions obtained above 
and discuss their role in the one- and two-electron spectrum. 

A. The c and si fermion velocities 

The c and si fermion group velocities derived from the energy dispersions are the elementary building blocks of the 
velocities associated with one- and two-electron excitations. They read, 

£i(g) = V>e al ($); 7 s ?(0=V f e° 1 (9); V s i (o) = -V*| A sl (q}\ ; V c (q h ) = v> e c (q h ) . (108) 

Such velocities can be expressed as, 

V sl (q) = V sl (q) e Mit) ; V s °i(<l) = V?i(S) ; V.f% = V£(fl ; V c (q h ) = V c (q h ) e^ h) , (109) 

where V s \{q), V^{(q), V s ^(q), V^(q), and V c (q h ) are the velocities absolute values and the unit vectors are obviously 
given by, 

_VsM. - _. __ V s ?(q) , ? _V c (q h ) 

e ^®-v^> e ^-vm' ^'v^' ^ h) ~wr (110) 

Within the notation used here a unit vector is such that its components are given in terms of the angle <f> that 
defines its direction as in Eq. (TT9"1) . 

The velocities V s \{q) and V s f{q) of Eq. (JTUHJ) can be expressed as, 

vJl(q) = ({A^-yvJi'Hqo); v„t% = ({A^-ytfHqo), (HI) 

where is the 2x2 matrix appearing in Eq. (|4U1) . [AjJ -1 its inverse matrix, ([^4fi] _1 ) T the transposition of the 
latter matrix, and the velocities V^'"(g ) and T^'"(g*o) are given by, 

V s iH?o) = V fo e°J(q ) ; V s ^(q Q ) = ~% I A lM)\ ■ (112) 
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Here e°J(q ) = ([A^]- 1 q) is the auxiliary energy dispersion of Eq. §7$ and (A^^o)! = lA^^]" 1 ^! the 
auxiliary gap function of Eq. (j9"2")l . 

For si band momenta q at or near the si boundary line and hole concentrations in the range x c \ < x < x c2 one 
found above in Section III that A d sl = A F , [A^]' 1 = A F d , and ([AfJ^ 1 ) 7 = A d F where A d F is the 2 x 2 matrix given 
in Eq. (|44j) so that expressions (|111[) simplify to, 



v?M = A d F KiHqo) ; v s ?(o) = A d F v s V l (qo) . 

The si velocities of Eq. (|112[) and the c fermion velocity have Cartesian components given by, 



(113) 



V£ '"(<fo) = -sgn{e s i(9 0a;i ) - e s i(q 0<B3 )}^p- [v sl (q 0xi ),-v s i(q 0x2 )} , 

VVsl 

V s °i(qo) = [v sl (q 0xi ),v sl (q 0x2 )}; V c (q h ) = [v c (q^ ), v c (q% 2 )] , 



IAI 



(114) 



where e 7 (g) and Vj(q) with 7 = c, si are the elementary functions and elementary velocities, respectively, of Eq. (|69l) . 
|A| one half the energy parameter defined in Eq. (|104[) . and W s i the energy bandwidth of the si fermion auxiliary 
dispersion provided in Eq. (|571) . 

The following velocities vanish: V c {tz) = V c (0) = V^°(0) = 0. The velocities Vj±(q) and V^(q6) have the same 
absolute value Vji = V s ^(q), which can be written as, 



V£ (q) = |^ G sl (q) ; G sl (q) = (<fo) = ^ 



Vsiiqox!) 2 + v sl (q 0x2 



where the function G'|J 1 (go) is such that, 



Gl(q ) 6(0,1); G^(g Bsl ) = | sin2«>| ; G^(gf sl ) = l; G» 1 (^[)=0. 
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(115) 



(116) 



Here the expression G^^qBsi) — \ sin2</>| is exact for < x <C 1 and is expected to be a good approximation for 
< x < x* and the auxiliary momenta q Bs i and q B ^[ are particular cases of the general auxiliary momentum q Bs \ 
defined in Eq. (|55j) and point in the nodal and anti-nodal directions, respectively. The corresponding momenta 
q B ^i and q B ^ N are particular cases of the general si boundary-line momentum q Bsl of Eq. (j45|) . If as discussed in 
Section III one assumes that the hole concentration Xh above which the Fermi line is particle like obeys the inequality 
x C 2 < Xh < x* (rather than Xh > x*), for the hole-concentration range Xh < x < 2* for which the Fermi line is particle 
like that the range of that the expression G^^qBsi) = | sin2</>| refers to is given in Eq. (|37|) so that the minimum 
magnitude of Gjj 1 (qs s i) rather than vanishing reads Gl^qssi) ~ %4>an where <Pan {%) is small. Such a minimum 
magnitude is reached at (j) — c/)an and <p = tt/2 — 4>an instead of at <f> = and (f> — 7r/2, respectively. 
At the si boundary line the velocity V s i(q) reads, 
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where the second expression is valid for x c \ < x < x c2 . To reach such an expression we used the relations (|113p 
and that A d ± = A d F for momenta at the si boundary line and x c i < x < x C 2- The equality V s i(q Bsl ) — V s ^(q Bsl ) 
implies that the velocity V s ^(q Bsl ) plays an important role in the quantum liquid physics. The corresponding auxiliary 



velocity V^'" (<fssi) and the absolute value V s ^(q Bsl ) of both V^(q Bsl ) and V^'" (q Bs i) are given by, 
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V^{q Bs i) = sgn{e sl (q Bslxi ) - e sl (q BslX2 )}-^- [v sl (q Bslxi ), 

W s i 



-v sl {q BslX2 )] ; V£(q£ sl ) = U Gl(q Bsl ) , (118) 



V f e al (q) simplifies to V sl (q) = V^(q) 



respectively. 

For the si band nodal directions such that qo Xl — ±qo X2 the velocity V s i(q) 
so that it is not well defined for the four momenta of the si boundary line pointing in the nodal directions. Hence in 
that case we consider instead the velocity at the si boundary-line momenta q B ^ whose auxiliary momenta q Bsl given 
in Eq. (|55|) point in the infinitesimal vicinity of the nodal directions and hence for the quadrant where qo Xl = qa X2 < 
have Cartesian components q Bsl = [— q Bslxi (1 + Sq), —q Bslx (1 + Sq)) where Sq — > 0. (According to Eq. (|99|) one has 
for small values of x and U/At > u Q that q Bslxi — q Bs i Xo = — [ 7r /2](l — x).) For such nodal momentum values the 
corresponding auxiliary velocity reads, 



Kt'teiH-sgn{<M 



w sl 



v s i(q Bs i X1 ) 
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Finally, the use of the energy-dispersion expression given in Eq. (IBTj) in the c fcrmion velocity expression provided 
in Eq. (fTT4|) leads to, 



V c (q h ) = -^, (120) 
m* 



for x < 1 and C//4t > u . 



B. The angle <f>, Fermi-line one-electron spectrum, and corresponding Fermi velocity 

Complementarity to the results on the Fermi line momentum values reported in Subsection III-C, here we study 
several one-electron physical quantities. Unfortunately, except for the expression e c (q) = — U/2 + 2tcosq valid for 
U/At ^> 1 and that of e s \(q) given in Eq. (|Cip of Appendix C for small values of the hole concentration and U/At > uq, 
one does not know accurately the form of the elementary functions e c (q) and e s i(q) that control the dependence on 
the momentum Cartesian components of the c and si fermion energy dispersions given in Eqs. (|67p and (|76p and 
corresponding velocities of Eq. (1108p . However, one is aware of some properties and symmetries of such functions 
as for instance those behind the limiting values reported in Eqs. ([55]) and (111 6[) for the related functions F s i(qB S i) 
and G s i(qBsi) of Eqs. (|M| and (|115l) . respectively. Those imply that the following expression is exact for small hole 
concentration values, m = 0, and U/At > uq and is expected to be a good approximation for finite hole concentration 
values < x < x* provided that uq < U/At < u„ 

2cj) = sgn I ] arctan ( G U<lBsi) \ = ( kkA ( V2Kfffii) \ (121) 

1 "' ' \Fl{q Bsl )) \k h Fx J \ |A s i(g| sl )| J 




Here 4> is the angle of Eq. f|37[) . which defines the direction of the hole Fermi momentum whose expression is given in 
that equation and in Eq. ((38]). 

The dependence on </> of the angles ^ = 4>fi(4>) and § d c = (j) d ((j)) provided in Eqs. and ([ST))) , respectively, is 
valid only for the isotropic Fermi- velocity range x c \ < x < x C 2- Such angles define the relative directions of the c 
fermion hole momentum and si fermion momentum on the right-rand side of Eq. Q38p for the hole Fermi momentum, 
which for that range are perpendicular to each other. In turn, expression (|121j) for 2<j) does not involve any relation 
between c fermion and si fermion quantities. Indeed, it involves the ratio V^(q^ sl )/\A s i(q^ sl )\ that refers only to si 
fermion quantities. This is why Eq. (| 12 1[) is a good approximation for the whole range < x < x*. If as discussed in 
Section III the hole concentration Xh above which the Fermi line is particle like belongs to the range x c i < Xh < x» , 
then for Xh < x < x* the angle </> of Eq. (|12ip belongs to the ranges given in Eq. (|37]). (The angle 4*an appearing in 
such ranges vanishes for x < xt and is small for x S (xh,x„).) 

Consistently with Eq. (|12ip . the absolute value V^ sl = V s ^(qg sl ) of both the velocities Vji (<fj sl ) and (qBsi) 
of Eq. (TlT7l) reads, 

V^ 1 = ^f(^ 1 ) = J^| S ia2^|. (122) 

Since the angle 4>an is small, for the hole-concentration range Xh < x < x* we consider corrections up to first order 
in (pAN, so that max-j^E^} ~ |A|. Moreover, the minimum magnitude of the velocity (|122p reads min{V^ sl } ps 
[<j>AN V2]\A\. In turn, for < x < Xh its minimum magnitude is zero for momenta q^fi 1 whose auxiliary momenta 
point in the anti-nodal directions. 

Hence the above two velocities can be written as, 

Vg sl = V* sl e a,, ; ? s ? - 11 (q Bsl ) = e . A .„ . (123) 

r osl ^Bsl 

Here the angles (j)gfi = ^(^l) and 0Bai = ^fi (&si) read : 

^flil = <t>Bal + <t>F i B ai = ~ arCtan 1 \ I ( 124 ) 

\Vsl\qBslX! ) / 

where the expression for </>g^ is valid for x c i < x < x C 2 and follows from Eq. (|117|) . 

For < x < the velocity (|122l) is a function of the angle cf> with the following limiting behaviors, 

V£ sl = 1 -^, <t> = n/A; V£ sl =0 7 <f> = 0, tt/2 . (125) 
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For Xh < x < the maximum magnitude remains the same whereas Vg sl rs [4>an \/2] | A | for 4> — 4>an and 
4> = [tt/2 — 0ajv]- That for < x < x^ the velocity V^ sl vanishes for auxiliary momenta pointing in the anti- nodal 
directions together with the si band remaining full for x > and m = ground states has important physical 
consequences. Indeed, finite-energy excitations involving creation of si fermion holes whose auxiliary momenta point 
in the anti-nodal directions refer to real-space states. The energy of such excitations equals the pseudogap energy 
scale. As further discussed in Ref.— , such pseudogap real-space states are observed in experiments on hole-doped 
superconductors^. 

In turn, the velocity V c {q h ) of Eq. (jl20[) is for q h — qfi^ given by, 

V c (Qh d ) = -V Fc e> ; V Fc = « ^ , (126) 

m* m* 

where the angle is for x € (x c i, x C 2) provided in Eq. (1501) and the mass m* = l/2r c t is that of the c fermion energy 
dispersion given in Eq. It is associated with the charge mass ratio r c of Eq. 

According to the c Fermi-line definition of Eq. (|68l) and following the form of the general c fermion energy dispersion 
expression provided in Eq. (|76p . the latter energy vanishes for c band momentum values belonging to the c Fermi 
line. The use of the c and si fermion energy dispersions given in Eqs. (|67p and (|93l) and related expressions provided 
in Eqs. (1M1) - (|107[) in the general energy functional of Eq. (|55|) leads for m < U/At < and hole concentrations 
< x < x* to the following general expression for the one-electron energy spectrum at the Fermi line whose hole 
Fermi momenta are given in Eq. (|38l) . 

E F = ft + 5E F (<l>)] SE F (^) = -e sl (q£ sl ) = \ A sl = |A|| cos20| . (127) 

Here the zero-temperature chemical potential [i is given in Eqs. (jT2"j) and (|82p. it and 5Ep(4>) are the isotropic and 
anisotropic Fermi-energy terms, respectively, and the equality |A s i(<fJ sl )| = |A||cos20| follows from Eqs. (|121[) and 
()123|) . The Fermi-line anisotropic and gapped one-electron spectrum of excitation momentum and energy k F and Ep, 
respectively, has a d-wave-symmetry like structure. Its relation to the physics observed in the hole-doped cupratesi - — 
is discussed in Refj^. 

The value of the momentum q^ sl appearing in the energy spectrum (1127[) refers according to Eq. (|38[) to exactly 
one value of the Fermi hole momentum kp . The inverse relation is two- valued, with each Fermi hole momentum k F 
corresponding to the momenta q£ s i and q^ s i associated with two alternative and degenerate one-electron excited 
states with the same excitation momentum k F and energy Ep given in Eqs. (|55[) and (|127p . respectively. Let us 
confirm that such excited states refer to electrons at the same point of the Fermi line but with different Fermi velocities. 
This is equivalent to show that the electronic velocity at the Fermi hole momentum kp depends on the doublicity 
d = ±1, unlike the corresponding energy Ep and momentum kp. 

Since following the relation between quantum overlaps and the transformation laws of the si fermion operators 
discussed in Ref^ and according to the numbers and number deviations of Table IIVI of Appendix A creation of one 
electron involves creation of one c fermion and one si fermion hole, the Fermi velocity is straightforwardly given by, 

% d = \Vc{q^) - t^Li)] = V F d e K , (128) 

where 

V$ = V Fc V /l-2r A cos(^-<0+'i 5 

Here is the velocity ratio of Eq. (|51?|) and <^ sl and (\> d c the angles that define the direction of the si and c fermion 
velocities of Eqs. (|123p and (|126l) . respectively. 

The form of the expressions given in Eq. (| 129[) confirms that both the absolute value and angle of the Fermi velocity 
depend indeed on the doublicity d = ±1. Note however that r\ = Vg sl /Vp c — t)a \ sin2</>| <C 1 for x > x c \ so that the 
velocity Vp 1 of Eq. (|129[) becomes independent of both d and and given approximately by Vp ~ Vf = V Fc . 

C. Charge and spin excitations 

The Fermi velocity of Eqs. (|128|) - (I129[) involves both the charge c fermion velocity V c (q h ) and spin si fermion 
velocity Vj± (q). In contrast, the velocities corresponding to the charge and spin excitations involve only the velocities 



j / sin^J? — r a sin^n''? \ 

4 F = arctan ^ ^ . (129) 

V cos0° - r A cos(f> B ' sl l 
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V c (q h ) and V s i(q), respectively. For two-electron charge and spin excitations the matrices and A d F are the 2x2 
unit matrix so that the si boundary line is not deformed. Therefore, <fj sl — qBsi for the si boundary line of such 
states, so that the si fermion velocity Vf^{qBsi) and c fermion velocity V c (qp c ) are independent and read, 

V,i = Vtx e A ; V c (qh) = ~V F c e 0c , (130) 

respectively. Here the velocity absolute values Vg sl and Vj? c are those provided in Eqs. f|122[) and (|126[) . respectively, 
the angle <^ equals that given in Eq. (|124j) and thus reads, 

4>7i = - arctan ^ , (131) 

and 4> c = <p. 

The c fermions describe the charge degrees of freedom of the rotated electrons that singly occupy sites and thus 
carry the electronic charge — e, which remains invariant under the electron - rotated-electron unitary transformation. 
Hence they couple to external charge potentials and carry the elementary charge currents. Since in contrast to the si 
fermions the c fermions do not emerge from a Jordan- Wigner transformation, their direct c - c fermion interactions 
vanish or are very weak and the c fermion elementary charge current reads, 

3c(q h ) = -eauV c (q h ); J c (g£ c ) = -e % e 0c+7r « -e a^2l g . au = . (i 32 ) 

m c me rric 

Here m p is a renormalized transport mass and q Fc ~ y/xn 2 is a good approximation for hole concentrations x G (0, x c i) 
and is expected to be a reasonably good approximation for x £ (x c i,x C 2) provided that U/At > Uo- Since for finite 
values of J7/4t the model ([1]) does not commute with the charge current operator, it follows that m p < = l/2t 
for U/At > uq where the equality refers to the limit U/At — > oo. Indeed, the Drude peak exhausts the conductivity 
sum-rule both for U/At — > and U/At — > oo and otherwise its spectral weight is smaller than that associated with 
such a sum- rule. For U/At > uq one then estimates, 

^=2tmgwr c ; = U/At > u , (133) 

mg° nic r c r£ 

where r c is the charge mass ratio given in Eq. (1831) . For instance, for U/At w u* = 1.525 one finds m p c /ra* c — \/oljj ps 
= (ttO.27) 2 w 0.72, whereas mg/m* = 1 for [7/4t -> oo. 

For the Hubbard model on the ID lattice the general description of Refi^ also applies and since for that model 
there is only c and si fermion zero-momentum forward-scattering, the transport mass m p c can be evaluated explicitly 
by suitable use of the ID exact solution and is given in Eq. (139) of Ref.— . 

Finally, we consider the general spin spectrum. According to the results of Ref. 9 and numbers and number deviations 
of Table ITVl of Appendix A, both the spin-singlet and spin-triplet excitations relative to the m = ground state involve 
creation of two si fermion holes. The use of Eqs. (|65|) and (|66|) then leads to the following general spin spectrum, 

SE spm = -e sl (q)-e sl (q'); SP = Sq° -q-q'; 6q?=Tr,x = 0; 5$ = , x > , (134) 

where the momentum shift 5cQ refers to the momentum q^ of Eq. (|34[) for 7 = c. For spin excitations of the square- 
lattice quantum liquid one has that 5q® = 7? = [w, tc] at x = and <5<f c ° = for x > for initial m = ground states. 
Derivative of the excitation energy 8E sp i n of Eq. (|134|) relative to the corresponding momentum 5P confirms that 
the spin velocity is |V s i(g) + V s i(q')] and hence involves only the si fermion velocity. 



VI. RELATION TO THE EXACT SOLUTION OF THE ID MODEL, AGREEMENT OF THE 
SQUARE-LATTICE QUANTUM LIQUID WITH RESULTS OBTAINED BY THE STANDARD 
FORMALISM OF MANY-BODY PHYSICS, AND RELATION TO OTHER SCHEMES 

The general quantum liquid of c and si fermions introduced in Refi 9 and further studied in this paper refers to the 
Hubbard model in the one- and two-electron subspace on a ID or square lattice. For both such lattices that quantum 
problem is non-perturbative in terms of electron operators so that rewriting the theory in terms of the standard 
formalism of many-electron physics is an extremely complex problem. 

In this section we discuss the use of the c and si fermion description to construct a dynamical theory for the 
ID Hubbard model that provides correlation-function expressions both at low and finite energy and recovers the 
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well-known low-enery behavior of that ID quantum liquid. Furthermore, we show that the predictions of the square- 
lattice quantum liquid approach concerning the spin spectrum at half filling agree both with experiments on LCO and 
results obtained by the standard formalism of many-body physics. Finally, we discuss the relation of the square-lattice 
quantum liquid of c and si fermions to other schemes. 

A. The c and si fermion description of the ID Hubbard model 

For the ID model there is an exact solution^—. However, it does not provide correlation-function expressions at 
finite energy. The relation to that solution of the general c and si fermion description and its underlying theory has 
been clarified 9 -^ and one can profit from such a relation to further develop and test the theory. That allows going 
beyond the range of the exact solution alone and derive one- and two-electron spectral-weights and corresponding 
spectral- and correlation-function expressions for both low and finite energ y 58 ' 59 . Fortunately and as discussed in 
the following, in the low-energy limit such a general theory recovers the so called Tomonaga-Luttinger-liquid (TTL) 
behavior and anomalous scaling dimension of spin and charge correlations^ -. 

The general rotated-electron description introduced in Ref£ for the Hubbard model on the square lattice has already 
been considered in Refi 5 ^ for the same model on the ID lattice. Within the notation of the latter reference, the r\- 
spinons, spinons, c fermions, si fermions, independent 77-spinons, and independent spinons considered here are called 
in Ref. 9 holons, spinons, c pseudoparticles, si pseudoparticles, Yang holons, and HL spinons, respectively. (HL stands 
for Heilmann and Lieb 5 ^.) Moreover, the 2^-ry-spinon composite r\v fermions and 2^-spinon sv composite fermions 
considered in Ref. 9 are called in Ref^i cv pseudoparticles and sv pseudoparticles, respectively. The relation to the 
rotated-electron occupancy configurations of all such objects that generate the energy eigenstates of the ID model 
are m Refs£ and 5 ! the same. (The studies of the former reference extended the description to the Hubbard model on 
the square lattice.) 

The c and si fermion scheme refers to the limit N a 3> 1 where the thermodynamic Bethe-ansatz equations of 
Takahashi apply^ 5 .. Exploring the relation of such a scheme to the exact solution has allowed the construction of a ID 
c and si fermion dynamical theory, which provides general expressions for the finite-energy correlation and spectral 
dynamical functions 5 - 8 -, beyond those provided by bosonization and conformal-field theory. Such a dynamical theory 
was recently successfully used to study the finite-energy singular features in photoemission of the organic compound 
tetrathiafulvalene-tetracyanoquinodimethane (TTF-TCNQ) metallic phased. The c and si fermion dynamical theory 
refers to the ID Hubbard model. More recently, other methods for the study of finite-energy spectral and dynamical 
functions of ID correlated systems were introduced 6 -^—. Both the finite-energy spectral-weight distributions studied 
by the c and si fermion dynamical theory for the ID Hubbard model and the methods of Refs^ir— for other 
ID correlated problems include power-law singularities near well-defined finite-energy branch lines with exponents 
depending on the interaction strength and the excitation momentum. 

As discussed in Reb£2, the study of the singularities of spectral functions using models such as the ID Hubbard 
model is not limited to low energies or to weak interactions and the idea of extracting exponents of finite-energy 
spectral functions from the Bethe-ansatz exact solution by combining it with the present c and si fermion description 
has appeared in the pseudofermion dynamical theory for the ID Hubbard model of Ref. 58 . In turn, according to 
the discussions of Ref. 64 , our method relies in part on the integrability of the ID Hubbard model, whereas the 
phenomenology developed in that reference does not require any special property of the underlying microscopic 
interaction. The solvability of the ID Hubbard model is in the N a ^> 1 limit the c and si fermion theory refers to 
associated with the occurrence of an infinite set of conservation law a 36 ' 42 . As discussed in ReL 9 -, within our description 
such laws follow from the commutation of the Hamiltonian with each of the infinite av translation generators q av in 
the presence of the fictitious magnetic fields B au (fj) considered in that reference whose eigenvalues are the components 
of the microscopic momenta of the av fermions. Such microscopic momenta are for the ID model good quantum 
numbers for the whole Hilbcrt space. For the model on the square lattice and x > the c and si fermions undergo 
inelastic collision o 28 ' 29 . In contrast, due to the occurrence of the above conservation laws in ID they have only zero- 
momentum forward scattering, associated with two- fermion phase shifts 58 . Importantly, the studies of Reb£°- reveal 
that when both momenta in the argument of such phase shifts lie at or near the c Fermi line and/or si boundary 
line, their square in units of 2tt fully determines the anomalous dimensions of the primary fields of conformal-field 
theory. The low-energy physics of correlated ID problems has some universal properties described by the TLL 6 - 5 -. 
The investigations of Reb£°- confirm that the ID c and si fermion dynamical theory gives a correct description of the 
TLL behavior and anomalous scaling dimensions of spin and charge correlations. Such a consistency provides further 
evidence of the validity of our general method. Indeed, the finite-energy c and si fermion dynamical theory is shown 
in that reference to recover in the limit of low energy the usual low-energy TLL results. It is an improved ID version 
of the c and si fermion scheme used in our studies for the model on the square lattice, which also refers to the general 
rotated-electron description of Refsi^^. 
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The construction of the c and si fermion dynamical theory of Refi^ profits from the relation of the quantum objects 
of our description to the quantum numbers of the Bethe-ansatz solution. The derivation of the general correlation- 
function and dynamical-function expressions by that dynamical theory takes into account implicitly that the in the 
case of the spin-neutral two-spinon si fermion branch such quantum numbers are the eigenvalues of the si translation 
generator q sl in the presence of the fictitious magnetic field B s \(fj) of Eq. ([2~Tj). The consistency of our description 
concerning the isomorphism between the c and av fermion microscopic momentum values and the quantum numbers 
of the exact solution is addressed in Ref^. 



B. Agreement of the square-lattice quantum liquid with experiments and results obtained by the standard 

formalism of many-body physics 

It is desirable that the results of the square-lattice quantum liquid studied in this paper are compared with those 
of the standard formalism of many-body physics. Unfortunately, such a quantum liquid is non-perturbative in terms 
of electron operators so that, in contrast to a 3D isotropic Fermi liqui d 39 ' 40 , rewriting the theory in terms of it is an 
extremely complex problem. 

In spite of the lack of an exact solution for the model on a square lattice and the non-perturbative character of 
the quantum problem in terms of electrons, in this subsection results achieved by the square-lattice quantum liquid 
description are compared with those obtained by methods relying on the standard formalism of many-body physics. 
Unfortunately, there are not many controlled results for the Hubbard model on the square lattice from approximations 
relying on that formalism. Here we consider the interesting problem of the spin-excitation spectrum of the half-filling 
Hubbard model on the square lattice. 

Within the present description and as discussed below, since the c and si fermion momentum values are for the 
model in the one- and two-electron subspace good quantum numbers, that problem refers to an effectively non- 
interacting limit whereas in terms of electrons it is an involved many-body problem. Fortunately, there are reliable 
results on that particular problem obtained by controlled approximations of the standard formalism of many-electron 
physics: within such approximations its solution requires summing up an infinite set of ladder diagrams, to find the 
spin-wave dispersion of the half-filled Hubbard model on the square lattice in a spin-density-wave-broken symmetry 
ground stated. In turn, within our description the spin spectrum is that of Eq. (|134p . which involves the creation of 
two holes in the si band whose energy dispersion e s ±(q) is given in Eq. ([551) . 

Agreement between the two methods is both a further checking of the validity of our description and a confirmation 
that the c and si fermion interactions are indeed residual and their momentum values good quantum numbers for the 
model on the square lattice in the one- and two-electron subspace. (The results of Refs^^ confirm that for ID these 
momentum values are good quantum numbers for the whole Hilbert space.) As a side result, in this subsection we 
find the magnitudes for U/4t s» 1.525 of the energy scales /i° and 2Aq whose limiting values are given in Eqs. (|A12|) 
of Appendix A and below in Table III, respectively, and of Eq. (fT3"|) . That value of U /At corresponds to the 
critical hole concentration w 0.27 and is that appropriate to the description by the theory of the properties of the 
parent compound LCO. In reference^ 9 , strong evidence is provided that such a U /At value is also that appropriate to 
the hole-doped cuprates with superconducting zero-temperature critical hole concentrations x ~ 0.05 and x ~ 0.27. 

The spin-triplet excitations relative to the x = and m — absolute ground state involve creation of two holes 
in the si band along with a shift tt/N^ of all discrete momentum values of the full c band so that the general spin 
spectrum (|134p reads, 

w(fc) = [-e sl (<?)-e sl (?')]; k = [t? - q- q '} , (135) 

where n = ±[w, ±7r] and the si fermion energy dispersion e s i(q) is given in Eq. (1891) . For x = and m = both 
the c and si bands are full for the initial ground state and since the c band remains full for the excited states one 
can ignore the si - si and si - c fermion interactions. Indeed, then the residual fermion interactions studied in 
Ref. 28 have little effect on the occupancy configurations of the two holes created in the si band upon the two-electron 
spin-triplet excitations. This is consistent with the lack of a c Fermi line for the initial ground state and the lack of 
si band holes other than the two holes created upon the spin-triplet excitation so that in spite of the si - si fermion 
long-range interactions associated with the effective vector potential of Eq. (|2Tj) the exclusion principle, phase-space 
restrictions, and momentum and energy conservation drastically limit the number of available momentum occupancy 
configurations of the final excited states. 

The excitation spectrum (|135[) refers both to coherent and incoherent spin spectral weight. In contrast to the ID 
case where as mentioned in the previous subsection a suitable c and si fermion dynamical theory is availabl e 58 ' 59 , 
for the square-lattice quantum liquid there are within the present status of the theory no suitable tools to calculate 
matrix elements between the ground state and one- and two-electron excited states. Hence, one cannot calculate 
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explicitly spin-spin correlation functions. Within the ID c and si fermion dynamical theory, the sharp features of the 
spin two-electron spectral weight distributions result from processes where one of the two si fermion holes is created 
at the si boundary line. 

The coherent spin spectral weight is here associated with a Goldstone-mode-like gapless spin-wave spectrum. It 
consists of sharp <5-peaks having as background the incoherent spectral- weight distribution. From comparison with the 
results of Ref4^, we have confirmed that such spectral weight is generated by processes corresponding to well-defined 
values of the momenta q and q ' of the general spectrum (|135|1 such that one hole is created at a momentum pointing 
in the nodal directions of the si band and the other hole at a momentum belonging the si band boundary line, as 
expected from analogy with the ID spectral-weight distributions. The incoherent part corresponds to the remaining 
values of q and q' of the excitation spectrum (|135[) . The occurrence of the Goldstone-mode-like gapless spin- wave 
spectrum follows from the long-range antiferromagnetic order of the initial x — and m = ground state. In turn, 
for the x > short-range spin ordered phase the spin weight distribution associated with the general spin spectrum 
provided in Eq. (| 1 34[) has no coherent part. In ID it has not coherent part both for x = and x > 0, due to the lack 
of a ground-state long-range antiferromagnetic order. 

As mentioned above, for the original electrons the problem is highly correlated and involves an infinite set of ladder 
diagrams and no simple analytical expression was found for the spin-wave energy spectrum^. In contrast, for the si 
fermion description it is effectively non interacting and described by simple analytical expressions. Let us profit from 
symmetry and limit our analysis to the sector k x £ (0, 7r) and k y £ (0, k x ) of the (k,ui) space. Within the description 
of the quantum problem used here, for 1/ks T — > oo the coherent spin-spectral- weight distribution derived in Ref4^ 
corresponds to a surface of energy and momentum given by, 
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This is a particular case of the general spin spectrum given in Eq. (|135[) . which corresponds to the above-mentioned 
specific processes associated with the following choices for fr, q, and q' , 
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for the sub-sector such that k x £ (0, tt), k y £ (0, k x ) for k x < tt/2, and k y £ (0, tt — k x ) for k x > tt/2 and, 
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for k y < 7r/2, and k x £ (k y ,Tt) for k y > tt/2, respectively. 



Note that as mentioned above, the components of the si band momenta q appearing in Eqs. (|137l) and (|138[) are such 
that q Xl — q X2 = — an d thus belong to the half-filling si boundary line defined by Eq. (|86l) . whereas those of the 
momenta q' in the same equations obey the relation q' = —q' X2 so that point in the nodal directions. 

Next let us consider the high symmetry directions in the Brillouin zone. Indeed, the curves plotted in Fig. 5 of 
Ref4i refer only to such directions and the use of our above general expressions leads for U /At = 1.525 to an excellent 
agreement with such curves. These directions correspond also to those measured by high-resolution inelastic neutron 
scattering in LCO, as plotted in Fig. 3 (A) of RefJ£. We denote such symmetry directions by MO, TO, XM, FX, 
and XO. They connect the momentum-space points M = [tt,tt}, O = [tt/2, tt/2], T — [0,0], and X = [ir,0] of the 
general spin- wave spectrum provided in Eq. (|136l) . The use of that equation reveals that the spin- wave excitation 
spectrum is in such symmetry directions given by, 



^ro(fc) = Ysin(fej), 

k = [tt, -tt] - [tt/2 - k t ,~Tr/2 - hi] - [tt/2, -tt/2] 
— [ki, ki] ; ki = k x = k y £ (0, tt/2) , 



(139) 
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FIG. 1: The theoretical spin spectra (|139p - (ll43[l (solid lines) plotted in the second Brillouin zone for fi° — 565.6 meV and 
W°i = 49.6 meV and the experimental data of Re&« (circles) in meV. Such theoretical magnitudes correspond to t ~ 0.295 eV 
and U w 1.800 eV, so that U/4t ~ 1.525. The momentum is given in units of 2n. The solid lines plotted here correspond to the 
simple analytical expressions provided in Eqs. (|139[) - (|143[) . Those plotted in Fig. 5 of Ref.— are very similar yet are obtained 
within the standard formalism of many-body physics by summing up an infinite number of ladder diagrams. Experimental 
points from Ref.—. 



k = [tt, tt] - [?r/2 - h, 3tt/2 - h] - [tt/2, -tt/2] 
= [ki,ki] ; ki = k x = k v e {tt/2, tt) , 
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k 



2 



y y si 



sin(k x /2) , 



= [tt, -tt] - [tt/2 - k x /2, -tt/2 - k x /2] 
- [tt/2 - k x /2, -tt/2 + k x /2] 

= [k x ,o]\ k x e (o,7r), 



(141) 



k 



cos(k y /2) , 



[tt, tt] - [-k y /2, tt - k y /2] - [ky/2, -ky/2] 



[tt, k y ] ; k y € (0, tt) . 



(142) 



= y + ^ s icos(fc,), 

fc = [tt, -tt] - [0, — 7r] - [tt - k x , -tt + k x ] 

= [tt,tt] — [0,7r] — [k y) — k y ] 

= [k x ,Tt-k x ]; k x e(ir/2,TT) 

= [7T-k y ,ky]; k y e(0,Tr/2). (143) 

The theoretical spin excitation spectra (|139|) - (|143|) are plotted in Fig. 1 (solid line) for fi° = 565.6 meV and 
Wg X = 49.6 meV together with the experimental results (circles) for T = 10 K. This gives a Mott-Hubbard gap 
2/i° = 1131.2 meV. The spin-spectrum expressions provided in Eqs. (|139p ~ (|143[) refer to the first Brillouin zone. 
In Fig. 1 we plot them in the second Brillouin zone, alike in Fig. 3 (A) of Ref^ and Fig. 5 of Ref. 41 . An 
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U/At 


A 




^°/2 


< 1 




m 


16* 


MO 


t / 7v w 0.318* 


e L t/-K « 0.865 * 


e L t/ , K w 0.865 * 




0.285 1 


0.673 i 


0.959 1 


> 19 


n [Atf/U 


7r[4t]7C 


[(7/4 - 2*] 



TABLE III: Approximate magnitudes of the energy scales 2Ao, 4Wji, and /i°/2 for different values of U/At. The energy 

parameter A interpolates between fi°/2 w 16* e ~ n V 4t/u for (7/4i < 1 and 4W°j. w tt [4*] 2 /(7 for C//4* > 19 and goes through 
a maximum magnitude max{Ao} pa t/ir at the [7/4* value U/At — uo ~ 1.302 at which jtt°/2 ~ 4W°i. The magnitudes of Ao, 
4Wji, and /i°/2 found here for [7/4* = u* = 1.525 are also given. 

excellent quantitative agreement is reached for these magnitudes of the involved energy scales, which according to the 
complementary results of Ref4^ correspond to U/At pa 1.525 and t pa 295 meV. Moreover, from the use of the relation 
A ~ r s 4W°x of Eq. (|102j) we find that A pa 84 mcV for [7/4* pa = 1.525. The magnitudes of the energy scales 
Ao, 4147°, and fi°/2 found here for [7/4* pa 1.525 are given in Table III together with those for C7/4* < 1, U/it = u , 
and U/At ^> 19. The magnitude of energy parameter Ao interpolates between that of /x°/2 for C7 /4* <C 1 and 4W S ° 1 
for [7/4* 3> 1. It reaches a maximum magnitude approximately given by Ao ~ */vr ~ 0.318* at U/At = uo« 1.302. 

For the intermediate value U /At pa it* = 1.525 € (tto, ui) corresponding to Fig. 1 our investigations find the relation 
W%i pa 0.168*, which for a constant value of * provides a magnitude about twelve times smaller than that found 
by the use of the limiting expression Wjh ~ J ~ Ant 2 /U pa 2.060* for that [J/4* value. Hence for the model on the 
square lattice the usual energy scale J pa Ant 2 /U controls the physics for a smaller [7/4* range than in ID, which 
corresponds to very large U/At ^> 19 values. As discussed in Ref£, the value [7/4* pa 19 is that at which the energy 
scale J pa Ant 2 /U reads J = 0.168* and thus has the magnitude that reaches at U/At pa 1.525. This is why the 
large-[7 /At limiting expressions pa J and Ao pa 4 J are valid for U /At 3> 19. It follows that the intermediate- U /At 
range plays a major role in the physics of the square-lattice quantum liquid studied in this paper and in Refs i 28 i 29 . 

The above results are consistent with the momentum values of the c and si fermions being good quantum numbers. 
Furthermore, they confirm that the predictions of the square-lattice quantum liquid theory concerning the spin 
spectrum at half filling agree both with experiments on the parent compound LCO and results obtained by the 
standard formalism of many-body physics. 

C. Relation to other schemes where the fermions emerge from individual spin-1/2 spins or spinons 

For most previous studies on the large-[7 Hubbard model and * — J model on a square lattice involving for instance 
the slave particle formalis m 4 ' 24 i 25 or Jordan- Wigner transformations 1 ^ the spinless fermions arise from individual 
spin-1/2 spins or spinons. In contrast, within our si fcrmion extended Jordan- Wigner transformation the si fermions 
emerge from spin-neutral two-spinon composite si bond particles. Indeed, here the hard-core bosons refer to spin- 
neutral two-spinon si bond operators. Often in previous related schemes involving spin-liquid mean-field theories, 
which refer in general to individual spin-1/2 spinons or spinless fermions that emerge from them, such objects are 
coupled to a gauge field 4 -" 2 ^. Since in 2D gauge theory is confining, in some of these liquids individual spin-1/2 spinons 
are not observable, as in for instance in the simple dimer state 2 ^. 

Here the original hard-core particles themselves are spin-neutral two-spinon objects. Furthermore, since the spin- 
1/2 spinons have been constructed to inherently referring here to the rotated-electron singly occupied sites, the 
one-rotated-electron-per-site constraint is naturally fulfilled for [7/4* > so that the fluctuations needed to enforce 
such a constraint are for the square-lattice model incorporated in the electron - rotated-electron unitary operator V. 
Related resonating- valence-bond pictures for spin-singlet occupancy configurations of ground states were considered 
long ago 25 ' 66 , yet often problems arise in the construction of energy eigenstates due to bond states being overcomplete 
and non-orthogonal. 

Motivated by previous studies where the spinless fermions arise from individual spin-1/2 spins or spinons, we have 
also performed a Jordan- Wigner transformation for the model both on the square and ID lattices whose spinless 
fermions emerge from the spin-1/2 spinon operators si of Eq. (fTT]) . According to Eqs. (fT3"|) and (JT3J) such operators 
have a hard-core character and the corresponding spinons are well-defined for U/At > 0. One then finds out that for 
ID the momentum values of the emerging spinless fermions do not coincide with the quantum numbers of the exact 
solution. In contrast, the discrete momentum values of the si fermions emerging from spin-neutral two-spinon si 
bond particles correspond to such quantum numbers, as confirmed in Refi^. 

For the model on the square lattice in the one- and two-electron subspace considered in this paper and of interest 
for the further investigations on real materials of Refsi 28 i 29 , our description occupancy configurations of the c fermions 
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and si fermions generate according to the analysis of Rcf. 9 energy cigenstates. This is in contrast to the alternative 
scheme in terms of spinless fermions emerging from the spin-1/2 spinon operators s-. of Eq. (jlip. which generates 
a complete set of states in that subspace, which however are not in general energy eigenstates. The same applies 
to the ID model in the whole Hilbert space. Following the results of Refj£, the momentum values of the c and si 
fermions are good quantum numbers for the model on the square lattice in the one- and two-electron subspace. Our 
scheme uses several approximations to access the shape of the si boundary line, and the c and si fcrmion energy 
dispersions. Our corresponding results are not exact, yet they rely on a general subspace structure and correct 
symmetries, which are compatible with the unknown exact expressions of the quantities under consideration. We 
expect that our approximations provide a good description of the physics for intermediate values of U /At, consistently 
with the results of the previous subsection concerning the half-filling spin spectrum. That spin spectrum is found 
above to lead for the model on the square lattice to analytical expressions for the spin- wave spectrum that agree with 
the corresponding numerical curves obtained by a different controlled approximation involving the sum of an infinite 
set of electronic ladder diagrams^. 

For the model on the square lattice the above alternative scheme leads to a vector potential of the same general 
form as that given in Eq. (|2~4")1 . but with the operator n?. iS i of Eq. (|2li)) replaced by [1/2 + s^ 3 ] and the real-space 
variable fj referring to the spin effective lattice rather than to the si effective lattice. Except that such a scheme 
refers here to the spins of the rotated electrons that singly occupy sites and is valid for U /At > 0, it is identical to the 
spin problem and corresponding Jordan- Wigner transformation considered in Refj^. 

A qualitative difference, which has impact on the physics, is that ai>0 and m = ground state is described here 
by a full si momentum band whereas for the spinless fermions emerging from the spin-1/2 spinon operators si of 

Eq. (fTTj) it is described by a half-filled spin band. This is alike the Wigner- Jordan (WJ) fermions associated with 
the WJ representation of quantum spins with antiferromagnetic interaction of RefJ^, which at zero temperature fill 
exactly half of the total states. In turn, within our c and si fermion description the si momentum band is full for 
the initial x > and m = ground state and addition or removal of one electron leads to the emergence of a single 
si fermion hole in that band. Indeed, following the transformation laws under the electron - rotated-electron unitary 
transformation of the objects whose occupancy configurations generate the energy eigenstates of the square-lattice 
model in the one- and two-electron subspace, excited states with three or a larger odd number of si fermion holes 
have nearly no overlap with one-electron excitations 9 . Furthermore, symmetry selection rules 9 impose that such an 
overlap vanishes exactly for excited states with an even number of si fermion holes. 

The emergence of a single hole in the si band upon creation or annihilation of one electron plays a major role in the 
strong phase-space, exclusion-principle, and energy and momentum conservation restrictions of the si - si fermion 
and si - c fermion scattering. Within the emergence of a single si fermion hole there are obviously no available 
final channels associated with si - si fermion scattering. The above restrictions, together with the isotropic and 
anisotropic character of the c Fermi line and si boundary line, respectively, leads for x > to unusual one-electron 
scattering properties 29 . This includes an anisotropic one-electron inverse lifetime and scattering rate linear in the 
excitation energy lo, in excellent quantitative agreement with that observed in recent angle- resolved photoemission 
studies of the optimally doped high-temperature superconductor LSC O 30 i 31 . Such an inverse lifetime is different from 
that characteristic of an isotropic conventional 2D electron liquid, which is proportional to (u>/Ep) 2 hx{Ep/(jS) where 
Ep denotes the Fermi energy^. 

In turn, for the spinless fermions emerging from the spin-1/2 spinon operators si given in Eq. (jlll) . that their 
momentum band is half filled implies that the ground-state configurations contain as many spinless fermions as spinless 
fermion holes. It follows that in contrast to the si fermions, upon one-electron excitations there are available final 
channels for both spinless fermion - spinless fermion and spinless fermion - c fermion inelastic collisions. In addition, 
since the momenta carried by such spinless fermions are not good quantum numbers, their interactions are not residual 
so that in contrast to the c and si fermion description one cannot use Fermi's golden rule in terms of a collision integral 
as that associated with the c - si fermion interactions to calculate the one-electron lifetime 2 ^. In principle the final 
expressions should be the same independently of the representation. However, the one-electron scattering problem is 
much more complex when expressed in terms of the interactions of the spinless fermions emerging from the spin-1/2 
spinon operators si and so far we could not solve it. 

3 

Finally, only within our description do the si fermion occupancy configurations that generate the spin degrees of 
freedom of x > and m — ground states of the ID model become in the U /At — > oc limit those of the spins of 
the well-known spin-charge factorized wave function introduced independently by Woynarovich 3 ^ and by Ogata and 
Shiba 3 -, respectively. Moreover, only within that description become such ground-state configurations of the model 
on a square lattice in the U/At — > oo limit and within the suitable mean- field approximation (|29[) for the fictitious 
magnetic field B s i of Eq. ([2~Tj) those of a v s \ — 1 full lowest Landau level with TV 2 = N/2 one-sl-fermion degenerate 
states of the 2D QHE. 
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VII. CONCLUDING REMARKS 

The operator description introduced in Refi^ for the Hubbard model on the square lattice in the one- and two- 
electron subspace as defined in that rteference has been constructed to inherently the discrete momentum values of 
the c and si fermions beeing good quantum numbers. Complemcntarily, in this paper we derive approximate results 
and expressions for the si band boundary line, c and si fermion energy dispersions, and corresponding c and si 
fermion velocities of the square-lattice quantum liquid introduced in that reference. 

That the si band remains full alike for for x > and m = ground states alike for x = is related to the degree 
of localization in real space of the si fermion holes, which upon one- and two-electron excitations emerge near the 
si boundary line at momenta pointing in or near the anti-nodal directions. Indeed, the si fermion velocity studied 
in Section V vanishes for si boundary-line momenta pointing in such directions. The corresponding finite-energy 
real-space excitations occur near the pseudogap energy and break translational symmetry locally, due to the square 
si effective lattice spacing being given by a s \ = y/2/(l — x) a rather than by the original square-lattice spacing a. 
In turn, for x > the c band has for the ground state an isotropic c Fermi line associated with a finite c fermion 
velocity so that the low-energy one- and two-electron excitations involving creation of c fermions or c fermion holes 
near the c Fermi line and si fermion holes with momenta pointing in and near the nodal directions have a delocalized 
character. Consistently, in contrast to the isotropic c Fermi line it is found in this paper that the si boundary line is 
rather anisotropic, the corresponding si fermion velocity vanishing and reaching its maximum magnitude for momenta 
pointing in the anti-nodal and nodal directions, respectively. 

Such properties of the c Fermi line and si boundary line are behind the anisotropy of the Fermi velocity, which 
is larger for smaller x. They remain valid for the square-lattice quantum liquid weakly perturbed by the effects of 
3D anisotropy investigated in Refsi 28 i 29 . The results of that reference are fully consistent with the two-gap scenario 
observed in the hole-doped cuprates superconductors^. For instance, there is a direct connection to the recent 
experimental studies of RefA Indeed, the excitations involving coherent virtual-electron pairs whose spinless charge- 
— 2e zero-momentum c fermion pairs are assisted by spin-singlet two-spinon composite si fermions of momenta pointing 
in or near the nodal directions correspond to the superconducting superfluid of delocalized Cooper pairs in momentum 
space observed in such experiments. In turn, the finite-energy excitations involving si fermions whose momenta point 
in and near the anti-nodal directions refer to the locally translational-breaking pseudogap states in real space observed 
in the same experiments. 

Furthermore, the investigations of Refs . 28 ' 29 provide evidence that accounting for the interplay of the electronic 
correlations described by the square-lattice quantum liquid of c and si fermions studied in this paper with the weak 
effects of 3D anisotropy and intrinsic disorder leads to a successful theoretical description of the unusual properties of 
the hole-doped cuprates. Indeed, the corresponding scheme introduced in Ref* 2 - 9 -, leads to quantitative agreement with 
the universal properties of the hole-doped cuprates with superconducting zero-temperature critical hole concentrations 
x c ~ 0.05 and ~ 0.27. For instance, the isotropic character of the c Fermi line and the anisotropy of the si boundary 
line found in this paper combined with the si band being full and having a single hole for x > and m — ground 
states and their one-electron excited states, respectively, leads to unusual one-electron scattering properties controlled 
by the c - si fermion inelastic collisions. The investigations of that reference profit from such a c and si fermion 
description of the one-electron scattering problem and find a strongly anisotropic one-electron inverse lifetime and 
scattering rate in excellent quantitative agreement with that observed in recent angle-resolved photoemission studies 
of the optimally doped high-temperature superconductor LSCO 30 ' 31 . The studies of Ref^S also successfully address 
the normal-state linear- T resistivity in that and other hole doped cuprate o 32 i 33 , neutron resonance mode observed 
in families of hole doped cuprates with critical temperature T c s» 95 Ki, hole-concentration x dependence of several 
physical quantities such as T c and the superfluid density, and the LSCO neutron-scattering low-energy incommensurate 
peaks^. 

In conclusion, the results of Refsi 28 ' 29 confirm that addition to the square-lattice quantum liquid of c and si 
fermions studied in this paper of small 3D anisotropy and intrinsic disorder effects leads to a successful description of 
the unusual properties observed in the hole-doped cuprates. 

Acknowledgments 

I thank Miguel A. N. Araujo, Daniel Arovas, Miguel A. Cazalilla, Karlo Penc, Nuno M. R. Peres, Pedro D. 
Sacramento, and Maria J. Sampaio for discussions and the support of the ESF Science Program INSTANS and grant 
PTDC/FIS/64926/2006. 



56 



Appendix A: The one- and two-electron subspace 

The square-lattice quantum liquid of c and si fermions studied in this paper refers to the one- and two-electron 
subspace of the Hubbard model on the square lattice as defined in Ref.=. Here we provide some basic information 
about such a subspace needed for the studies of this paper. The results summarized in the following are derived and 
discussed in that paper and in Ref»i&. 

For hole concentrations < x < 1 and subspaces spanned by energy eigenstates with constant eigenvalue S c of the 
generator of the global U(l) symmetry of Ref.— there is a vacuum \0 VS }, which remains invariant under the electron 
- rotated-electron unitary transformation. For x = 1 such a vacuum coincides with the electronic vacuum whereas for 
< x < 1 it refers to the m = (1 — x) fully polarized ground state. For the ID model one can derive a closed form 
analytical expression for the critical magnetic field H c for fully polarized ferromagnetism as function of the electronic 
density n = (1 — x) and u = U/At. It is given in Eq. (2) of Ref.— and plotted in Fig. 1 of that reference. Alike 
in ID, it is expected that for the model on the square lattice the magnitude H c of the in-plane critical field for fully 
polarized ferromagnetism is for all values of U /At a decreasing function of the hole concentration x, which vanishes 
for x 1. It is expected as well that for x < 1 it is a decreasing function of U /At, which vanishes for U /At — > oo. 
Such a vacuum reads, 

|0„.) = |0„; <) x %-NZ) x \GS C - 2S C ) , (Al) 

where the ry-spin SU{2) vacuum \Q V \N^) associated with independent +1/2 77-spinons, the spin SU{2) vacuum 

|0 S ; N®) with independent +1/2 spinons, and the c U(l) vacuum \GS C ;2S C ) with N c = 2S C c fermions remain 
invariant under the electron - rotated-electron unitary transformation. For such a fully polarized problem the on-site 
electronic repulsion has no effects so that for m — > (1 — x) and H — > H c the rotated electrons become electrons. In 
that limit the N c — 2S C = N c fermions are the non-interacting spinless fermions, which describe the charge degrees 
of freedom of the N electrons. The c U(l) vacuum \GS C ; 2S C ) can be written as, 

|GS C ;2S C ) = JJ/^ C |GS C ;0). (A2) 

j'=i 

Here \GS C ; 0) refers to the x — 1 and N = electronic vacuum, 

|(W) = |0„; Nq) x |0 S ; 0) x \GS C ; 0) . (A3) 

The state (|A2|) is a particular case of the state of the same form appearing in Eq. (J32J . In turn, the vacuum (|A3j) 
corresponds to a limiting case of the vacuum (|A1|) associated with the values — , Na s = 0> an d 25 c = 0. 

The c fermion, av fermion, and independent 77-spinon and independent spinon description of Ref^ refers to a 
complete set of A N ^ momentum eigenstates | ( &c//4t) = V'^oo) generated from application onto the corresponding 
U/At —> 00 states |$oo) of the electron - rotated-electron unitary operator V* . Such states can be generated from 
their LWS as, 

\$ u/u )= [] ( J _' 1/2 \$ LWS .,v/4t) i C a = S La _ 1/2 ^ + n l[L a + l-l], (A4) 

a— 77, s V a 

where L a _ 1 / 2 is the number of 77-spin-projection —1/2 independent ry-spinons (a = rf) or spin-projection —1/2 

independent spinons (a — s). The corresponding LWS \^LWS-,u/it) = *^\^LWS;oo) an d U/At — > 00 LWS belonging 
to the same V tower read, 

\®LWS;U/At) = [UUUfh,^N a n)}[]lfl c \GS c -,0)} ; j\, ^ av = V^^^V; 4 = V^t c V, (A5) 

a v q' q 

and 

\*lws;oo) = lUUU4'^;N a n)][]jTl c \GS c ;0)], (A6) 

ol v q' q 

respectively. Here T^ c and J 7 !, au are the c fermion and av fermion creation operators for U/At — > oo, respectively. 
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For the ID Hubbard model the states (|A4[) and (|A5[) are both energy and momentum eigenstates. In turn, since 
for the model on the square lattice the set of av translation generators q av in the presence of the fictitious magnetic 
fields B al , considered in Ref.— whose eigenvalues are the components of the av band momenta do not in general 
commute with the Hamiltonian and the corresponding set of av fermion numbers {N av \ are not in general conserved, 
such states are not energy eigenstates and the microscopic momenta carried by the av fermions are not in general 
good quantum numbers. However, such states have been constructed to inherently being S v , $J, S s , S*, S c , and 

momentum eigenstates and each energy eigenstate \^u/it) — ^l^oo) can be expressed as a suitable superposition of 
a well-defined set of \$u/4t) = V |^°o) momentum eigenstates (|A4|) with the same momentum eigenvalue and values 
of S v , Sf, S s , Sg, S c , C v =J2 V v Nr)v, and C s =^2 V vN sv and the same c fermion momentum distribution function 
N c (q)S.. As discussed in this paper, for the model on the square lattice in the one- and two-electron subspace the form 
of the states (|A4[) and (|A5|) simplifies to that provided in Eq. (|32|) and as justified in Refi^ such states are energy 
eigenstates so that the microscopic momenta carried by the c and si fermions are good quantum numbers 

For spin densities to = and m = (1 — x) and U/At — ¥ oo and U/At > 0, respectively, the c fermions are for 
< x < 1 non-interacting and for the model on the square lattice have an energy dispersion given by, 

e ° c (q) = - 2 * [cos(gxi) + cos(gx 2 ) - cos(k Fxi ) - cos(k Fx2 )} . (A7) 

Here kp Xl and kp X2 are the Cartesian coordinates of the Fermi momentum hp, 

hp = kp{<j))e ( t,; ; = arctan ( ) ; e (0,2tt) . (A8) 



For such values of m and U/At the corresponding Fermi line encloses a momentum area, 

^[fc F (0)] 2 = (l + .x)47r 2 . (A9) 

It is twice the area enclosed by the Fermi line at to = and U/At = 0. 

Two energy scales that play a key role in the square-lattice quantum liquid are the chemical potential fi and the 
maximum si fermion spinon-pairing energy 2|A| of Eq. (|104[) . which in the limit x — > equals the energy scale 2Ao 
whose magnitudes are provided in Table III. Within the LWS representation, the convention that /i has the same sign 
as the hole concentration x is used. For < x < 1 and to = the chemical potential is an increasing function of the 
hole concentration x whose values belong to the range, 



fj° < n{x) < n 1 ; < x < 1 , to = 0, (A10) 



where /i 1 = lim,^! /i. /i 1 reads, 



/i 1 = [U/2 + 2Dt] ; D = l,2. (All) 
Moreover, /x° has the following approximate limiting behaviors^, 

' e' 2 ^) 1/D , U/4t<l; /i° « [U/2 - 2Dt] , U/Af»l, D=l,2, (A12) 

so that fjP — y as U/At — y whereas fjP — > oo for U/At ^ 1 for both the model on the ID and square lattices. As 
discussed in previous sections, 2yP is the half-filling Mott-Hubbard gap. 

Another energy scale that plays an important role is the energy e c for addition onto ai>0 and m = ground 
state of one c fermion and the corresponding energy energy —e c for removal from that state of one c fermion, which 
obey the inequalities, 

< e c < [ADt - WP] ; < -e c < W? ; D = l,2, (A13) 

respectively. In turn, the energy —e s \ for addition to a x > and to = ground state of one si fermion hole obeys 
the inequality, 

0< -e sl <max{W s i,(£»-l)|A|}, £ = 1,2. (A14) 

In such inequalities W]! and W c ' 1 = [ADt — W^} are the ground-state c fermion and c fermion-hole energy-dispersion 
bandwidths, respectively, = [ADt-Wf>] e (0,ADt) increases monotonously upon increasing the hole concentration 
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within the range x € (0, 1), W s \ is the si fermion auxiliary energy-dispersion bandwidth of Eq. (|104l) . and the energy 
parameter |A| corresponds for the model on the square lattice to the maximum si fermion pairing energy per spinon 
of Eq. (|107j) . The energy = \im x ^.Q W s i has the following approximate limiting behaviors^, 

W° sl = 2Dt; U/4± = 0, 

« 2D%t 2 /U; U/At > 1 . (A15) 

For the model on the square lattice the large-[//4i expression is expected to be a good approximation for very large 
values U /At ^S> 19, as justified in this paper and Refi^. 

The one- and two-electron subspace is spanned by a given initial x > and m = ground state and a well-defined 
set of its excited energy eigenstates of excitation energy to < 2/x for x > and lu < /i° for x = 0. Such states have no 
— 1/2 77-spinons, 2z/-7/-spinon composite r\v fermions, and 2z/-spinon composite sv' fermions with v 1 > 3 spinon pairs 
so that N vv = and N su > = for v 1 > 1, whereas according to the studies of Ref£ the numbers of independent ±1/2 
spinons and that of s2 fermions are restricted to the following ranges, 

L s ,±l/ 2 =0,1; N s2 =0, Af = l; 2S S + 2N s2 = 0, 1, 2 , M = 2. (A16) 

Here J\f = 1, 2 refers to the AA-electron operators whose application onto the ground state leads to the excited states 
under consideration. Furthermore, the numbers of c and si fermions read N c = N = (1 — x) and N s i = 
[N/2 - 2N s2 - S.] = (1 - as) N? /2 - [2N s2 + S s ], respectively. 

From analysis of the ground-state occupancies found in Refi^, for the one- and two-electron subspace there occur 
the fourteen classes of elementary excited states relative to an initial x > and m — ground state whose numbers 
and number deviations are given in Table IIVI The table provides the values of the deviations SN^ and numbers 
of such excited states, corresponding electron number deviations SN^- and 5N±, and independent-spinon numbers 
L s .+i/2 an d L s> -1/2 and s2 fermion numbers N s2 restricted to the value ranges of that subspace. The spin S s and 
deviations SS C , SN s i = [SS C — S s — 2N s2 ], and 5N asl = [5S C + S s ] of each excitation are also given. 

In reference^ it is found from analysis of the transformation laws under the electron - rotated-electron unitary 
transformation of the objects whose occupancy configurations generate the energy eigenstates of the Hubbard model 
that nearly the whole one- and two-electron spectral weight is contained in the one- and two-electron subspace whose 
numbers are given in Table IV. If in addition we restrict our considerations to the LWS-subspace of the one- and 
two-electron subspace, then L s ,_i/2 = in Eq. (|A16I) . whereas the values L s ,+i/2 = 0, 1 for 7V S 2 = and Af = 1 
remain valid and in 2^ + 2N s2 = 0, 1, 2 one has that 2S S = L s +1 / 2 for J\f = 2 in that equation. 

For the one- and two-electron subspace of the Hamiltonian ((TJ) only the c fermions and si fermions play an active role. 
The numbers and N® ± of sites of the c and si effective lattices, respectively, equal the number of corresponding 
c and si band discrete momentum values. The number , that of si fermions N s i, and that of both unoccupied 
sites of the si effective lattice and si band holes N' s \ are for the one- and two-electron subspace given by, 

N^ x = N sl + N% = N°J2 +S S = S C + S S ; N? m = (1 - x) N? , 

N,i = N° /2-S s + 2N s2 =S C -S S - 2N s2 ; = 2S S + 2N s2 = 0, 1, 2 , (A17) 

respectively, where = (1 — x)N® is the number of sites of the related spin effective lattice. 

It is found in RefJ^ that for the one- and two-electron subspace and Nj^ ^> 1, the operators g^.si (and gt sl ) 
related to the annihilation and creation of local sl fermion operators through Eq. (|16|) . which annihilate (and create) 
a sl bond particle at a site of the sl effective lattice of real-space coordinate fj, have the following general form both 
for the model on the ID and square lattices, 

N.i/2D-1 D 

9? j ,sl= X! h g a r j ,sl,g\ gl jtS i = (g? s , sl) Ofj ,sl,g = ^ ^ &F 3 ,+?° , ,sl,d,l,g ■ ( A18 ) 

3=0 d=ll=±l 

Hence the expression of the operator gl sl involves the operators, 

4 jlS l,g = ( a ri,sl,g) ] b l, s l,d,Ug = ( b W,d,l,g) j ■ ( Al9 ) 

The operators at. sl g and ap j<s i,g create and annihilate, respectively, a superposition of 2D = 2,4 two-site bonds of 

the type studied in RefJ^ and bt sl d ; and b?^ s \^d,i,g are two-site one-bond operators whose expression is given in 
that reference. 
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TABLE IV: The deviations SN^ = —25S C and numbers N al = [2S S + 2iV S 2] for the fourteen classes of one- and two-electron 
excited states of the x > and m — ground state that span the one- and two-electron subspace as defined in Ref. ? , 
corresponding electron number deviations SN-f and 5N±, and independent-spinon numbers L St +1/2 and L Sj _i/2 and s2 fermion 
numbers N S 2 restricted to value ranges of that subspace. The spin S s and deviations SS C , SN a i — [SS C — S a — 2iV S 2], and 
8Na sl = [SSc + Ss] of each excitation are also provided. 



Appendix B: Confirmation of the gauge structure of the underlying si effective lattice 

The si bond-particle description of Ref*i2 involves a change of gauge structur e) 10 ! 24 so that the real-space coordinates 
of the sites of the si effective lattice correspond to one of the two sub-lattices of the square spin effective lattice. Here 
it is confirmed that for the limit 1 that the description used in our studies refers to and for hole concentrations 

x such that the density n = (1 — x) is finite the two choices of si effective lattice lead to the same description. 

The initial ground states of the one- and two-electron subspace are x > 0, m = 0, and = states. Here 
Ng X = [N?? — N s i] is the number of unoccupied sites of the si effective lattice. Its expression is for the one- and two- 
electron subspace given in Eq. (IA17[) of Appendix A along with that of the number N s i of si fermions. We emphasize 
that for the present limit ^> 1 the transformation of Eq. (|25l) leads precisely to the same operators ft sl 

and fq jtS i independently on whether the real-space coordinates ry of the corresponding operators /1 ; sl and fp.,, a i, 

respectively, are those of any of the two sub-lattices of the spin effective lattice. Indeed, the two choices of real-space 
variables of the sl effective lattice considered in Ref.— are for = states related both by the transformations 
ry — > Tji +a s i e Xl and ry — ► ry +a s i e X2 . If in the expressions of Eq. (l25l we replace ry by ry +a s i e Xl or ry +a s i e X2 
and transform the summations in integrals, for the periodic boundary conditions considered in Ref.™ we obtain within 
the limit N® sl — > 00 exactly the same operators ft sl and fq j>8 i- 

This confirms that the change of gauge structur e) 10 ! 24 used in the construction of the sl effective lattice, which 
consists in choosing one of the two sub- lattices of the square spin effective lattice to play the role of sl effective 
lattice, does not affect the momentum values of the sl fermions: One reaches the same momentum values and thus 
the same physics for the two possible choices of sl effective-lattice real-space coordinates. 

That change of gauge structure was made in Refi^ - for the N^i = configuration state that generates the spin 
degrees of freedom of the x > and m = ground states. As discussed in that paper, all A^-finite configuration 
states associated with the excited states generated by application of one- and two-electron operators onto such ground 
states can be straightforwardly constructed from the Nj^ = configuration state. The point is that the momentum 
sl fermion operators obtained from the transformation of Eq. (|25p are precisely the same for any finite number Njh 
of sl band holes, independently on whether the corresponding A^\-finite configuration states are generated from the 
= configuration state when described by occupancy of one or the other sublattice sites. This confirms that for 
the limit N®^ ^> 1 the two alternative choices of real-space coordinates of the sl effective lattice refer indeed to a 
gauge structure, as stated in Refill. 

That such a change of gauge structure leads to the same sl fermions independently of the choice of the real-space 
coordinates of the sl effective lattice is a result that applies both to the model on the square and ID lattices. 

Appendix C: The x dependence of the elementary function e s i(g) for 0<i<1 and U/4t > «o 

The goal of this Appendix is the evaluation up to first order in x of the elementary function e s \{q) that for < x <C 1 
and U/4t > uo leads to the expressions qg sl ~ v / 2[i'/2](l — x) and ~ vr — \J x 2ir for the absolute values of the 
auxiliary nodal and anti- nodal momenta of Eq. (|99|) . respectively, energy scale 2| A| = 2Aq(1 — x/x®) and parameter 
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2r s /n given in Eq. (|104p . Specifically, such expressions are reached provided that e s \{q) 



[W^/2] cosq 



remains a good approximation for < x -C 1 and U/4t> uq except for q values near both zero and q^i and is given 
by, 



e s i(g) 



w°, 

v v sl 

2 

v y sl 

2 



cosqu ; q € (0,qu) , 
cosg; qE (qu,n~ q v ) . 



M 2 



qe {Tr~qu,q B sl) 



qu 



(l-r s ); 0<x<l, U/4t>u . 



(CI) 



The above expressions are then consistently obtained for < x <C 1 and U /At > uq by expanding up to first order in 
x the solution of the equation, 



0,||/.*AjV\_ W sl 
e sl \1Bs1 ) — ^~ 



COS qu + COS I 7T - \ l [TT - gg^ 



^2 j. ^2 



= ; < x < 1 , E//4f > m . 



(C2) 



In turn, use in Eq. (|101[) of the expressions provided in Eq. (ICip for the sl elementary function leads to the 
following x dependence for the energy scale 2|A|, 



2|A| = A 



cos qu - cos 7T - J [tt - ggfj 2 + [gc/] 



, < x < 1 , f//4t > u . 



(C3) 



Taking into account that [tt — ?s^i] 2 ~ x2tt and expanding within the < i « 1 limit the right-hand side of Eq. 
(|C3j) up to first order in x then leads to the x dependence 2|A| = (1 — x/x®)2A for U/At> uq provided in Eq. (|104j) 
where x® = 2r s /n. 

The expression given in Eq. (|C1[) has only physical meaning up to first order in x and up to that order the 
function [e s i(0) + e s i(g^)]/l / l /r ° 1 = —xtt/2 is independent of U/it so that the solution of Eq. (|C2p is given by 
[tt - q£%] 2 x2tt. In turn, [e s i(0) - e,!^^)]/^ depends on U/4t and through Eq. JC3| is behind the U/4t 
dependence of the parameter x° c — 2r s /ir of Eq. (|104[) for the range U/At > uq, which the present derivation refers to. 
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